Three-dimensional reconstruction of the flow in a human left heart by using magnetic resonance phase velocity encoding

CT-Based Simulation of Left Ventricular Hemodynamics: A Pilot Study in Mitral Regurgitation and Left Ventricle Aneurysm Patients

Background

Cardiac CT (CCT) is well suited for a detailed analysis of heart structures due to its high spatial resolution, but in contrast to MRI and echocardiography, CCT does not allow an assessment of intracardiac flow. Computational fluid dynamics (CFD) can complement this shortcoming. It enables the computation of hemodynamics at a high spatio-temporal resolution based on medical images. The aim of this proposed study is to establish a CCT-based CFD methodology for the analysis of left ventricle (LV) hemodynamics and to assess the usability of the computational framework for clinical practice.

Materials and Methods

The methodology is demonstrated by means of four cases selected from a cohort of 125 multiphase CCT examinations of heart failure patients. These cases represent subcohorts of patients with and without LV aneurysm and with severe and no mitral regurgitation (MR). All selected LVs are dilated and characterized by a reduced ejection fraction (EF). End-diastolic and end-systolic image data was used to reconstruct LV geometries with 2D valves as well as the ventricular movement. The intraventricular hemodynamics were computed with a prescribed-motion CFD approach and evaluated in terms of large-scale flow patterns, energetic behavior, and intraventricular washout.

Results

In the MR patients, a disrupted E-wave jet, a fragmentary diastolic vortex formation and an increased specific energy dissipation in systole are observed. In all cases, regions with an impaired washout are visible. The results furthermore indicate that considering several cycles might provide a more detailed view of the washout process. The pre-processing times and computational expenses are in reach of clinical feasibility.

Conclusion

The proposed CCT-based CFD method allows to compute patient-specific intraventricular hemodynamics and thus complements the informative value of CCT. The method can be applied to any CCT data of common quality and represents a fair balance between model accuracy and overall expenses. With further model enhancements, the computational framework has the potential to be embedded in clinical routine workflows, to support clinical decision making and treatment planning.

The evolution of intraventricular vortex during ejection studied by using vector flow mapping

AIMS

The purpose of this study was to assess the evolution of intraventricular vortex during left ventricular (LV) ejection.

METHODS

Vector flow mapping was performed in 51 patients with coronary artery disease and LV ejection fraction (EF) >50%, 70 patients with EF <50% (13 with coronary artery disease and 57 with dilated cardiomyopathy), and 62 healthy volunteers.

RESULTS

In normals and patients with EF >50%, the intraventricular vortex dissipated quickly during early ejection. In patients with EF <50%, the vortex stayed mainly at apex and persisted for a significantly longer time. The evolution of vortex during ejection was significantly correlated with QRS width, EF, fractional shortening, LV outflow velocity time integral, wall motion score index (WMSI), LV dimensions, left atrial diameter, and diastolic mitral annular velocities. LV end-diastolic short diameter and WMSI were the independent determinants of the duration of vortex (R(2) = 0.482, P < 0.001). End-systolic short diameter and apical WMSI were the independent determinants of duration of vortex corrected for ejection time (R(2) = 0.565, P < 0.001). End-systolic short diameter was the independent determinant of percentage change in vortex area during early ejection (R(2) = 0.355, P < 0.001). End-systolic short diameter and ejection time were the independent determinants of percentage change in vortex flow volume (R(2) = 0.415, P < 0.001).

CONCLUSIONS

In patients with LV systolic dysfunction, the vortex persists during ejection and stays mainly at apex. The vortex evolution during ejection is closely associated with LV dimensions and functions.

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.

Semi-automatic quantification of 4D left ventricular blood flow

BACKGROUND

The beating heart is the generator of blood flow through the cardiovascular system. Within the heart's own chambers, normal complex blood flow patterns can be disturbed by diseases. Methods for the quantification of intra-cardiac blood flow, with its 4D (3D+time) nature, are lacking. We sought to develop and validate a novel semi-automatic analysis approach that integrates flow and morphological data.

METHOD

