Echocardiographic Particle Image Velocimetry: A Novel Technique for Quantification of Left Ventricular Blood Vorticity Pattern

A novel left heart simulator for the multi-modality characterization of native mitral valve geometry and fluid mechanics

Numerical models of the mitral valve have been used to elucidate mitral valve function and mechanics. These models have evolved from simple two-dimensional approximations to complex three-dimensional fully coupled fluid structure interaction models. However, to date these models lack direct one-to-one experimental validation. As computational solvers vary considerably, experimental benchmark data are critically important to ensure model accuracy. In this study, a novel left heart simulator was designed specifically for the validation of numerical mitral valve models. Several distinct experimental techniques were collectively performed to resolve mitral valve geometry and hemodynamics. In particular, micro-computed tomography was used to obtain accurate and high-resolution (39 μm voxel) native valvular anatomy, which included the mitral leaflets, chordae tendinae, and papillary muscles. Three-dimensional echocardiography was used to obtain systolic leaflet geometry. Stereoscopic digital particle image velocimetry provided all three components of fluid velocity through the mitral valve, resolved every 25 ms in the cardiac cycle. A strong central filling jet (V ~ 0.6 m/s) was observed during peak systole with minimal out-of-plane velocities. In addition, physiologic hemodynamic boundary conditions were defined and all data were synchronously acquired through a central trigger. Finally, the simulator is a precisely controlled environment, in which flow conditions and geometry can be systematically prescribed and resultant valvular function and hemodynamics assessed. Thus, this work represents the first comprehensive database of high fidelity experimental data, critical for extensive validation of mitral valve fluid structure interaction simulations.

Left ventricular vortex formation is unaffected by diastolic impairment

Normal left ventricular (LV) filling occurs rapidly early in diastole caused by a progressive pressure gradient within the ventricle and with a low left atrial pressure. This normal diastolic function is altered in patients with heart failure. Such impairment of diastolic filling is manifested as an abrupt deceleration of the early filling wave velocity. Although variations within the early filling wave have been observed previously, the underlying hydrodynamic mechanisms are not well understood. Previously, it was proposed that the mitral annulus vortex ring formation time was the total duration of early diastolic filling and provided a measure of the efficiency of diastolic filling. However, we found that the favorable LV pressure difference driving early diastolic filling becomes zero simultaneously with the deceleration of the early filling wave propagation velocity and pinch-off of the LV vortex ring. Thus we calculated the vortex ring formation time using the duration of the early diastolic filling wave from its initiation to the time of the early filling wave propagation velocity deceleration when pinch-off occurs. This formation time does not vary with decreasing intraventricular pressure difference or with degree of diastolic dysfunction. Thus we conclude the vortex ring pinch-off occurs before the completion of early diastole, and its formation time remains invariant to changes of diastolic function.

On the three-dimensional vortical structure of early diastolic flow in a patient-specific left ventricle

We study the formation of the mitral vortex ring during early diastolic filling in a patient-specific left ventricle using direct numerical simulation. The geometry of the left ventricle is reconstructed from Magnetic Resonance Imaging (MRI). The heart wall motion is modeled by a cell-based activation methodology, which yields physiologic kinematics with heart rate equal to 52 beats per minute. We show that the structure of the mitral vortex ring consists of the main vortex ring and trailing vortex tubes, which originate at the heart wall. The trailing vortex tubes play an important role in exciting twisting circumferential instability modes of the mitral vortex ring. At the end of diastole, the vortex ring impinges on the wall and the intraventricular flow transitions to a weak turbulent state. Our results can be used to help interprete and analyze three-dimensional in-vivo flow measurements obtained with MRI.

Novel methods for assessment of right heart structure and function in pulmonary hypertension

Long-term increases in pulmonary vascular resistance and pulmonary arterial pressure resulting from structural alterations and abnormal vasoreactivity of the pulmonary vasculature may lead to right ventricular (RV) remodeling. Conventional methods of assessment of RV structure and function do not provide sensitive markers of RV remodeling for prognostic information. Advances in cardiac imaging have provided the capability to obtain quantitative information on the RV structure and function. This article reviews the clinical conditions that result in PH and discusses the novel and emerging methods for the assessment of right heart structure and function in PH in infants and children.

