Academic literature on the topic 'Fluid structure interaction, mechanical heart valves, vortex formation'

The work is focused on calculating hemodynamically negative effects of a flow through bileaflet mechanical heart valves (BMHV). Open-source FOAM-extend and cfMesh libraries were used for numerical simulation, the leaflet movement was solved as a fluid-structure interaction. A real model of the Sorin Bicarbon heart valve was employed as the default geometry for the following shape improvement. The unsteady boundary conditions correspond to physiological data of a cardiac cycle. It is shown how the modification of the shape of the original valve geometry positively affected the size of backflow areas. Based on numerical results, a significant reduction of shear stress magnitude is shown. The outcome of a direct numerical simulation (DNS) of transient flow was compared with results of low-Reynolds URANS model k-ω SST. Despite the limits of the two-dimensional solution and Newtonian fluid model, the suitability of models frequently used in literature was reviewed. Use of URANS models can suppress the formation of some relevant vortex structures which may affect the BMHV’s dynamics. The results of this analysis can find use in optimizing the design of the mechanical valve that would cause less damage to the blood cells and lower risk of thrombus formation.

We examine the development of a thermal plume originating from a localized heat source using direct numerical simulation. The Reynolds number of the plume, based on source diameter and the characteristic buoyancy velocity, is chosen to be 7700, which is sufficiently large so that the flow turns to a fully turbulent state. A highly resolved grid of 622 million points is used to capture the entire range of turbulent scales in the plume. Here at the source, only heat has been added with no mass or momentum addition and accordingly the vertical evolution of the mass, momentum and buoyancy fluxes computed from the simulation have been verified to follow those of a pure thermal plume. The computed vertical evolution of the time-averaged centreline velocity and temperature are in good agreement with available experimental measurements. Investigation of the time evolution of the plume shows periodic formation of vortex ring structure surrounding the main ascending column of hot fluid. The vortex ring forms very close to the heat source and even at formation it is three-dimensional. The vortex ring ascends with the plume and at an elevation of about two diameters it strongly interacts with and destabilizes the central column and subsequently a complex turbulent flow arises. Thus, relatively laminar, transitional and fully turbulent regimes of the plume evolution can be identified. In the fully turbulent regime, complex three-dimensional hairpin-like vortex structures are observed; but vestiges of the coherent vortex rolls that form close to the source can be observed in the turbulent statistics. It is shown that local entrainment consists of contraction and expulsion phases. Such instantaneous mechanisms drive the entrainment process, and the instantaneous entrainment coefficient shows large variation in both time and space with local values up to three times higher than the average entrainment level. Such findings support the view that entrainment mechanisms in plumes should be considered from an unsteady point of view. Movies are available with the online version of the paper.

Prosthetic heart valves are commonly used as a treatment for aortic valve deficiencies. The performance of these prosthetic valves should be in accordance with the natural heart valve with respect to opening and closing, blood flow, and vortex formation. These performance parameters depend on the design of leaflets and overall geometrical parameters of the valve. To better understand the effects of leaflet design on the performance of the valve, we have carried out fully coupled fluid–structure interaction analyses of opening and closing of prosthetic heart valves with various leaflet designs. Maximum stress, valve opening, and flow stream pattern are obtained for different valve designs and used to assess the performance of the valves. The results show that the stress and the valve opening depend on the curvature and the inclination of the leaflets. A 3D model is designed based on the obtained results, and a full FSI analysis is performed to assess its performance. The results show that the presented design gives better values for valve opening area and leaflet stresses than that in the published data.

Abstract Bileaflet mechanical heart valves (BMHV) are widely implanted to replace diseased heart valves. Despite many improvements in design, these valves still suffer from various complications, such as valve dysfunction, tissue overgrowth, hemolysis, and thromboembolism. Thrombosis and thromboembolism are believed to be initiated by platelet activation due to contact with foreign surfaces and nonphysiological flow patterns. The implantation of the valve causes nonphysiological patterns of vortex shedding behind the leaflets. This study signifies the importance of vorticity in platelet activation and aggregation in BMHV implants. A two-phase model with the first Eulerian phase for blood and the second discrete phase for platelets is used here. The generalized cross model of viscosity has been used to simulate the non-Newtonian viscosity of blood. A fluid–structure-interaction model has been used to simulate the motion of leaflets. This study has also estimated the platelet activation state (PAS), which is the mathematical estimation of the degree of activation of platelets due to flow-induced shear stresses that cause thrombus formation. The regions in the fluid domain with a higher vorticity field have been found to contain platelets with relatively higher PAS than regions with relatively lower vorticity fields. Also, this study has quantitatively reported the effect of vorticity on platelet aggregation. Platelet densities in fluid areas with higher vorticity are higher than densities in fluid areas with lower vorticity, indicating that highly activated platelets aggregate in areas with stronger vorticity.