In six healthy subjects and three patients with dilated cardiomyopathy, three-directional, three-dimensional cine phase-contrast cardiovascular magnetic resonance (CMR) velocity data and balanced steady-state free-precession long- and short-axis images were acquired. The LV endocardium was segmented from the short-axis images at the times of isovolumetric contraction (IVC) and isovolumetric relaxation (IVR). At the time of IVC, pathlines were emitted from the IVC LV blood volume and traced forwards and backwards in time until IVR, thus including the entire cardiac cycle. The IVR volume was used to determine if and where the pathlines left the LV. This information was used to automatically separate the pathlines into four different components of flow: Direct Flow, Retained Inflow, Delayed Ejection Flow and Residual Volume. Blood volumes were calculated for every component by multiplying the number of pathlines with the blood volume represented by each pathline. The accuracy and inter- and intra-observer reproducibility of the approach were evaluated by analyzing volumes of LV inflow and outflow, the four flow components, and the end-diastolic volume.

RESULTS

The volume and distribution of the LV flow components were determined in all subjects. The calculated LV outflow volumes [ml] (67 +/- 13) appeared to fall in between those obtained by through-plane phase-contrast CMR (77 +/- 16) and Doppler ultrasound (58 +/- 10), respectively. Calculated volumes of LV inflow (68 +/- 11) and outflow (67 +/- 13) were well matched (NS). Low inter- and intra-observer variability for the assessment of the volumes of the flow components was obtained.

CONCLUSIONS

This semi-automatic analysis approach for the quantification of 4D blood flow resulted in accurate LV inflow and outflow volumes and a high reproducibility for the assessment of LV flow components.

Three dimensional flow in the human left atrium

BACKGROUND

Abnormal flow patterns in the left atrium in atrial fibrillation or mitral stenosis are associated with an increased risk of thrombosis and systemic embolisation; the characteristics of normal atrial flow that avoid stasis have not been well defined.

OBJECTIVES

To present a three dimensional particle trace visualisation of normal left atrial flow in vivo, constructed from flow velocities in three dimensional space.

METHODS

Particle trace visualisation of time resolved three dimensional magnetic resonance imaging velocity measurements was used to provide a display of intracardiac flow without the limitations of angle sensitivity or restriction to imaging planes. Global flow patterns of the left atrium were studied in 11 healthy volunteers.

RESULTS

In all subjects vortical flow was observed in the atrium during systole and diastolic diastasis (mean (SD) duration of systolic vortex, 280 (77) ms; and of diastolic vortex, 256 (118) ms). The volume incorporated and recirculated within the vortices originated predominantly from the left pulmonary veins. Inflow from the right veins passed along the vortex periphery, constrained between the vortex and the atrial wall.

CONCLUSIONS

Global left atrial flow in the normal human heart comprises consistent patterns specific to the phase of the cardiac cycle. Separate paths of left and right pulmonary venous inflow and vortex formation may have beneficial effects in avoiding left atrial stasis in the normal subject in sinus rhythm.

Noninvasive estimation of transmitral pressure drop across the normal mitral valve in humans: importance of convective and inertial forces during left ventricular filling

OBJECTIVES

We hypothesized that color M-mode (CMM) images could be used to solve the Euler equation, yielding regional pressure gradients along the scanline, which could then be integrated to yield the unsteady Bernoulli equation and estimate noninvasively both the convective and inertial components of the transmitral pressure difference.

BACKGROUND

Pulsed and continuous wave Doppler velocity measurements are routinely used clinically to assess severity of stenotic and regurgitant valves. However, only the convective component of the pressure gradient is measured, thereby neglecting the contribution of inertial forces, which may be significant, particularly for nonstenotic valves. Color M-mode provides a spatiotemporal representation of flow across the mitral valve.

METHODS

In eight patients undergoing coronary artery bypass grafting, high-fidelity left atrial and ventricular pressure measurements were obtained synchronously with transmitral CMM digital recordings. The instantaneous diastolic transmitral pressure difference was computed from the M-mode spatiotemporal velocity distribution using the unsteady flow form of the Bernoulli equation and was compared to the catheter measurements.

RESULTS