Assessment of left ventricular hemodynamics and function of patients with uremia by vortex formation using vector flow mapping

A novel echocardiographic method, vector flow mapping (VFM), acquires velocity vector from color Doppler velocity data. The purpose of this study was to evaluate whether VFM could provide useful information on intracardiac flow and helpful to evaluate left ventricular (LV) function. Thirty-eight patients with uremia undergoing hemodialysis and 30 healthy volunteers were enrolled. The maximum vector velocity, maximum diameter and duration of the intracardiac vortex were measured using VFM software during systole and diastole. The maximum vector velocity of the vortex and the peak velocities at the basal septum and lateral mitral annulus measured by tissue Doppler imaging (TDI) were correlated. The maximum diameter and duration of vortex formation were significantly higher in uremic patients compared with the control group during the ejection phase (40.6 ± 7.9 cm/sec vs. 28.1 ± 3.9 cm/sec; 297.1 ± 22.1 msec vs. 145.4 ± 19.3 msec, all P < 0.001). The maximal diameters of the vortex were higher in uremic patients compared with the control group during diastole (25.6 ± 3.4 mm vs. 16.4 ± 2.1 mm; 34.3 ± 3.1 mm vs. 26.8 ± 3.9 mm; 37.5 ± 2.4 mm vs. 20.9 ± 2.1 mm; all P < 0.001). The maximum vector velocities were lower in mid-diastole and late diastole (23.6 ± 2.3 cm/sec vs. 45.2 ± 3.7 cm/sec; 31.9 ± 2.9 cm/sec vs. 54.7 ± 3.2 cm/sec, all P < 0.001). There was a correlation between the maximum vector velocity of the vortex in mid-diastole and E'/A' at the septum and lateral mitral annulus (r = 0.70, r = 0.76, P < 0.001). Vortex can be utilized to provide intracardiac dynamic information using VFM and it may be a good supplement for evaluating LV function.

Can echocardiographic particle image velocimetry correctly detect motion patterns as they occur in blood inside heart chambers? A validation study using moving phantoms

Aims

To validate Echo Particle Image Velocimetry (PIV)

Methods

High fidelity string and rotating phantoms moving with different speed patterns were imaged with different high-end ultrasound systems at varying insonation angles and frame rates. Images were analyzed for velocity and direction and for complex motion patterns of blood flow with dedicated software. Post-processing was done with MATLAB-based tools (Dflow, JUV, University Leuven).

Results

Velocity estimation was accurate up to a velocity of 42 cm/s (r = 0.99, p < 0.001, mean difference 0.4 ± 2 cm/s). Maximally detectable velocity, however, was strongly dependent on frame rate and insonation angle and reached 42 cm/s under optimal conditions. At higher velocities estimates became random. Direction estimates did depend less on velocity and were accurate in 80-90%. In-plane motion patterns were correctly identified with three ultrasound systems.

Conclusion

Echo-PIV appears feasible. Velocity estimates are accurate, but the maximal detectable velocity depends strongly on acquisition parameters. Direction estimation works sufficiently, even at higher velocities. Echo-PIV appears to be a promising technical approach to investigate flow patterns by echocardiography.

Medical image diagnostics based on computer-aided flow analysis using magnetic resonance images

Most of the cardiac abnormalities have an implication on hemodynamics and affect cardiovascular health. Diagnostic imaging modalities such as computed tomography and magnetic resonance imaging provide excellent anatomical information on myocardial structures, but fail to show the cardiac flow and detect heart defects in vivo condition. The computerized technique for fluid motion estimation by pixel intensity tracking based on magnetic resonance signals represents a promising technique for functional assessment of cardiovascular disease, as it can provide functional information of the heart in addition to analysis of its anatomy. Cardiovascular flow characteristics can be measured in both normal controls and patients with cardiac abnormalities such as atrial septal defect, thus, enabling identification of the underlying causes of these flow phenomena. This review paper focuses on an overview of a flow analysis scheme based on computer-aided evaluation of magnetic resonance intensity images, in comparison with other commonly used medical imaging modalities. Details of the proposed technique are provided with validations being conducted at selected abnormal cardiovascular patients. It is expected that this new technique can potentially extend applications for characterizing cardiovascular defects and their hemodynamic behavior.

