A numerical simulation of mechanical heart valve closure fluid dynamics

Research on fatigue optimization simulation of polymeric heart valve based on the iterative sub-regional thickened method

Prosthetic polymeric heart valves (PHVs) have the potential to overcome the inherent material and design limitations of traditional valves in the treatment of valvular heart disease; however, their durability remains limited. Optimal design of the valve structure is necessary to improve their durability. This study aimed to enhance the fatigue resistance of PHVs by improving the stress distribution. Iterative subregional thickening of the leaflets was used, and the mechanical stress distribution and hemodynamics of these polymeric tri-leaflet valves were characterized using a fluid-structure interaction approach. Subregional thickening led to a reduction in stress concentration on the leaflet, with the effective orifice area still meeting ISO 5840-3 and the regurgitant volume achieving a similar value to those in previous studies. The maximum stress in the final iteration was reduced by 28% compared with that of the prototype. The proposed method shows potential for analyzing the stress distribution and hemodynamic performance of subregional thickened valves and can further improve the durability of PHVs.

Multiscale model for blood flow after a bileaflet artificial aortic valve implantation

Cardiovascular diseases are the leading cause of mortality in the world, mainly due to atherosclerosis and its consequences. The article presents the numerical model of the blood flow through artificial aortic valve. The overset mesh approach was applied to simulate the valve leaflets motion and to realize the moving mesh, in the aortic arch and the main branches of cardiovascular system. To capture the cardiac system's response and the effect of vessel compliance on the outlet pressure, the lumped parameter model has been also included within the solution procedure. Three different turbulence modeling approaches were used and compared - the laminar, k-ϵ and k-ω model. The simulation results were also compared with the model excluding the moving valve geometry and the importance of the lumped parameter model for the outlet boundary condition was analyzed. Proposed numerical model and protocol was found as suitable for performing the virtual operations on the real patient vasculature geometry. The time-efficient turbulence model and overall solving procedure allows to support the clinicians in making decisions about the patient treatment and to predict the results of the future surgery.

Development of an adaptive pulmonary simulator for in vitro analysis of patient populations and patient-specific data

BACKGROUND AND OBJECTIVE

Patient-specific modeling (PSM) is gaining more attention from researchers due to its ability to potentially improve diagnostic capabilities, guide the design of intervention procedures, and optimize clinical management by predicting the outcome of a particular treatment and/or surgical intervention. Due to the hemodynamic diversity of specific patients, an adaptive pulmonary simulator (PS) would be essential for analyzing the possible impact of external factors on the safety, performance, and reliability of a cardiac assist device within a mock circulatory system (MCS). In order to accurately and precisely replicate the conditions within the pulmonary system, a PS should not only account for the ability of the pulmonary system to supply blood flow at specific pressures, but similarly consider systemic outflow dynamics. This would provide an accurate pressure and flow rate return supply back into the left ventricular section of the MCS (i.e. the initial conditions of the left heart).

METHODS

Employing an embedded Windkessel model, a control system model was developed utilizing MathWorks' Simulink® Simscape™. Following a verification and validation (V&V) analysis approach, a PI-controlled closed-loop hydraulic system was developed using Simscape™. This physical system modeling tool was used to (1) develop and control the in silico system during verification studies and (2) simulate pulmonary performance for validation of this control architecture.

RESULTS

The pulmonary Windkessel model developed is capable of generating the left atrial pressure (LAP) waveform from given pulmonary factors, aortic conditions, and systemic variables. Verification of the adaptive PS's performance and validation of this control architecture support this modeling methodology as an effective means of reproducing pulmonary pressure waveforms and systemic outflow conditions, unique to a particular patient. Adult and geriatric with and without Heart Failure and a Normal Ejection Fraction (HFNEF) are presented.

CONCLUSIONS

The adaptability of this modelling approach allows for the simulation of pulmonary conditions without the limitations of a dedicated hardware platform for use in in vitro investigations.

Computational methods for the aortic heart valve and its replacements

INTRODUCTION

Replacement with a prosthetic device remains a major treatment option for the patients suffering from heart valve disease, with prevalence growing resulting from an ageing population. While the most popular replacement heart valve continues to be the bioprosthetic heart valve (BHV), its durability remains limited. There is thus a continued need to develop a general understanding of the underlying mechanisms limiting BHV durability to facilitate development of a more durable prosthesis. In this regard, computational models can play a pivotal role as they can evaluate our understanding of the underlying mechanisms and be used to optimize designs that may not always be intuitive. Areas covered: This review covers recent progress in computational models for the simulation of BHV, with a focus on aortic valve (AV) replacement. Recent contributions in valve geometry, leaflet material models, novel methods for numerical simulation, and applications to BHV optimization are discussed. This information should serve not only to infer reliable and dependable BHV function, but also to establish guidelines and insight for the design of future prosthetic valves by analyzing the influence of design, hemodynamics and tissue mechanics. Expert commentary: The paradigm of predictive modeling of heart valve prosthesis are becoming a reality which can simultaneously improve clinical outcomes and reduce costs. It can also lead to patient-specific valve design.