From 56 beats in 16 hemodynamic stages, inclusion of the inertial term ([deltapI]max = 1.78+/-1.30 mm Hg) in the noninvasive pressure difference calculation significantly increased the temporal correlation with catheter-based measurement (r = 0.35+/-0.24 vs. 0.81+/-0.15, p< 0.0001). It also allowed an accurate approximation of the peak pressure difference ([deltapc+I]max = 0.95 [delta(p)cathh]max + 0.24, r = 0.96, p<0.001, error = 0.08+/-0.54 mm Hg).

CONCLUSIONS

Inertial forces are significant components of the maximal pressure drop across the normal mitral valve. These can be accurately estimated noninvasively using CMM recordings of transmitral flow, which should improve the understanding of diastolic filling and function of the heart.

Three-dimensional computational model of left heart diastolic function with fluid-structure interaction

Aided by advancements in computer speed and modeling techniques, computational modeling of cardiac function has continued to develop over the past twenty years. The goal of the current study was to develop a computational model that provides blood-tissue interaction under physiologic flow conditions, and apply it to a thin-walled model of the left heart. To accomplish this goal, the Immersed Boundary Method was used to study the interaction of the tissue and blood in response to fluid forces and changes in tissue pathophysiology. The fluid mass and momentum conservation equations were solved using Patankar's Semi-Implicit Method for Pressure Linked Equations (SIMPLE). A left heart model was developed to examine diastolic function, and consisted of the left ventricle, left atrium, and pulmonary flow. The input functions for the model included the pulmonary driving pressure and time-dependent relationship for changes in chamber tissue properties during the simulation. The results obtained from the left heart model were compared to clinically observed diastolic flow conditions for validation. The inflow velocities through the mitral valve corresponded with clinical values (E-wave = 74.4 cm/s, A-wave = 43 cm/s, and E/A = 1.73). The pressure traces for the atrium and ventricle, and the appearance of the ventricular flow fields throughout filling, agreed with those observed in the heart. In addition, the atrial flow fields could be observed in this model and showed the conduit and pump functions that current theory suggests. The ability to examine atrial function in the present model is something not described previously in computational simulations of cardiac function.

Age-dependent changes in spatial and temporal blood velocity distribution of early left ventricular filling

This study describes early diastolic inflow dynamics based on three-directional magnetic resonance velocity data and investigates age-dependent changes in early diastolic inflow characteristics. We examined 26 young healthy volunteers age 25 (3) years (mean, SD), and 23 healthy older volunteers age 63 (8) years. Three-directional magnetic resonance velocity mapping was performed in a long axis plane through the heart. Transverse velocity profiles were read in five different positions in the early diastolic inflow stream of the left ventricle. The size and timing of the maximum velocities at each level were recorded and the repeatability of the method was tested. Compared with the younger group, the older group was characterized by: 1) lower maximum velocity in all positions, 2) increased deceleration of blood downstream from the mitral leaflet tips, and 3) delayed velocity propagation. The described method was repeatable and enabled detection of the age-dependent differences between groups of normal subjects. In conclusion, the early diastolic inflow pattern changes with age, probably reflecting changes in the diastolic function of the myocardium.

Particle trace visualization of intracardiac flow using time-resolved 3D phase contrast MRI

The flow patterns in the human heart are complex and difficult to visualize using conventional two-dimensional (2D) modalities, whether they depict a single velocity component (Doppler echocardiography) or all three components in a few slices (2D phase contrast MRI). To avoid these shortcomings, a temporally resolved 3D phase contrast technique was used to derive data describing the intracardiac velocity fields in normal volunteers. The MRI data were corrected for phase shifts caused by eddy currents and concomitant gradient fields, with improvement in the accuracy of subsequent flow visualizations. Pathlines describing the blood pathways through the heart were generated from the temporally resolved velocity data, starting from user-specified locations and time frames. Flow trajectories were displayed as 3D particle traces, with simultaneous demonstration of morphologic 2D slices. This type of visualization is intuitive and interactive and may extend our understanding of dynamic and previously unrecognized patterns of intracardiac flow.