Describing the highly three dimensional Right Ventricle flow

Visualization of the three-dimensional flow within the Right Ventricle (RV) is a challenging issue due to the fully three-dimensional geometry of the ventricular cavity. To date proper characterization and quantification of the RV flow still remains incomplete, and techniques that can be easily applied to current medical imaging are scarce. A method for simulating the highly complex, multi directional flow within the RV is presented by coupling 4D echocardiography imaging with numerical simulations based on the Immersed Boundaries Method (IBM). A novel formulation for accurately computing the space-time distribution of the blood residence time inside the cavity is introduced. Results showed an initial compact vortex forming past the tricuspid orifice at early diastole that quickly breaks into a weakly turbulent flow pattern and rearranges, during systole, into a peculiar stream-wise vortex spinning out towards the pulmonary orifice. This arrangement is maintained when the Ejection Fraction (EF) is reduced from 58 to 32%. The average blood transit time is found to scale almost inversely proportional to the EF. A careful analysis of the residence time permitted to assess the relative significance of the different flow components (from the direct flow, with a residence time less than one heartbeat, to the residual volume, that stagnates in the ventricle) and their distribution in space.

In vivo measurements of blood flow in a rat using X-ray imaging technique

To measure instantaneous velocity fields of venous blood flow in a rat using X-ray particle tracking method. Gold nanoparticles (AuNPs) incorporated chitosan microparticles were applied as biocompatible flow tracers. After intravenous injection of the AuNP-chitosan particles into 7- to 9-week-old male rat vein, X-ray images of particle movement inside the cranial vena cava were consecutively captured. Individual AuNP-chitosan particles in the venous blood flow were clearly observed, and the corresponding velocity vectors were successfully extracted. The measured velocity vectors are in good agreement with the theoretical velocity profile suggested by Casson. This is the first trial to measure blood flow in animals under in vivo conditions with X-ray imaging technique. The results show that X-ray particle tracking technique has a great potential for in vivo measurements of blood flow, which can extend to various biomedical applications related with the diagnosis of circulatory vascular diseases.

Visualization of the intracavitary blood flow in systemic ventricles of Fontan patients by contrast echocardiography using particle image velocimetry

BACKGROUND

Flow patterns in univentricular hearts may have clinical value. Therefore, it is our objective to asses and characterize vortex flow patterns with Fontan circulation in comparison with healthy controls.

METHODS

Twenty-three patients (8 Fontan and 15 normal patients) underwent echocardiography with intravenous contrast agent (Sonovue®) administration. Dedicated software was used to perform particle image velocimetry (PIV) and to visualize intracavitary flow in the systemic ventricles of the patients. Vortex parameters including vortex depth, length, width, and sphericity index were measured. Vortex pulsatility parameters including relative strength, vortex relative strength, and vortex pulsation correlation were also measured.

RESULTS

The data from this study show that it is feasible to perform particle velocimetry in Fontan patients. Vortex length (VL) was significantly lower (0.51 ± 0.09 vs 0.65 ± 0.12, P = 0.010) and vortex width (VW) (0.32 ± 0.06 vs 0.27 ± 0.04, p = 0.014), vortex pulsation correlation (VPC) (0.26 ± 0.25 vs -0.22 ± 0.87, p = 0.05) were significantly higher in Fontan patients. Sphericity index (SI) (1.66 ± 0.48 vs 2.42 ± 0.62, p = 0.005), relative strength (RS) (0.77 ± 0.33 vs 1.90 ± 0.47, p = 0.0001), vortex relative strength (VRS) (0.18 ± 0.13 vs 0.43 ± 0.14, p = 0.0001) were significantly lower in the Fontan patients group.

CONCLUSIONS

PIV using contrast echocardiography is feasible in Fontan patients. Fontan patients had aberrant flow patterns as compared to normal hearts in terms of position, shape and sphericity of the main vortices. The vortex from the Fontan group was consistently shorter, wider and rounder than in controls. Whether vortex characteristics are related with clinical outcome is subject to further investigation.

The left ventricular intracavitary vortex during the isovolumic contraction period as detected by vector flow mapping

AIMS