Effects of Leaflet Design on Transvalvular Gradients of Bioprosthetic Heart Valves

Bioprosthetic aortic valves (BAVs) are becoming the prostheses of choice in heart valve replacement. The objective of this paper is to assess the effects of leaflet geometry on the mechanics and hemodynamics of BAVs in a fluid structure interaction model. The curvature and angle of leaflets were varied in 10 case studies whereby the following design parameters were altered: a circular arch, a line, and a parabola for the radial curvature, and a circular arch, a spline, and a parabola for the circumferential curvature. Six different leaflet angles (representative of the inclination of the leaflets toward the surrounding aortic wall) were analyzed. The 3-dimensional geometry of the models were created using SolidWorks, Pointwise was used for meshing, and Comsol Multiphysics was used for implicit finite element calculations. Realistic loading was enforced by considering the time-dependent strongly-coupled interaction between blood flow and leaflets. Higher mean pressure gradients as well as von Mises stresses were obtained with a parabolic or circular curvature for radial curvature or a parabolic or spline curvature for the circumferential curvature. A smaller leaflet angle was associated with a lower pressure gradient, and, a lower von Mises stress. The leaflet curvature and angle noticeably affected the speed of valve opening, and closing. When a parabola was used for circumferential or radial curvature, leaflets displacements were asymmetric, and they opened and closed more slowly. A circular circumferential leaflet curvature, a linear leaflet radial curvature, and leaflet inclination toward the surrounding aortic wall were associated with superior BAVs mechanics.

Mathematical Modeling and Simulation of an Occlusion Device in a Blood Vessel

An occlusion device is placed in an abnormal opening of the heart or its surrounding vessels to regain normal blood flow. There are various occlusion devices available for treatment of various congenital heart defects like PDA, ASD, etc. However, they have limitations like residual shunting, erosion of tissue, displacement and breakage of device, thrombus formation and sudden death. To improve efficiency and to reduce failure of occlusion devices, it is important to simulate blood flow through defect before and after placement of device. It is also important to evaluate stresses and forces exerted by blood flow on device and by the device on the vessel wall. Contact friction between device and vessel wall plays a crucial role in anchoring the device. The objective is to develop a framework to determine conditions to restrict dislocation of device in terms of contact friction. Typical occlusion devices are porous initially and later due to thrombogenesis, their porosity reduces until eventually it acts as a natural permanent plug. Thus, a porous sponge is a good model for an occlusion device. The mathematical model developed here is for differential pressure causing incipient movement of device, and minimum value of contact friction for restricting movement of the device for two shapes, cylindrical and conical, in uncompressed as well as pre compressed forms. The model for differential pressure is fitted by conducting physical experiment with sponge. Mathematically, porosity is modeled using viscous resistance and inertial resistance which are calculated by experiment and simulation with ANSYS. We perform computer experiments (simulations) on a cylindrical device in a cylindrical vessel and on a conical device in a tapered vessel to determine the differential pressure across the device and hence contact friction with varied porosity under boundary conditions as in body. The contact friction required to retain device is lesser in case of conical device compared to cylindrical device. As compression of device increases, friction require to retain it decreases. Hence, lesser porosity results in larger differential pressure and lesser compression which will eventually need higher friction values to retain the device.

Decreased Hemolysis and Improved Hemodynamic Performance of Synchronized Bileaflet Mechanical Valve

PURPOSE

This study compared the performance of a newly patented synchronized valve (SV) with that of a commercially available (CAV) bileaflet mechanical heart valve.

DESCRIPTION

A high-speed camera was used to record the leaflet kinematics of the SV vs the CAV along the flow channel. Transvalvular energy loss, effective orifice area, and hemolysis ratios were obtained using a mock circulatory system at two fixed pulse rates and at various cardiac outputs with a fixed aortic pressure.

EVALUATION

The rotational radius and inertia of the SV was lower than that of the CAV during valve closure. For heart rates and at cardiac outputs of 7, 5, and 4 L/min, the ratio of total energy loss to effective energy of the SV was significantly less than the CAV, whereas the effective orifice area of the SV was significantly larger than that of CAV. The hemolysis ratio after 4 hours was significantly higher in the CAV than in the SV for both pulse rates.

CONCLUSIONS

