It is a widely accepted axiom that localized concentration of mechanical stress and large flexural deformation is closely related to the calcification and tissue degeneration in bioprosthetic heart valves (BHV). In order to investigate the complex BHV deformations and stress distributions throughout the cardiac cycle, it is necessary to perform an accurate dynamic analysis with a morphologically and physiologically realistic material specification for the leaflets. We have developed a stress resultant shell model for BHV leaflets incorporating a Fung-elastic constitutive model for in-plane and bending responses separately. Validation studies were performed by comparing the finite element predicted displacement and strain measures with the experimentally measured data under physiological pressure loads. Computed regions of stress concentration and large flexural deformation during the opening and closing phases of the cardiac cycle correlated with previously reported regions of calcification and/or mechanical damage on BHV leaflets. It is expected that the developed experimental and computational methodology will aid in the understanding of the complex dynamic behavior of native and bioprosthetic valves and in the development of tissue engineered valve substitutes.

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.

The wall shear stress induced by the leaflet motion during the valve-closing phase has been implicated with thrombus initiation with prosthetic valves. Detailed flow dynamic analysis in the vicinity of the leaflets and the housing during the valve-closure phase is of interest in understanding this relationship. A three-dimensional unsteady flow analysis past bileaflet valve prosthesis in the mitral position is presented incorporating a fluid-structure interaction algorithm for leaflet motion during the valve-closing phase. Arbitrary Lagrangian-Eulerian method is employed for incorporating the leaflet motion. The forces exerted by the fluid on the leaflets are computed and applied to the leaflet equation of motion to predict the leaflet position. Relatively large velocities are computed in the valve clearance region between the valve housing and the leaflet edge with the resulting relatively large wall shear stresses at the leaflet edge during the impact-rebound duration. Negative pressure transients are computed on the surface of the leaflets on the atrial side of the valve, with larger magnitudes at the leaflet edge during the closing and rebound as well. Vortical flow development is observed on the inflow (atrial) side during the valve impact-rebound phase in a location central to the leaflet and away from the clearance region where cavitation bubbles have been visualized in previously reported experimental studies.

No abstract

The relationships among vascular geometry, hemodynamics, and plaque development in the coronary arteries are complex and not yet well understood. This paper reports a methodology for the quantitative analysis of in vivo coronary morphology and hemodynamics, with particular emphasis placed on the critical issues of image segmentation and the automated classification of disease severity. We were motivated by the observation that plaque more often developed at the inner curvature of a vessel, presumably due to the relatively lower wall shear stress at these locations. The presented studies are based on our validated methodology for the three-dimensional fusion of intravascular ultrasound (IVUS) and X-ray angiography, introducing a novel approach for IVUS segmentation that incorporates a robust, knowledge-based cost function and a fully optimal, threedimensional segmentation algorithm. Our first study shows that circumferential plaque distribution depends on local vessel curvature in the majority of vessels. The second study analyzes the correlation between plaque distribution and wall shear stress in a set of 48 in vivo vessel segments. The results were conclusive for both studies, with a stronger correlation of circumferential plaque thickness with local curvature than with wall shear stress. The inverse relationship between local wall shear stress and plaque thickness was significantly more pronounced (p < 0.025) in vessel cross sections exhibiting compensatory enlargement (positive remodeling) without luminal narrowing than when the full spectrum of disease severity was considered. The inverse relationship was no longer observed in vessels where less than 35% of vessel cross sections remained without luminal narrowing. The findings of this study confirm, in vivo, the hypothesis that relatively lower wall shear stress is associated with early plaque development.