The purpose of this study was to characterize left ventricular (LV) intracavitary flow during the isovolumic contraction (IVC) period in humans using vector flow mapping.

METHODS

Color flow Doppler imaging was performed from the apical long-axis view in 61 patients with heart failure and 58 healthy volunteers. Doppler flow data obtained during IVC were analyzed offline with vector flow mapping.

RESULTS

A large vortex was formed from the LV inflow toward the outflow during IVC. In normal subjects, the area of the vortex was sustained, but the flow volume decreased significantly during IVC (P < 0.001). A significant apex-to-base flow velocity gradient was shown along the outflow axis on aortic valve opening. However, both the area and flow volume of the vortex decreased more severely during IVC in the patients (P < 0.001). The apex-to-base flow velocity gradient along the outflow axis disappeared and a reversed velocity gradient was observed at the basal-mid level on aortic valve opening. In multivariate models, a decreased LV ejection fraction was the only independent predictor of the percentage decrease in area of the vortex during the IVC (P < 0.001), and a larger QRS width (P = 0.028) and LV end-systolic long diameter (P = 0.002) were independent predictors of the percentage decrease in flow volume of the vortex.

CONCLUSIONS

The vortex across the LV inflow-outflow region during IVC facilitates the ejection of blood during early systole, and an unsustained vortex may be associated with impaired cardiac function.

Imaging of wall motion coupled with blood flow velocity in the heart and vessels in vivo: a feasibility study

The mechanical property and geometry changes as a result of cardiovascular disease affect both the wall motion and blood flow in the heart and vessels, whereas the latter two are also coupled and therefore continuously influence one another. Simultaneous and registered imaging of both cardiovascular wall motion and blood velocity may thus contribute to more complete computational models of cardiovascular mechanical and fluid dynamics as well as provide additional diagnostic information. The objective of this paper was to determine the feasibility of imaging cardiovascular wall motion coupled with blood flow in vivo. Normal (n = 6) and infarcted (n = 5) murine left ventricles, and normal (n = 5) and aneurysmal (n = 4) murine abdominal aortas, were imaged in longitudinal views with a 30-MHz ultrasound probe. Using electrocardiogram (ECG) gating, 2-D radio-frequency (RF) data were acquired at a frame rate of 8 kHz. The axial wall velocity and blood velocity were estimated using a speckle-tracking technique. Spatially and temporally registered imaging of both cardiovascular wall motion and blood flow was shown to be feasible. Reduced wall motion was detected in the infarcted region, whereas vortex flow patterns were imaged in diastolic phases of both normal and infarcted left ventricles. The myocardial wall motion and blood flow were found to be more synchronous in the normal heart, where the blood moves toward the anteroseptal wall after the mitral valve opens (i.e., rapid filling phase), and the anteroseptal wall simultaneously undergoes outward motion. In the infarcted heart, however, in the rapid filling phase, the basal anteroseptal wall starts moving about 20 ms before the mitral valve opens and the blood enters the left ventricle. In the normal aorta, the wall motion and blood velocity were uniform and synchronous. In the aneurysmal aorta, reduced and spatially varied wall motion and vortex flow patterns in the aneurysmal sac were found. The wall motion and blood velocity were thus less synchronous in the aneurysmal aorta. Cardiovascular wall motion and blood flow were both imaged in mice in vivo. This dual information may provide important insights for the diagnosis of cardiovascular disease as well as essential parameters for its biomechanical modeling.

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.

In vitro and preliminary in vivo validation of echo particle image velocimetry in carotid vascular imaging