The synchronized leaflet motion mitigated leaflet rebound and regurgitation during valve closure, which could decrease energy loss, increase the effective orifice area, and reduce hemolysis.

FLUID DYNAMIC CHARACTERISTICS OF SYSTOLIC BLOOD FLOW OF THE LEFT VENTRICLE

Ejection of blood from the left ventricle to the aorta is studied using two-dimensional Navier–Stokes equations, the work-energy equation and the magnetic resonance imaging of a normal ventricular motion. Vortex shedding in the sinuses of Valsalva is dominated by the aortic jet, flow acceleration and valve motion. Momentums produced by ventricular contraction are in concert with vortices in the ventricle for blood ejection. Shear stresses and net pressures on the aortic valve are calculated during valve opening and closing. The rate of work done by shear and the energy dissipation in the ventricle are small. The Bernoulli energy flux delivered to blood from ventricular contraction is practically balanced by energy flux at the aortic root and the rate change of kinetic energy in the ventricle.

An immersogeometric variational framework for fluid-structure interaction: application to bioprosthetic heart valves

In this paper, we develop a geometrically flexible technique for computational fluid-structure interaction (FSI). The motivating application is the simulation of tri-leaflet bioprosthetic heart valve function over the complete cardiac cycle. Due to the complex motion of the heart valve leaflets, the fluid domain undergoes large deformations, including changes of topology. The proposed method directly analyzes a spline-based surface representation of the structure by immersing it into a non-boundary-fitted discretization of the surrounding fluid domain. This places our method within an emerging class of computational techniques that aim to capture geometry on non-boundary-fitted analysis meshes. We introduce the term "immersogeometric analysis" to identify this paradigm. The framework starts with an augmented Lagrangian formulation for FSI that enforces kinematic constraints with a combination of Lagrange multipliers and penalty forces. For immersed volumetric objects, we formally eliminate the multiplier field by substituting a fluid-structure interface traction, arriving at Nitsche's method for enforcing Dirichlet boundary conditions on object surfaces. For immersed thin shell structures modeled geometrically as surfaces, the tractions from opposite sides cancel due to the continuity of the background fluid solution space, leaving a penalty method. Application to a bioprosthetic heart valve, where there is a large pressure jump across the leaflets, reveals shortcomings of the penalty approach. To counteract steep pressure gradients through the structure without the conditioning problems that accompany strong penalty forces, we resurrect the Lagrange multiplier field. Further, since the fluid discretization is not tailored to the structure geometry, there is a significant error in the approximation of pressure discontinuities across the shell. This error becomes especially troublesome in residual-based stabilized methods for incompressible flow, leading to problematic compressibility at practical levels of refinement. We modify existing stabilized methods to improve performance. To evaluate the accuracy of the proposed methods, we test them on benchmark problems and compare the results with those of established boundary-fitted techniques. Finally, we simulate the coupling of the bioprosthetic heart valve and the surrounding blood flow under physiological conditions, demonstrating the effectiveness of the proposed techniques in practical computations.

Bioengineering Strategies for Polymeric Scaffold for Tissue Engineering an Aortic Heart Valve: An Update

The occurrence of dysfunctional aortic valves is increasing every year, and current replacement heart valves, although having been shown to be clinically successful, are only short-term solutions and suffer from many agonizing long-term drawbacks. The tissue engineering of heart valves is recognized as one of the most promising answers for aortic valve disease therapy, but overcoming current shortcomings will require multidisciplinary efforts. The use of a polymeric scaffold to guide the growth of the tissue is the most common approach to generate a new tissue for an aortic heart valve. However, optimizing the design of the scaffold, in terms of biocompatibility, surface morphology for cell attachments and the correct rate of degradation is critical in creating a viable tissue-engineered aortic heart valve. This paper highlights the bioengineering strategies that need to be followed to construct a polymeric scaffold of sufficient mechanical integrity, with superior surface morphologies, that is capable of mimicking the valve dynamics in vivo. The current challenges and future directions of research for creating tissue-engineered aortic heart valves are also discussed.

A fluid-structure interaction model of the aortic valve with coaptation and compliant aortic root

While aortic valve root compliance and leaflet coaptation have significant influence on valve closure, their implications have not yet been fully evaluated. The present study developed a full fluid-structure interaction (FSI) model that is able to cope with arbitrary coaptation between the leaflets of the aortic valve during the closing phase. Two simplifications were also evaluated for the simulation of the closing phase only. One employs an FSI model with a rigid root and the other uses a "dry" (without flow) model. Numerical tests were performed to verify the model. New metrics were defined to process the results in terms of leaflet coaptation area and contact pressure. The axial displacement of the leaflets, closure time and coaptation parameters were similar in the two FSI models, whereas the dry model, with imposed uniform load on the leaflets, produced larger coaptation area and contact pressure, larger axial displacement and faster closure time compared with the FSI model. The differences were up to 30% in the coaptation area, 55% in the contact pressure and 170% in the closure time. Consequently, an FSI model should be used to accurately resolve the kinematics of the aortic valve and leaflet coaptation details during the end-closing stage.

