Three-dimensional diastolic blood flow in the left ventricle

Three-dimensional blood flow in a human left ventricle is studied via a computational analysis with magnetic resonance imaging of the cardiac motion. Formation, growth and decay of vortices during the myocardial dilation are analyzed with flow patterns on various diametric planes. They are dominated by momentum transfer during flow acceleration and deceleration through the mitral orifice. The posterior and anterior vortices form an asymmetric annular vortex at the mitral orifice, providing a smooth transition for the rapid inflow to the ventricle. The development of core vortex accommodates momentum for deceleration and for acceleration at end diastolic atrial contraction. The rate of energy dissipation and that of work done by viscous stresses are small; they are approximately balanced with each other. The kinetic energy flux and the rate of work done by pressure delivered to blood from ventricular dilation is well balanced by the total energy influx at the mitral orifice and the rate change of kinetic energy in the ventricle.     

Intraventricular flow patterns and stasis in the LVAD-assisted heart

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

Numerical Modeling of Intraventricular Flow during Diastole after Implantation of BMHV

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

Determinants of the occurrence of vortex rings in the left ventricle during diastole

This study employs classical inviscid fluid dynamics theory to investigate whether LV diastolic inflow volume and the size of the LV play a role in vortex ring formation. Fluid injection across an orifice into a large container results in the generation of a vortex ring having a constant size and speed. Relations between the vortex size and speed and the injection were obtained by applying conservation laws regarding kinetic energy, impulse and vorticity; the initial state was computed using a bolus injection model, and the final state by using the Kelvin vortex model. An important parameter in the equations is the relative injection length, i.e., the ratio of the length of the injected bolus and the radius of the orifice (L/R). Its estimated highest value in man, L/R = 15, produces a rather thick vortex ring (relative thickness 0.77). Comparable results following from the Hill vortex model convinced us that the Kelvin vortex model can be applied in the whole range of injection lengths in the human left ventricle. In an in vitro model it is shown experimentally that vortex rings can be generated for L/R in the range from 2 to 16. The measured traveling speed of the vortex ring is in fair agreement with the theory, as well as the ring radius for large injections. A vortex ring located in a narrow channel cannot reach its proper traveling speed. The method of images is used to estimate the speed reduction of vortex rings within a cylinder. It turns out that propagation of vortex rings is possible when the ratio of orifice to cylinder radius is less than about 0.5.(ABSTRACT TRUNCATED AT 250 WORDS)     

Computational analysis of blood flow in an integrated model of the left ventricle and the aorta

To study the effects of intraventricular flow dynamics on the aortic flow, we created an integrated model of the left ventricle and aorta and conducted a computer simulation of diastolic and systolic blood flow within this model. The results demonstrated that the velocity profile at the aortic annulus changed dynamically, and was influenced by the intraventricular flow dynamics. The profile was almost flat in early systole but became nonuniform as systole progressed, and was skewed toward the posterior side in midsystole and toward the anterior side in later systole. At a distance from the aortic annulus, a different velocity profile was induced by the twisting and torsion of the aorta. In the ascending aorta, the fastest flow was initially located in the posteromedial sector, and it moved to the posterior section along the circumference as systole progressed. The nonuniformity of the aortic inflow gave rise to a complex wall shear stress (WSS) distribution in the aorta. A comparison of the WSS distribution obtained in this integrated analysis with that obtained in flow calculations using an isolated aorta model with Poiseuille and flat inlet conditions showed that intraventricular flow affected the WSS distribution in the ascending aorta. These results address the importance of an integrated analysis of flow in the left ventricle and aorta.     

Particle flow behavior in models of branching vessels. II. Effects of branching angle and diameter ratio on flow patterns