The fluid dynamics during valve closure resulting in high shear flows and large residence times of particles has been implicated in platelet activation and thrombus formation in mechanical heart valves. Our previous studies with bileaflet valves have shown that large shear stresses induced in the gap between the leaflet edge and valve housing results in relatively high platelet activation levels, whereas flow between the leaflets results in shed vortices not conducive to platelet damage. In this study we compare the result of closing dynamics of a tilting disk valve with that of a bileaflet valve. The two-dimensional fluid-structure interaction analysis of a tilting disk valve closure mechanics is performed with a fixed grid Cartesian mesh flow solver with local mesh refinement, and a Lagrangian particle dynamic analysis for computation of potential for platelet activation. Throughout the simulation the flow remains in the laminar regime, and the flow through the gap width is marked by the development of a shear layer, which separates from the leaflet downstream of the valve. Zones of recirculation are observed in the gap between the leaflet edge and valve housing on the major orifice region of the tilting disk valve and are seen to be migrating toward the minor orifice region. Jet flow is observed at the minor orifice region and a vortex is formed, which sheds in the direction of fluid motion, as observed in experiments using PIV measurements. The activation parameter computed for the tilting disk valve at the time of closure was found to be 2.7 times greater than that of the bileaflet mechanical valve and was found to be in the vicinity of the minor orifice region, mainly due to the migration of vortical structures from the major to the minor orifice region during the leaflet rebound of the closing phase.

2009/2010
Fluid structure interaction (FSI) is one of fundamental phenomena encountered everywhere in nature or in industrial systems as well as one of the most studied and the most challenging topics in the fluid mechanics. Its research presents the core objective of this dissertation, along with experimental study of artificial heart devices. Better understanding of FSI could turn the still unexploited phenomenon into a powerful tool for resolving wealthy of multi-physics problems. Recently computational fluid dynamics community has been putting enormous efforts to uncover, make clear and answer yet numerous issues related to this developing topic. In addition, the FSI is often followed by the vortex formation, one more phenomena that could be both powerful driving force as well as distracting, disturbing occurrence. Consequently, this dissertation will begin with addressing some open issues related to the fluid-structure interaction associated with the simple system made of movable rigid leaflet and an unsteady viscous fluid flow. Such two-dimensional model, even if it appears extremely simple, is able to produce fairly rich flow features which deserve careful analytical and accurate numerical solution. Thus, we have performed a significant number of numerical experiments with the objective to uncover the role of the structure inertia in the overall behavior of the fluid-leaflet system, under the different flow recurrences. For that purpose, we have constructed a strong-coupling code and resolved the fluid and structure dynamics simultaneously, paying particular care of solution accuracy around the moving boundary. The complex problem of large fluid deformation in response to the rapid structure movements has been resolved by the time-dependent conformal mapping, exclusively developed for this specific physical arrangement. The numerical findings, even if theoretical in nature, allowed for the classification and characterization of body’s and fluid dynamics in functionality of different structure inertia and Strouhal numbers, which have been used as free parameters in all numerical experiments. The study is completed by a brief analysis of the more realistic system of actual prosthetic heart valves. Besides many problems that follow the performance of mechanical heart valve prosthesis, the complications related to the complex blood-leaflet interaction are a key factor. The intraventricular flow is characterized by large vortical structures, without significant turbulence, in a smooth circulatory pattern that, in presence of pathological conditions or mechanical devices, could be disturbed. Thus, among the criteria for the assessment of mitral valve functionality and mechanical valve design are the proper vortical features inside the left ventricle. Until nowadays the standard mechanical valves, designed originally for the aortic replacement and without exceptions symmetrical, have never satisfied the regularity of natural vortical dynamics. Thus, we have been motivated to investigate the flow features downstream of asymmetrical prototypes, exclusively designed for the mitral replacement with attempt to better mimic the natural intraventricular flow. Experimental outcomes allowed for preliminary conclusions that the break of symmetry in the novel prosthesis creates the asymmetrical vortical flow in the left ventricle, which is more similar to the natural one, although the concept introduced by this prototype has to undergo deeper testing and careful improvements before querying in the real hearts.
XXIII Ciclo
1982