FLOW DYNAMIC COMPARISON BETWEEN RECESSED HINGE AND OPEN PIVOT BI-LEAFLET HEART VALVE DESIGNS

The flow dynamics through the peripheral and hinge regions of a bi-leaflet mechanical heart valve are complex and result in abnormally high shear stresses particularly during the closing phase of the valve function. It has been observed that, the late stages of closure is more significant in the dynamics of platelet activation; therefore, the later stages of closure is simulated by solving the two-dimensional Navier-Stokes equations using an Eulerian Levelset based sharp interface Cartesian grid method. Using a fixed Cartesian mesh incorporating local mesh refinement for solution accuracy and efficiency, the flow through and within a recessed hinge design and an open pivot hinge design is compared. Platelets are modelled as point particles by Lagrangian particle tracking algorithm with one way coupling. A dilute particle flow is assumed and particle-particle interactions are neglected. It was observed that the hinge region of the open pivot valve indicated a lower potential for activation of platelets compared to that in valves with a recessed hinge design.

Prosthetic aortic heart valves: modeling and design

Although heart valve replacement is among the most common cardiovascular surgical procedures, their outcome is often difficult to predict. One of the reasons is the design and choice of the materials used for the fabrication of the prostheses. This review paper describes the use of modeling techniques in prosthetic heart valve (HV) design and aims at the justification and development of a polymer based trileaflet mechanical heart valve (MHV). The closing/opening phase behavior of the bileaflet MHV was investigated. The potential problem of valve failure due to crack propagation in the brittle pyrolytic carbon leaflet was also discussed. These studies suggest that although bileaflet MHV performs satisfactorily, there are justifications for improvement. Since the native aortic HV is trileaflet and made of anisotropic and hyperelastic tissue, one possible approach to a better MHV design is based on our ability to closely mimic the natural geometry and biomaterial properties.

Cavitation phenomena in mechanical heart valves: studied by using a physical impinging rod system

When studying mechanical heart valve cavitation, a physical model allows direct flow field and pressure measurements that are difficult to perform with actual valves, as well as separate testing of water hammer and squeeze flow effects. Movable rods of 5 and 10 mm diameter impinged same-sized stationary rods to simulate squeeze flow. A 24 mm piston within a tube simulated water hammer. Adding a 5 mm stationary rod within the tube generated both effects simultaneously. Charged-coupled device (CCD) laser displacement sensors, strobe lighting technique, laser Doppler velocimetry (LDV), particle image velocimetry (PIV) and high fidelity piezoelectric pressure transducers measured impact velocities, cavitation images, squeeze flow velocities, vortices, and pressure changes at impact, respectively. The movable rods created cavitation at critical impact velocities of 1.6 and 1.2 m/s; squeeze flow velocities were 2.8 and 4.64 m/s. The isolated water hammer created cavitation at 1.3 m/s piston speed. The combined piston and stationary rod created cavitation at an impact speed of 0.9 m/s and squeeze flow of 3.2 m/s. These results show squeeze flow alone caused cavitation, notably at lower impact velocity as contact area increased. Water hammer alone also caused cavitation with faster displacement. Both effects together were additive. The pressure change at the vortex center was only 150 mmHg, which cannot generate the magnitude of pressure drop required for cavitation bubble formation. Cavitation occurred at 3-5 m/s squeeze flow, significantly different from the 14 m/s derived by Bernoulli's equation; the temporal acceleration of unsteady flow requires further study.

On the biomechanics of heart valve function

Heart valves (HVs) are fluidic control components of the heart that ensure unidirectional blood flow during the cardiac cycle. However, this description does not adequately describe the biomechanical ramifications of their function in that their mechanics are multi-modal. Moreover, they must replicate their cyclic function over an entire lifetime, with an estimated total functional demand of least 3x10(9) cycles. The focus of the present review is on the functional biomechanics of heart valves. Thus, the focus of the present review is on functional biomechanics, referring primarily to biosolid as well as several key biofluid mechanical aspects underlying heart valve physiological function. Specifically, we refer to the mechanical behaviors of the extracellular matrix structural proteins, underlying cellular function, and their integrated relation to the major aspects of valvular hemodynamic function. While we focus on the work from the author's laboratories, relevant works of other investigators have been included whenever appropriate. We conclude with a summary of important future trends.