To further elucidate the role of fluid mechanical factors in the localization of atherogenesis and thrombogenesis, we have studied the 3-dimensional flow patterns in square T-junctions with branching angles theta from 30 degrees to 150 degrees and diameter ratios d/D (side: main tube) from 1.05/3.0 to 1.0. Cine films of the motions of tracer microspheres in dilute suspensions were taken at inflow Reynolds numbers from 15 to 400 and flow ratios (main: side tube) from 0.1 to 4.0. Flow patterns with suspension entering through the main tube were similar to those previously described in uniform 3 mm diameter T-junctions: paired vortices (spiral secondary flows) symmetrical about the common median plane formed at the entrances of the main and side daughter tubes. Particles circulated through the main vortex, some crossing above and below the mainstream into and through the side vortex. At the geometrical flow ratio, the main vortex became smaller and smaller as the branching angle (theta less than 90 degrees) and diameter ratio decreased, and was confined to a thin side tube was a minimum. In obtuse angle T-junctions the stagnation point shifted from the flow divider into the side tube, enhancing the flow disturbance there. The velocity distributions in main and side tubes were skewed towards the inner walls close to the flow divider. When flow entered through the side tube, a pair of recirculation zones formed in the main tube at the inner wall of the bend with a sharper angle.     

A new approach for the evaluation of the severity of coarctation of the aorta using Doppler velocity index and effective orifice area: in vitro validation and clinical implications

Early detection and accurate estimation of COA severity are the most important predictors of successful long-term outcome. However, current clinical parameters used for the evaluation of the severity of COA have several limitations and are flow dependent. The objectives of this study are to evaluate the limitations of current existing parameters for the evaluation of the severity of coarctation of the aorta (COA) and suggest two new parameters: COA Doppler velocity index and COA effective orifice area. Three different severities of COAs were tested in a mock flow circulation model under various flow conditions and in the presence of normal and stenotic aortic valves. Catheter trans-COA pressure gradients and Doppler echocardiographic trans-COA pressure gradients were evaluated. COA Doppler velocity index was defined as the ratio of pre-COA to post-COA peak velocities measured by Doppler echocardiography. COA Doppler effective orifice area was determined using continuity equation. The results show that peak-to-peak trans-COA pressure gradient significantly increased with flow rate (from 83% to 85%). Peak Doppler pressure gradient also significantly increased with flow rate (80-85%). A stenotic or bicuspid aortic valve increased peak Doppler pressure gradient by 20-50% for a COA severity of 75%. Both COA Doppler velocity index and COA effective orifice area did not demonstrate significant flow dependence or dependence upon aortic valve condition. As a conclusion, COA Doppler velocity index and COA effective orifice area are flow independent and do not depend on aortic valve conditions. They can, then, more accurately predict the severity of COA.     

Mechanisms of blood flow during pneumatic vest cardiopulmonary resuscitation

Mechanisms of blood flow during cardiopulmonary resuscitation (CPR) were studied in a canine model with implanted mitral and aortic flow probes and by use of cineangiography. Intrathoracic pressure (ITP) fluctuations were induced by a circumferential pneumatic vest, with and without simultaneous ventilation, and by use of positive-pressure ventilation alone. Vascular volume and compression rate were altered with each CPR mode. Antegrade mitral flow was interpreted as left ventricular (LV) inflow, and antegrade aortic flow was interpreted as LV outflow. The pneumatic vest was expected to elevate ITP uniformly and thus produce simultaneous LV inflow and LV outflow throughout compression. This pattern, the passive conduit of "thoracic pump" physiology, was unequivocally demonstrated only during ITP elevation with positive-pressure ventilation alone at slow rates. During vest CPR, LV outflow started promptly with the onset of compression, whereas LV inflow was delayed. At compression rates of 50 times/min and normal vascular filling pressures, the delay was sufficiently long that all LV filling occurred with release of compression. This is the pattern that would be expected with direct LV compression or "cardiac pump" physiology. During the early part of the compression phase, catheter tip transducer LV and left atrial pressure measurements demonstrated gradients necessitating mitral valve closure, while cineangiography showed dye droplets moving from the large pulmonary veins retrograde to the small pulmonary veins. When the compression rate was reduced and/or when intravascular pressures were raised with volume infusion, LV inflow was observed at some point during the compressive phase. Thus, under these conditions, features of both thoracic pump and cardiac pump physiology occurred within the same compression. Our findings are not explained by the conventional conceptions of either thoracic pump or cardiac compression CPR mechanisms alone.     