Surgical replacement, in the incidence of severely diseased heart valve, is vital in order to restore the normal heart function. Every year around 280,000 valve replacements occur around the world, half of them are bileaflet mechanical heart valves (BMHVs). Despite the remarkable improvement in valve design resulting in minimizing prosthetic valve complications (thromboembolic events or pannus formation), these complications are still possible with BMHV Implantation. As a consequence, an obstruction in one or both MHV leaflets could happen and threaten the patient life. In the present study, an obstructed bileaflet MHV with different percentages of malfunction was simulated assuming 3-D fluid structure interaction (FSI) adapting k-w turbulence as a robust model for the transitional flow using 2.5 million elements and creating a realistic aortic root model with three sinuses. Velocity contours for different percentages of malfunction were compared mainly at B-datum plane and the perpendicular plane to the B-Datum. Also, the development of coherent structures was investigated. Clinically, the maximum pressure gradients were estimated by mimicking the Echo Doppler assumptions (using the simplified Bernoulli equation).

Altered haemodynamics are implicated in the blood cells damage that leads to thromboembolic complications in presence of prosthetic cardiovascular devices, with platelet activation being the underlying mechanism for cardioemboli formation in blood flow past mechanical heart valves (MHVs). Platelet activation can be initiated and maintained by flow patterns arising from blood flowing through the MHV, and can lead to an enhancement in the aggregation of platelets, increasing the risk for thromboemboli formation. Hellums and colleagues compiled numerous experimental results to depict a locus of incipient shear related platelet activation on a shear stress – exposure time plane, commonly used as a standard for platelet activation threshold [1]. However, platelet activation and aggregation is significantly greater under pulsatile or dynamic condition relative to exposure to constant shear stress [2]. Previous studies do not allow to determine the relationship existing between the measured effect — the activation of a platelet, and the cause — the time-varying mechanical loading, and the time of exposure to it as might be expected in vivo when blood flows through the valve. The optimization of the thrombogenic performance of MHVs could be facilitated by formulating a robust numerical methodology with predictive capabilities of flow-induced platelet activation. To achieve this objective, it is essential (i) to quantify the link between realistic valve induced haemodynamics and platelet activation, and (ii) to integrate theoretical, numerical, and experimental approaches that allow for the estimation of the thrombogenic risk associated with a specific geometry and/or working conditions of the implantable device. In this work, a comprehensive analysis of the Lagrangian systolic dynamics of platelet trajectories and their shear histories in the flow through a bileaflet MHV is presented. This study uses information extracted from the numerical simulations performed to resolve the flow field through a realistic model of MHV by means of an experimentally validated fluid-structure interaction approach [3]. The potency of the device to mechanically induce activation/damage of platelets is evaluated using a Lagrangian-based blood damage cumulative model recently identified using in vitro platelet activity measurements [4,5].

Abstract Ventricular assist devices (VAD) are designed to provide circulatory support to patients suffering from advanced-stage heart failure. While not pulsatile, the advantages of continuous axial flow VADs include a compact size and low mechanical failure rate. However, being compact, they operate at high speed, resulting in adverse effects, such as hemolysis caused by high flow shear and thrombosis formation in stagnant regions are common and threaten the successful use of the device. While state-of-the-art computational fluid dynamics (CFD) is widely used in designing these devices, detailed high-resolution experimental measurements of the flow within them are not readily available in the literature. Such experimental data is crucial for understanding the flow inside the VAD and its interaction with blood cells as well as for validating the CFD predictions. The present study investigates the flow inside a 1:1 exact replica of a VAD – ReliantHeart HeartAssist5®. This 12mm diameter device consists of an inlet guide vane (IGV), a rotor and a stator. However, unlike in the real machine, the rotor is driven by a thin shaft that penetrates through the center of the IGV. Refractive index-matching is used to facilitate optical measurements. Hence, all the blades and housing of the pump are made of transparent acrylic. The working fluid is a mixture of water, sodium iodide and glycerin, which matches the refractive index of acrylic, and the kinematic viscosity of blood. Performance tests have been carried out at speeds ranging from 7000 to 9000 RPM. Trends of the results are consistent with those of the actual machine. While scaled data for the pressure rise across the rotor at different speeds collapse, the total head rise across the entire machine does not. High-resolution 2D PIV measurements with vector spacing of 30μm have been conducted in meridional planes of the rotor passage at several blade orientations. They have been performed at 8000RPM and flowrate of 4.5L/min, consistent with physiological requirements. They show that near the rotor front end, the flow in the tip region is dominated by the tip leakage vortex (TLV), associated blockage effects, and very high turbulence level. The upstream end of this domain remains aligned with the leading edge, implying that blades persistently dissect the remnants of a previous TLVs. Interaction of the hub boundary layer with the blades also generates secondary flows. Downstream of the helical section of the rotor, the axial flow is high near the hub and low in the outer part of the passage due to leakage flow-related blockage. Hub boundary layer separation occurs upstream of the stator, generating a large circumferential vortex that occupies nearly half of the rotor span. The complexity of the flow structure and turbulence in this machine would be a challenge to model.