Two-dimensional simulation of flow and platelet dynamics in the hinge region of a mechanical heart valve

The hinge region of a mechanical bileaflet valve is implicated in blood damage and initiation of thrombus formation. Detailed fluid dynamic analysis in the complex geometry of the hinge region during the closing phase of the bileaflet valve is the focus of this study to understand the effect of fluid-induced stresses on the activation of platelets. A fixed-grid Cartesian mesh flow solver is used to simulate the blood flow through a two-dimensional geometry of the hinge region of a bileaflet mechanical valve. Use of local mesh refinement algorithm provides mesh adaptation based on the gradients of flow in the constricted geometry of the hinge. Leaflet motion is specified from the fluid-structure interaction analysis of the leaflet dynamics during the closing phase from a previous study, which focused on the fluid mechanics at the gap between the leaflet edges and the valve housing. A Lagrangian particle tracking method is used to model and track the platelets and to compute the magnitude of the shear stress on the platelets as they pass through the hinge region. Results show that there is a boundary layer separation in the gaps between the leaflet ear and the constricted hinge geometry. Separated shear layers roll up into vortical structures that lead to high residence times combined with exposure to high-shear stresses for particles in the hinge region. Particles are preferentially entrained into this recirculation zone, presenting the possibility of platelet activation, aggregation, and initiation of thrombi.

Curvilinear Immersed Boundary Method for Simulating Fluid Structure Interaction with Complex 3D Rigid Bodies

The sharp-interface CURVIB approach of Ge and Sotiropoulos [L. Ge, F. Sotiropoulos, A Numerical Method for Solving the 3D Unsteady Incompressible Navier-Stokes Equations in Curvilinear Domains with Complex Immersed Boundaries, Journal of Computational Physics 225 (2007) 1782-1809] is extended to simulate fluid structure interaction (FSI) problems involving complex 3D rigid bodies undergoing large structural displacements. The FSI solver adopts the partitioned FSI solution approach and both loose and strong coupling strategies are implemented. The interfaces between immersed bodies and the fluid are discretized with a Lagrangian grid and tracked with an explicit front-tracking approach. An efficient ray-tracing algorithm is developed to quickly identify the relationship between the background grid and the moving bodies. Numerical experiments are carried out for two FSI problems: vortex induced vibration of elastically mounted cylinders and flow through a bileaflet mechanical heart valve at physiologic conditions. For both cases the computed results are in excellent agreement with benchmark simulations and experimental measurements. The numerical experiments suggest that both the properties of the structure (mass, geometry) and the local flow conditions can play an important role in determining the stability of the FSI algorithm. Under certain conditions unconditionally unstable iteration schemes result even when strong coupling FSI is employed. For such cases, however, combining the strong-coupling iteration with under-relaxation in conjunction with the Aitken's acceleration technique is shown to effectively resolve the stability problems. A theoretical analysis is presented to explain the findings of the numerical experiments. It is shown that the ratio of the added mass to the mass of the structure as well as the sign of the local time rate of change of the force or moment imparted on the structure by the fluid determine the stability and convergence of the FSI algorithm. The stabilizing role of under-relaxation is also clarified and an upper bound of the required for stability under-relaxation coefficient is derived.

Two-Dimensional Dynamic Simulation of Platelet Activation During Mechanical Heart Valve Closure

A major drawback in the operation of mechanical heart valve prostheses is thrombus formation in the near valve region. Detailed flow analysis in this region during the valve closure phase is of interest in understanding the relationship between shear stress and platelet activation. A fixed-grid Cartesian mesh flow solver is used to simulate the blood flow through a bi-leaflet mechanical valve employing a two-dimensional geometry of the leaflet with a pivot point representing the hinge region. A local mesh refinement algorithm allows efficient and fast flow computations with mesh adaptation based on the gradients of the flow field in the leaflet-housing gap at the instant of valve closure. Leaflet motion is calculated dynamically based on the fluid forces acting on it employing a fluid-structure interaction algorithm. Platelets are modeled and tracked as point particles by a Lagrangian particle tracking method which incorporates the hemodynamic forces on the particles. A platelet activation model is included to predict regions which are prone to platelet activation. Closure time of the leaflet is validated against experimental studies. Results show that the orientation of the jet flow through the gap between the housing and the leaflet causes the boundary layer from the valve housing to be drawn in by the shear layer separating from the leaflet. The interaction between the separating shear layers is seen to cause a region of intensely rotating flow with high shear stress and high residence time of particles leading to high likelihood of platelet activation in that region.

Characteristics of pulsatile blood flow through the curved bileaflet mechanical heart valve installed in two different types of blood vessels: velocity and pressure of blood flow