Numerical and experimental investigations of pulsatile blood flow pattern through a dysfunctional mechanical heart valve

Around 250,000 heart valve replacements are performed every year around the world. Due their higher durability, approximately 2/3 of these replacements use mechanical prosthetic heart valves (mainly bileaflet valves). Although very efficient, these valves can be subject to valve leaflet malfunctions. These malfunctions are usually the consequence of pannus ingrowth and/or thrombus formation and represent serious and potentially fatal complications. Hence, it is important to investigate the flow field downstream of a dysfunctional mechanical heart valve to better understand its impact on blood components (red blood cells, platelets and coagulation factors) and to improve the current diagnosis techniques. Therefore, the objective of this study will be to numerically and experimentally investigate the pulsatile turbulent flow downstream of a dysfunctional bileaflet mechanical heart valve in terms of velocity field, vortex formation and potential negative effect on blood components. The results show that the flow downstream of a dysfunctional valve was characterized by abnormally elevated velocities and shear stresses as well as large scale vortices. These characteristics can predispose to blood components damage. Furthermore, valve malfunction led to an underestimation of maximal transvalvular pressure gradient, using Doppler echocardiography, when compared to numerical results. This could be explained by the shifting of the maximal velocity towards the normally functioning leaflet. As a consequence, clinicians should try, when possible, to check the maximal velocity position not only at the central orifice but also through the lateral orifices. Finding the maximal velocity in the lateral orifice could be an indication of valve dysfunction.     

Direct numerical simulation of a 2D-stented aortic heart valve at physiological flow rates

We study the nonlinear interaction of an aortic heart valve, composed of hyperelastic corrugated leaflets of finite density attached to a stented vessel under physiological flow conditions. In our numerical simulations, we use a 2D idealised representation of this arrangement. Blood flow is caused by a time-varying pressure gradient that mimics that of the aortic valve and corresponds to a peak Reynolds number equal to 4050. Here, we fully account for the shear-thinning behaviour of the blood and large deformations and contact between the leaflets by solving the momentum and mass balances for blood and leaflets. The mixed finite element/Galerkin method along with linear discontinuous Lagrange multipliers for coupling the fluid and elastic domains is adopted. Moreover, a series of challenging numerical issues such as the finite length of the computational domain and the conditions that should be imposed on its inflow/outflow boundaries, the accurate time integration of the parabolic and hyperbolic momentum equations, the contact between the leaflets and the non-conforming mesh refinement in part of the domain are successfully resolved. Calculations for the velocity and the shear stress fields of the blood reveal that boundary layers appear on both sides of a leaflet. The one along the ventricular side transfers blood with high momentum from the core region of the vessel to the annulus or the sinusoidal expansion, causing the continuous development of flow instabilities. At peak systole, vortices are convected in the flow direction along the annulus of the vessel, whereas during the closure stage of the valve, an extremely large vortex develops in each half of the flow domain.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

Role of blood cell-wall interactions in thrombogenesis and atherogenesis: a microrheological study