Noninvasive, easy-to-use and accurate measurements of wall shear stress (WSS) in human blood vessels have always been challenging in clinical applications. Echo particle image velocimetry (Echo PIV) has shown promise for clinical measurements of local hemodynamics and wall shear rate. Thus far, however, the method has only been validated under simple flow conditions. In this study, we validated Echo PIV under in vitro and in vivo conditions. For in vitro validation, we used an anatomically correct, compliant carotid bifurcation flow phantom with pulsatile flow conditions, using optical particle image velocimetry (optical PIV) as the reference standard. For in vivo validation, we compared Echo PIV-derived 2-D velocity fields obtained at the carotid bifurcation in five normal subjects against phase-contrast magnetic resonance imaging (PC-MRI)-derived velocity measurements obtained at the same locations. For both studies, time-dependent, 2-D, two-component velocity vectors; peak/centerline velocity, flow rate and wall shear rate (WSR) waveforms at the common carotid artery (CCA), carotid bifurcation and distal internal carotid artery (ICA) were examined. Linear regression, correlation analysis and Bland-Altman analysis were used to quantify the agreement of different waveforms measured by the two techniques. In vitro results showed that Echo PIV produced good images of time-dependent velocity vector maps over the cardiac cycle with excellent temporal (up to 0.7 ms) and spatial (∼0.5 mm) resolutions and quality, comparable with optical PIV results. Further, good agreement was found between Echo PIV and optical PIV results for velocity and WSR measurements. In vivo results also showed good agreement between Echo PIV velocities and phase contrast MRI velocities. We conclude that Echo PIV provides accurate velocity vector and WSR measurements in the carotid bifurcation and has significant potential as a clinical tool for cardiovascular hemodynamics evaluation.

Two-dimensional intraventricular flow mapping by digital processing conventional color-Doppler echocardiography images

Doppler echocardiography remains the most extended clinical modality for the evaluation of left ventricular (LV) function. Current Doppler ultrasound methods, however, are limited to the representation of a single flow velocity component. We thus developed a novel technique to construct 2D time-resolved (2D+t) LV velocity fields from conventional transthoracic clinical acquisitions. Combining color-Doppler velocities with LV wall positions, the cross-beam blood velocities were calculated using the continuity equation under a planar flow assumption. To validate the algorithm, 2D Doppler flow mapping and laser particle image velocimetry (PIV) measurements were carried out in an atrio-ventricular duplicator. Phase-contrast magnetic resonance (MR) acquisitions were used to measure in vivo the error due to the 2D flow assumption and to potential scan-plane misalignment. Finally, the applicability of the Doppler technique was tested in the clinical setting. In vitro experiments demonstrated that the new method yields an accurate quantitative description of the main vortex that forms during the cardiac cycle (mean error for vortex radius, position and circulation). MR image analysis evidenced that the error due to the planar flow assumption is close to 15% and does not preclude the characterization of major vortex properties neither in the normal nor in the dilated LV. These results are yet to be confirmed by a head-to-head clinical validation study. Clinical Doppler studies showed that the method is readily applicable and that a single large anterograde vortex develops in the healthy ventricle while supplementary retrograde swirling structures may appear in the diseased heart. The proposed echocardiographic method based on the continuity equation is fast, clinically-compliant and does not require complex training. This technique will potentially enable investigators to study of additional quantitative aspects of intraventricular flow dynamics in the clinical setting by high-throughput processing conventional color-Doppler images.

Modeling radial viscoelastic behavior of left ventricle based on MRI tissue phase mapping

The viscoelastic behavior of myocardial tissue is a measure that has recently found to be a deterministic factor in quality of contraction. Parameters imposing the viscoelastic behavior of the heart are influenced in part by sarcomere function and myocardial composition. Despite the overall agreement on significance of cardiac viscoelasticity, a practical model that can measure and characterize the viscoelastic behavior of the myocardial segments does not yet exist. Pressure-Volume (P-V) curves are currently the only measure for stiffness/compliance of the left ventricle. However, obtaining P-V curves requires invasive cardiac catheterization, and only provides qualitative information on how pressure and volume change with respect to each other. For accurate assessment of myocardial mechanical behavior, it is required to obtain quantitative measures for viscoelasticity. In this work, we have devised a model that yields myocardial elastic and viscous damping coefficient functions through the cardiac cycle. The required inputs for this model are kinematic information with respect to changes in LV short axes that were obtained by Magnetic Resonance Imaging (MRI) using a tissue phase mapping (TPM) pulse sequence. We evaluated viscoelastic coefficients of LV myocardium in two different age groups of 20-40 and greater than 60. We found that the magnitude of stiffness coefficients is noticeably greater in the older subjects. Additionally, we found that slope of viscous damping functions follow similar patterns for each individual age group. This method may shed light on dynamics of contraction through MRI in conditions where composition of myocardium is changed such as in aging, adverse remodeling, and cardiomyopathies.