The aim of this study was to investigate the flow fields of blood flowing through the curved bileaflet mechanical heart valve. A numerical analysis was carried out with the fluid-structure interaction between the blood flow and the motion of leaflets in two different types of blood vessels (type A, with sinus blood vessel, and type B, without sinus blood vessel). When the leaflet was fully opened, a fluttering phenomenon was detected in association with the blood flow, and recirculation flows were observed in the sinus region of the blood vessel for type A. During the closing phase, regurgitation was formed between the ring and the edge of the each leaflet for both types. When the leaflet came into contact with the valve ring at the end of the closing phase, rebound of the leaflet occurred. In consideration of the entire domain, the pressure drop occurs mainly in the valve region. The present results showed tendencies similar to those obtained by previous experiments for blood flow and contribute to the development of the curved bileaflet mechanical heart valve prostheses.

The hemodynamics of the Berlin pulsatile VAD and the role of its MHV configuration

The 3D flow in a model of the Berlin ventricular assist device (VAD) chamber with monoleaflet valves placed in S-shape conduits was simulated numerically. The blood flow dynamics were described in terms of flow patterns, velocity, pressure, and shear stress. The hemodynamic properties and the VAD's potential risk for thrombosis were evaluated in terms of mixing and washout properties, and global estimations of platelet level of activation (LOA). In order to evaluate the role of valves on the flow in the chamber, the flow in a model with bileaflet valves in straight conduits was simulated and compared with the original case. The results showed that in both models a large rotating flow was developed in the chamber during filling. This vortex filled the entire chamber and moved constantly up to the peak ejection phase, resulting in relatively low shear stress (up to 0.4 Pa) and no lasting stagnation regions. Significant shear stresses were found near the valves with higher values near the outlet valve in both models. The configuration of valves and conduits had a large effect on VAD washout and mixing properties, with advantage to the bileaflet model. However, since the bileaflet valves exhibited higher shear stresses, higher LOA were found for the bileaflet model.

Flow in prosthetic heart valves: state-of-the-art and future directions

Since the first successful implantation of a prosthetic heart valve four decades ago, over 50 different designs have been developed including both mechanical and bioprosthetic valves. Today, the most widely implanted design is the mechanical bileaflet, with over 170,000 implants worldwide each year. Several different mechanical valves are currently available and many of them have good bulk forward flow hemodynamics, with lower transvalvular pressure drops, larger effective orifice areas, and fewer regions of forward flow stasis than their earlier-generation counterparts such as the ball-and-cage and tilting-disc valves. However, mechanical valve implants suffer from complications resulting from thrombus deposition and patients implanted with these valves need to be under long-term anti-coagulant therapy. In general, blood thinners are not needed with bioprosthetic implants, but tissue valves suffer from structural failure with, an average life-time of 10-12 years, before replacement is needed. Flow-induced stresses on the formed elements in blood have been implicated in thrombus initiation within the mechanical valve prostheses. Regions of stress concentration on the leaflets during the complex motion of the leaflets have been implicated with structural failure of the leaflets with bioprosthetic valves. In vivo and in vitro experimental studies have yielded valuable information on the relationship between hemodynamic stresses and the problems associated with the implants. More recently, Computational Fluid Dynamics (CFD) has emerged as a promising tool, which, alongside experimentation, can yield insights of unprecedented detail into the hemodynamics of prosthetic heart valves. For CFD to realize its full potential, however, it must rely on numerical techniques that can handle the enormous geometrical complexities of prosthetic devices with spatial and temporal resolution sufficiently high to accurately capture all hemodynamically relevant scales of motion. Such algorithms do not exist today and their development should be a major research priority. For CFD to further gain the confidence of valve designers and medical practitioners it must also undergo comprehensive validation with experimental data. Such validation requires the use of high-resolution flow measuring tools and techniques and the integration of experimental studies with CFD modeling.

Time-dependent analysis of leaflets in mechanical aortic bileaflet heart valves in closing phase using the finite strip method

BACKGROUND AND AIMS OF THE STUDY

Mechanical heart valves (MHV) are widely used to replace dysfunctional and failed heart valves. The bileaflet MHV design is very popular due to its superior hemodynamics. Since their introduction in 1977, the hemodynamics of bileaflet prostheses has been extensively studied. In this study the dynamic behaviour during the closing phase of a bileaflet MHV under normal physiological conditions has been investigated.

METHODS

Fluid analysis is based on the control volume with moving boundaries in the vicinity of the occluder. Unsteady continuity equation, unsteady momentum equation on the control volume and unsteady Bernoulli's equation have been used to calculate velocity of blood flow and force on the occluder tip. To solve the governing equations for the calculation of pressure and the related force, the finite strips method has been implemented. Only 32 strips are sufficient to calculate the force due to pressure on the leaflets. The equations of motion have been solved using the Runge-Kutta method in the fourth order.