Tissue engineered heart valves (TEHV) have been observed to respond to mechanical conditioning in vitro by expression of activated myofibroblast phenotypes followed by improvements in tissue maturation. In separate studies, cyclic flexure, stretch, and flow (FSF) have been demonstrated to exhibit both independent and coupled stimulatory effects. Synthesis of these observations into a rational framework for TEHV mechanical conditioning has been limited, however, due to the functional complexity of trileaflet valves and the inherent differences of separate bioreactor systems. Toward quantifying the effects of individual mechanical stimuli similar to those that occur during normal valve function, a novel bioreactor was developed in which FSF mechanical stimuli can be applied to engineered heart valve tissues independently or in combination. The FSF bioreactor consists of two identically equipped chambers, each having the capacity to hold up to 12 rectangular tissue specimens (25 × 7.5 × 1 mm) via a novel "spiralbound" technique. Specimens can be subjected to changes-in-curvature up to 50 mm −1 and uniaxial tensile strains up to 75%. Steady laminar flow can be applied by a magnetically coupled paddlewheel system. Computational fluid dynamic (CFD) simulations were conducted and experimentally validated by particle image velocimetry (PIV). Tissue specimen wall shear stress profiles were predicted as a function of paddlewheel speed, culture medium viscosity, and the quasi-static state of specimen deformation (i.e., either undeformed or completely flexed). Velocity profiles predicted by 2D CFD simulations of the paddlewheel mechanism compared well with PIV measurements, and were used to determine boundary conditions in localized 3D simulations. For undeformed specimens, predicted inter-specimen variations in wall shear stress were on average ±7%, with an average wall shear stress of 1.145 dyne/cm 2 predicted at a paddlewheel speed of 2000 rpm and standard culture conditions. In contrast, while the average wall shear stress predicted for specimens in the quasi-static flexed state was ~59% higher (1.821 dyne/cm 2 ), flexed specimens exhibited a broad intra-specimen wall shear stress distribution between the convex and concave sides that correlated with specimen curvature, with peak wall shear stresses of ~10 dyne/ cm 2 . This result suggests that by utilizing simple flexed geometric configurations, the present system can also be used to study the effects of spatially varying shear stresses. We conclude that the present design provides a robust tool for the study of mechanical stimuli on in vitro engineered heart valve tissue formation.

Platelet activation, adhesion, and aggregation on the blood vessel and implants result in the formation of mural thrombi. Platelet dynamics in blood flow is influenced by the far more numerous erythrocytes (RBCs). This is particularly the case in the smaller blood vessels (arterioles) and in constricted regions of blood flow (such as in valve leakage and hinge regions) where the dimensions of formed elements of blood become comparable with that of the flow geometry. In such regions, models to predict platelet motion, activation, aggregation and adhesion must account for platelet-RBC interactions. This paper studies platelet-RBC interactions in shear flows by performing simulations of micro-scale dynamics using a computational fluid dynamics (CFD) model. A level-set sharp-interface immersed boundary method is employed in the computations in which RBC and platelet boundaries are tracked on a two-dimensional Cartesian grid. The RBCs are assumed to have an elliptical shape and to deform elastically under fluid forces while the platelets are assumed to behave as rigid particles of circular shape. Forces and torques between colliding blood cells are modeled using an extension of the soft-sphere model for elliptical particles. RBCs and platelets are transported under the forces and torques induced by fluid flow and cell-cell and cell-platelet collisions. The simulations show that platelet migration toward the wall is enhanced with increasing hematocrit, in agreement with past experimental observations. This margination is seen to occur due to hydrodynamic forces rather than collisional forces or volumetric exclusion effects. The effect of fluid shear forces on the platelets increases exponentially as a function of hematocrit for the range of parameters covered in this study. The micro-scale analysis can be potentially employed to obtain a deterministic relationship between fluid forces and platelet activation and aggregation in blood flow past cardiovascular implants.

Most surgical procedures for patients with mitral regurgitation (MR) focus on optimization of annular dimension and shape utilizing ring annuloplasty to restore normal annular geometry, increase leaflet coaptation, and reduce regurgitation. Computational studies may provide insight on the effect of annular motion on mitral valve (MV) function through the incorporation of patient-specific MV apparatus geometry from clinical imaging modalities such as echocardiography. In the present study, we have developed a novel algorithm for modeling patient-specific annular motion across the cardiac cycle to further improve our virtual MV modeling and simulation strategy. The MV apparatus including the leaflets, annulus, and location of papillary muscle tips was identified using patient 3D echocardiography data at end diastole and peak systole and converted to virtual MV model. Dynamic annular motion was modeled by incorporating the ECG-gated time-varying scaled annular displacement across the cardiac cycle. We performed dynamic finite element (FE) simulation of two sets of patient data with respect to the presence of MR. Annular morphology, stress distribution across the leaflets and annulus, and contact stress distribution were determined to assess the effect of annular motion on MV function and leaflet coaptation. The effect of dynamic annular motion clearly demonstrated reduced regions with large stress values and provided an improved accuracy in determining the location of improper leaflet coaptation. This strategy has the potential to better quantitate the extent of pathologic MV and better evaluate functional restoration following MV repair.