The relationship between blood flow and the localization of thrombosis and atherosclerosis in vivo was investigated using the approach and techniques of microrheology. The flow patterns and wall-adhesion of platelets were studied in the captive annular vortex formed at a sudden tubular expansion at various hematocrits in steady and pulsatile flow. The adhesion density exhibited a peak within the vortex and just downstream of the reattachment point, which is also a stagnation point. The peaks flattened out with increasing Reynolds number in steady flow and also in pulsatile flow. Platelet adhesion increased markedly with increasing hematocrit. The localization of adhesion peaks was explained by curvature of the streamlines carrying platelets to the wall on either side of the reattachment point. The relevance of these results to the circulation is that stagnation points are found in regions of disturbed flow at various sites in the arterial and venous circulations. This was shown in experiments using a technique whereby flow was visualized in isolated transparent natural blood vessels prepared from dogs and humans postmortem. In dog saphenous vein bileaflet valves, there was a large primary spiral vortex as well as a smaller secondary vortex, the latter acting as a trap and generator of thrombi. Recirculation zones also existed in the dog aorta at T-junctions of the celiac, cranial mesenteric and renal arteries. Finally, in the human carotid bifurcation, a large standing recirculation zone consisting of spiral secondary flows formed in the carotid sinus at physiological flow conditions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effect of systolic flow rate on the prediction of effective prosthetic valve orifice area

The Gorlin equation for the hemodynamic assessment of valve area is commonly used in cardiac catheterization laboratories. A study was performed to test the prediction capabilities of the Gorlin formula as well as the Aaslid and Gabbay formula for the effective orifice area of prosthetic heart valves. Pressure gradient, flow, and valve opening area measurements were performed on four 27 mm valve prostheses (two mechanical bileaflet designs, St. Jude and Edwards-Duromedics, an Edwards pericardial tissue valve, and a trileaflet polyurethane valve) each mounted in the aortic position of an in vitro pulse duplicator. With the known valve orifice area, a different discharge coefficient was computed for each of the four valves and three orifice area formulas. After some theoretical considerations, it was proposed that the discharge coefficient would be a function of the flow rate through the valve. All discharge coefficients were observed to increase with increasing systolic flow rate. An empirical relationship of discharge coefficient as a linear function of systolic flow rate was determined through a regression analysis, with a different relationship for each valve and each orifice area formula. Using this relationship in the orifice area formulas improved the accuracy of the prediction of the effective orifice area with all three formulas performing equally well.     

Interactions between the right ventricle and pulmonary vasculature in the fetus.

A midsystolic plateau differentiates the pattern of fetal pulmonary trunk blood flow from aortic flow. To determine whether this plateau arises from interactions between the left (LV) and right ventricle (RV) via the ductus arteriosus or from interactions between the RV and the lung vasculature, we measured blood flows and pressures in the pulmonary trunk and aorta of eight anesthetized (ketamine and alpha-chloralose) fetal lambs. Wave-intensity analysis revealed waves of energy traveling forward, away from the LV and the RV early in systole. During midsystole, a wave of energy traveling back toward the RV decreased blood flow velocity from the RV and produced the plateau in blood flow. Calculations revealed that this backward-traveling wave originated as a forward-traveling wave generated by the RV that was reflected from the lung vasculature back toward the heart and not as a forward-traveling wave generated by the LV that crossed the ductus arteriosus. Elimination of this backward-traveling wave and its associated effect on RV flow may be an important component of the increase in RV output that accompanies birth.     

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.     

Abnormal cardiac inflow patterns during postnatal development in a mouse model of Holt-Oram syndrome