RESULTS

The maximum velocity of the leakage flow in the closing phase falls within the range of 3.5-4.4 m/s. The maximum velocity of the occluder tip is in the range of 2.4-3.2 m/s. The backflow also exhibits oscillation similar to that of the occluder with net backflow rate in the range of 9.7-12.3 ml/beat. The impact force between occluder and its housing is in the range of 80-140 N and impact occurs during 33.1-41.0 ms and the leaflets are completely settled at 108-115 ms in the closing phase.

CONCLUSION

The finite strip method was implemented to study the closing phase of a bileaflet MHV. Results are consistent with the previous experimental data. This method is of general applicability to study dynamic behaviour of MHVs.

A concentrated parameter model for the human cardiovascular system including heart valve dynamics and atrioventricular interaction

Numerical modeling of the human cardiovascular system has always been an active research direction since the 19th century. In the past, various simulation models of different complexities were proposed for different research purposes. In this paper, an improved numerical model to study the dynamic function of the human circulation system is proposed. In the development of the mathematical model, the heart chambers are described with a variable elastance model. The systemic and pulmonary loops are described based on the resistance-compliance-inertia concept by considering local effects of flow friction, elasticity of blood vessels and inertia of blood in different segments of the blood vessels. As an advancement from previous models, heart valve dynamics and atrioventricular interaction, including atrial contraction and motion of the annulus fibrosus, are specifically modeled. With these improvements the developed model can predict several important features that were missing in previous numerical models, including regurgitant flow on heart valve closure, the value of E/A velocity ratio in mitral flow, the motion of the annulus fibrosus (called the KG diaphragm pumping action), etc. These features have important clinical meaning and their changes are often related to cardiovascular diseases. Successful simulation of these features enhances the accuracy of simulations of cardiovascular dynamics, and helps in clinical studies of cardiac function.

Numerical simulation of cardiovascular dynamics with healthy and diseased heart valves

This paper presents a new concentrated parameter model for cardiovascular dynamics that includes an innovative model of heart valve dynamics, which is embedded in the overall model of the four chambers of the heart and the systemic and pulmonary circulation loops. The heart chambers are described with a variable elastance model, and the systemic and pulmonary loops are described with modified Windkessel models. In modelling the heart valve dynamics, the various factors that influence the valve motion are examined, and the governing differential equation for valve motion is derived. The heart valve model includes the influence of the blood pressure effect, the friction effect from the tissue, and from blood motion. As improvement from previous works, the contribution of the blood vortex effect in the vicinity of the valve leaflets to valve motion is specially considered. The proposed model is then used in simulation of healthy and certain pathological conditions such as mitral valve stenosis and aortic regurgitation. The predicted results agree well with results illustrated in cardiology textbooks.

Mechanism for Cavitation in the Mechanical Heart Valve With an Artificial Heart: Nuclei and Viscosity Dependence

Until now, we have estimated cavitation for mechanical heart valves (MHV) mounted in an electrohydraulic total artificial heart (EHTAH) with tap water. However, tap water at room temperature is not a proper substitute for blood at 37 degrees C. We therefore investigated fluid characterization in studies of MHV cavitation associated with the viscosity and nuclei content of a testing fluid. We used the Medtronic Hall valve mounted in the mitral position of the EHTAH. As testing fluids, tap water, distilled water, and glycerin solution were used. The valve-closing velocity, pressure-drop measurements, and a high-speed video camera were employed to determine the cavitation intensity in MHV. Most of the cavitation bubbles were observed at the edge of the valve stop. Our analysis of the results indicates that squeeze flow is the major cause of cavitation in the Medtronic Hall valve. The cavitation intensity increased with increases in the fluid viscosity and the valve-closing velocity. Even if cavitation intensity in glycerin solution was greater, the cavitation occurrence probability was less in glycerin solution than in tap water. Our results suggest that tap water contains particles that cause an increase in the cavitation occurrence probability. We conclude that cavitation intensity is greatly affected by the nuclei concentration in the fluid and the fluid viscosity.

Fluid Mechanics of Heart Valves

Valvular heart disease is a life-threatening disease that afflicts millions of people worldwide and leads to approximately 250,000 valve repairs and/or replacements each year. Malfunction of a native valve impairs its efficient fluid mechanic/hemodynamic performance. Artificial heart valves have been used since 1960 to replace diseased native valves and have saved millions of lives. Unfortunately, despite four decades of use, these devices are less than ideal and lead to many complications. Many of these complications/problems are directly related to the fluid mechanics associated with the various mechanical and bioprosthetic valve designs. This review focuses on the state-of-the-art experimental and computational fluid mechanics of native and prosthetic heart valves in current clinical use. The fluid dynamic performance characteristics of caged-ball, tilting-disc, bileaflet mechanical valves and porcine and pericardial stented and nonstented bioprostheic valves are reviewed. Other issues related to heart valve performance, such as biomaterials, solid mechanics, tissue mechanics, and durability, are not addressed in this review.