Tbx5(del/+) mice provide a model of human Holt-Oram syndrome. In this study, the cardiac functional phenotypes of this mouse model were investigated with 30-MHz ultrasound by comparing 12 Tbx5(del/+) mice with 12 wild-type littermates at 1, 2, 4, and 8 wk of age. Cardiac dimensions were measured with two-dimensional and M-mode imaging. The flow patterns in the left and right ventricular inflow channels were evaluated with Doppler flow sampling. Compared with wild-type littermates, Tbx5(del/+) mice showed significant changes in the mitral flow pattern, including decreased peak velocity of the left ventricular (LV) early filling wave (E wave), increased peak velocity of the late filling wave (A wave), and decreased or even reversed peak E-to-A ratio. The prolongation of LV isovolumic relaxation time was detected in Tbx5(del/+) neonates as early as 1 wk of age. In Tbx5(del/+) mice, LV wall thickness appeared normal but LV chamber dimension was significantly reduced. LV systolic function did not differ from that in wild-type littermates. In contrast, the Doppler flow spectrum in the enlarged tricuspid orifice of Tbx5(del/+) mice demonstrated increased peak velocities of both E and A waves and increased total time-velocity integral but unchanged peak E/A. In another 13 mice (7 Tbx5(del/+), 6 wild-type) at 2 wk of age, significant correlation was found between Tbx5 gene expression level in ventricular myocardium and LV filling parameters. In conclusion, the LV diastolic function of Tbx5(del/+) mice is significantly deteriorated, whereas the systolic function remains normal.     

Diastolic right ventricular filling vortex in normal and volume overload states

Functional imaging computational fluid dynamics simulations of right ventricular (RV) inflow fields were obtained by comprehensive software using individual animal-specific dynamic imaging data input from three-dimensional (3-D) real-time echocardiography (RT3D) on a CRAY T-90 supercomputer. Chronically instrumented, lightly sedated awake dogs (n = 7) with normal wall motion (NWM) at control and normal or diastolic paradoxical septal motion (PSM) during RV volume overload were investigated. Up to the E-wave peak, instantaneous inflow streamlines extended from the tricuspid orifice to the RV endocardial surface in an expanding fanlike pattern. During the descending limb of the E-wave, large-scale (macroscopic or global) vortical motions ensued within the filling RV chamber. Both at control and during RV volume overload (with or without PSM), blood streams rolled up from regions near the walls toward the base. The extent and strength of the ring vortex surrounding the main stream were reduced with chamber dilatation. A hypothesis is proposed for a facilitatory role of the diastolic vortex for ventricular filling. The filling vortex supports filling by shunting inflow kinetic energy, which would otherwise contribute to an inflow-impeding convective pressure rise between inflow orifice and the large endocardial surface of the expanding chamber, into the rotational kinetic energy of the vortical motion that is destined to be dissipated as heat. The basic information presented should improve application and interpretation of noninvasive (Doppler color flow mapping, velocity-encoded cine magnetic resonance imaging, etc.) diastolic diagnostic studies and lead to improved understanding and recognition of subtle, flow-associated abnormalities in ventricular dilatation and remodeling.     

Quantification of presystolic blood flow organization and energetics in the human left ventricle

Intracardiac blood flow patterns are potentially important to cardiac pumping efficiency. However, these complex flow patterns remain incompletely characterized both in health and disease. We hypothesized that normal left ventricular (LV) blood flow patterns would preferentially optimize a portion of the end-diastolic volume (LVEDV) for effective and rapid systolic ejection by virtue of location near and motion towards the LV outflow tract (LVOT). Three-dimensional cine velocity and morphological data were acquired in 12 healthy persons and 1 patient with dilated cardiomyopathy using MRI. A previously validated method was used for analysis in which the LVEDV was separated into four functional flow components based on the blood's locations at the beginning and end of the cardiac cycle. Each component's volume, kinetic energy (KE), site, direction, and linear momentum relative to the LVOT were calculated. Of the four components, the LV inflow that passes directly to outflow in a single cardiac cycle (Direct Flow) had the largest volume. At the time of isovolumic contraction, Direct Flow had the greatest amount of KE and the most favorable combination of distance, angle, and linear momentum relative to the LVOT. Atrial contraction boosted the late diastolic KE of the ejected components. We conclude that normal diastolic LV flow creates favorable conditions for ensuing ejection, defined by proximity and energetics, for the Direct Flow, and that atrial contraction augments the end-diastolic KE of the ejection volume. The correlation of Direct Flow characteristics with ejection efficiency might be a relevant investigative target in cardiac failure.