Thrombogenicity of Mechanical Aortic Valves in an Animal Model: Site Specific Testing Is Crucial

We evaluated a new trileaflet prosthesis and a control bileaflet prosthesis in the mitral and aortic positions in 27 calves. The prototype trileaflet valve (TV1) functioned satisfactorily in the mitral position (TV1m, n = 7) but later yielded thrombogenic complications in the aortic position (TV1a, n = 4). The valve was redesigned (TV2) and retested in the mitral (TV2m n = 4) and aortic (TV2a, n = 5) positions, along with control valves (Cm, n = 4; Ca, n = 3). At necropsy, the valves were graded on a scale of 0 (no visible thrombi) to 4 (thrombi greater than 5 mm and/or obstructed leaflets). The TV1m, TV2m, and Cm animals, respectively, had implant durations of 215+/-112, 140+/-63, and 159+/-89 days and thrombus grades of 0.71+/-0.76, 0.33+/-0.58, and 1.50+/-0.58. The TV1a, TV2a, and Ca animals had implant durations of 18+/-12, 159+/-61, and 108+/-62 days and thrombus grades of 2.75+/-1.00, 0.50+/-0.58, and 0.67+/-0.58 (p < .005; TV2a vs. TV1a). Three TV1a calves died early of valve related complications. A design irregularity, undetected in the mitral position but revealed in the aortic position, caused a high early mortality in the TV1a animals. Redesigning the prosthesis eliminated valve related mortality and significantly reduced the thrombus grade. Because satisfactory performance in the mitral position does not guarantee safety or efficacy in the aortic position, site specific preclinical testing is crucial for mechanical heart valves.

Numerical simulation of flow in mechanical heart valves: grid resolution and the assumption of flow symmetry

A numerical method is developed for simulating unsteady, 3-D, laminar flow through a bileaflet mechanical heart valve with the leaflets fixed. The method employs a dual-time-stepping artificial-compressibility approach together with overset (Chimera) grids and is second-order accurate in space and time. Calculations are carried out for the full 3-D valve geometry under steady inflow conditions on meshes with a total number of nodes ranging from 4 x 10(5) to 1.6 x 10(6). The computed results show that downstream of the leaflets the flow is dominated by two pairs of counter-rotating vortices, which originate on either side of the central orifice in the aortic sinus and rotate such that the common flow of each pair is directed away from the aortic wall. These vortices intensify with Reynolds number, and at a Reynolds number of approximately 1200 their complex interaction leads to the onset of unsteady flow and the break of symmetry with respect to both geometric planes of symmetry. Our results show the highly 3-D structure of the flow; question the validity of computationally expedient assumptions of flow symmetry; and demonstrate the need for highly resolved, fully 3-D simulations if computational fluid dynamics is to accurately predict the flow in prosthetic mechanical heart valves.

Cavitation phenomena in mechanical heart valves: the role of squeeze flow velocity and contact area on cavitation initiation between two impinging rods

In this study, the closing dynamics of two impinging rods were experimentally analyzed to simulate the cavitation phenomena associated with mechanical heart valve closure. The purpose of this study was to investigate the cavitation phenomena with respect to squeeze flow between two impinging surfaces and the parameter that influences cavitation inception. High-speed flow imaging was employed to visualize and identify regions of cavitation. The images obtained favored squeeze flow as an important mechanism in cavitation inception. A correlation study of the effects of impact velocities, contact areas and squeeze flow velocity on cavitation inception showed that increasing impact velocities results in an increase in the risk of cavitation. It was also shown that for similar impact velocities, regions near the point of impact were found to cavitate later for those with smaller contact areas. It was found that the decrease in contact areas and squeeze flow velocities would delay the onset and reduce the intensity of cavitation. It is also interesting to note that the squeeze flow velocity alone does not provide an indication if cavitation inception will occur. This is corroborated by the wide range of published critical squeeze flow velocity required for cavitation inception. It should be noted that the temporal acceleration of fluid, often neglected in the literature, can also play an important role on cavitation inception for unsteady flow phenomenon. This is especially true in mechanical heart valves, where for the same leaflet closing velocity, valves with a seat stop were observed to cavitate earlier. Based on these results, important inferences may be made to the design of mechanical heart valves with regards to cavitation inception.