Image-Based Computational Fluid Dynamics in Blood Vessel Models: Toward Developing a Prognostic Tool to Assess Cardiovascular Function Changes in Prolonged Space Flights

One of NASA's objectives is to be able to perform a complete, pre-flight, evaluation of cardiovascular changes in astronauts scheduled for prolonged space missions. Computational fluid dynamics (CFD) has shown promise as a method for estimating cardiovascular function during reduced gravity conditions. For this purpose, MRI can provide geometrical information, to reconstruct vessel geometries, and measure all spatial velocity components, providing location specific boundary conditions. The objective of this study was to investigate the reliability of MRI-based model reconstruction and measured boundary conditions for CFD simulations. An aortic arch model and a carotid bifurcation model were scanned in a 1.5T Siemens MRI scanner. Axial MRI acquisitions provided images for geometry reconstruction (slice thickness 3 and 5 mm; pixel size 1x1 and 0.5x0.5 square millimeters). Velocity acquisitions provided measured inlet boundary conditions and localized three-directional steady-flow velocity data (0.7-3.0 L/min). The vessel walls were isolated using NIH provided software (ImageJ) and lofted to form the geometric surface. Constructed and idealized geometries were imported into a commercial CFD code for meshing and simulation. Contour and vector plots of the velocity showed identical features between the MRI velocity data, the MRI-based CFD data, and the idealized-geometry CFD data, with less than 10% differences in the local velocity values. CFD results on models reconstructed from different MRI resolution settings showed insignificant differences (less than 5%). This study illustrated, quantitatively, that reliable CFD simulations can be performed with MRI reconstructed models and gives evidence that a future, subject-specific, computational evaluation of the cardiovascular system alteration during space travel is feasible

Relationship between the extent of non-viable myocardium and regional left ventricular function in chronic ischemic heart disease

Purpose. To define the relationship between left ventricular (LV) regional contractile function and the extent of myocardial scar in patients with chronic ischemic heart disease and multi-vessel coronary artery disease. Methods. Twenty-three patients with chronic ischemic heart disease and 5 healthy volunteers underwent magnetic resonance imaging (MRI). In patients, the relative area ( Percent Scar) and transmural extent (Transmurality) of myocardial infarction were computed from short-axis delayed enhancement images. In each image, myocardial segments were categorized based on the extent of infarction they contained, with 6 categories each for Percent Scar and Transmurality: normal, from healthy volunteers; and 0%; 1–25%, 26–50%, 51–75%, and \u3e 76% from patients. In patients and volunteers, regional LV function was quantified by absolute systolic wall thickening from cine images and midwall circumferential strain using tagged images. Results. Compared to normal segments, regional LV function in patients was significantly diminished in all scar extent intervals, with wall thickening=-8% for all categories. Systolic wall thickening was reduced significantly in all categories above 50% Percent Scar and above 25% Transmurality in patients, relative to corresponding 0% categories. Circumferential strain was significantly reduced above 25% Percent Scar and above 25% Transmurality. Conclusions. In patients with chronic ischemic heart disease and multivessel coronary artery disease, wall thickening was more sensitive to changes in scar Transmurality than to changes in Percent Scar. However, circumferential strain was equally sensitive to both indices. In general, circumferential strain was more sensitive than wall thickening to increases in scar extent

Segmentation of Non-viable Myocardium in Delayed Enhancement Magnetic Resonance Images

Purpose: To evaluate six algorithms for segmenting non-viable left ventricular (LV) myocardium in delayed enhancement (DE) magnetic resonance imaging (MRI). Methods: Twenty-three patients with known chronic ischemic heart disease underwent DE-MRI. DE images were first manually thresholded using an interactive region-filling tool to isolate non-viable myocardium. Then, six thresholding algorithms, based on the image intensity characteristics of either LV blood pool (BP), viable LV myocardium, or both, were applied to each image. For the Mean−2SDBP algorithm, thresholds were equal to the mean BP intensity minus twice its standard deviation. For the Mean+2SDSemi, Mean+3SDSemi, Mean+2SDAuto, and Mean+3SDAuto algorithms, thresholds equaled the mean intensity of viable myocardium plus twice (or thrice, as denoted by the name) the standard deviation of intensity (subscripts denote how these values were determined: automatic or semi-automatic). For the Minimum Intensity algorithm, the threshold equaled the minimum intensity between the BP and LV myocardium mean intensities. Percent Scar was defined as the ratio of non-viable to total myocardial pixels in each image. Agreement between each algorithm and manual thresholding was assessed using Bland–Altman analysis. Results: Mean Percent Scar was 25 ± 16% by manual thresholding. Five of the six algorithms demonstrated mean bias within ±3% (all except Mean+2SDAuto); however, limits of agreement (LoA) were large in general (range 12–36%). The best overall agreement was demonstrated by the Mean+2SDSemi (bias, 0%; LoA, 12%) and Mean+3SDSemi(bias, −3%; LoA, 14%) algorithms. Conclusion: On average, five of the six algorithms proved satisfactory for clinical implementation; however, in some images, manual correction of automatic results was necessary

Relationship between the extent of non-viable myocardium and regional left ventricular function in chronic ischemic heart disease

Purpose. To define the relationship between left ventricular (LV) regional contractile function and the extent of myocardial scar in patients with chronic ischemic heart disease and multi-vessel coronary artery disease. Methods. Twenty-three patients with chronic ischemic heart disease and 5 healthy volunteers underwent magnetic resonance imaging (MRI). In patients, the relative area ( Percent Scar) and transmural extent (Transmurality) of myocardial infarction were computed from short-axis delayed enhancement images. In each image, myocardial segments were categorized based on the extent of infarction they contained, with 6 categories each for Percent Scar and Transmurality: normal, from healthy volunteers; and 0%; 1–25%, 26–50%, 51–75%, and \u3e 76% from patients. In patients and volunteers, regional LV function was quantified by absolute systolic wall thickening from cine images and midwall circumferential strain using tagged images. Results. Compared to normal segments, regional LV function in patients was significantly diminished in all scar extent intervals, with wall thickening=-8% for all categories. Systolic wall thickening was reduced significantly in all categories above 50% Percent Scar and above 25% Transmurality in patients, relative to corresponding 0% categories. Circumferential strain was significantly reduced above 25% Percent Scar and above 25% Transmurality. Conclusions. In patients with chronic ischemic heart disease and multivessel coronary artery disease, wall thickening was more sensitive to changes in scar Transmurality than to changes in Percent Scar. However, circumferential strain was equally sensitive to both indices. In general, circumferential strain was more sensitive than wall thickening to increases in scar extent

Noninvasive Quantification of Fluid Mechanical Energy Losses in the Total Cavopulmonary Connection with Magnetic Resonance Phase Velocity Mapping

A major determinant of the success of surgical vascular modifications, such as the total cavopulmonary connection (TCPC), is the energetic efficiency that is assessed by calculating the mechanical energy loss of blood flow through the new connection. Currently, however, to determine the energy loss, invasive pressure measurements are necessary. Therefore, this study evaluated the feasibility of the viscous dissipation (VD) method, which has the potential to provide the energy loss without the need for invasive pressure measurements. Two experimental phantoms, a U-shaped tube and a glass TCPC, were scanned in a magnetic resonance (MR) imaging scanner and the images were used to construct computational models of both geometries. MR phase velocity mapping (PVM) acquisitions of all three spatial components of the fluid velocity were made in both phantoms and the VD was calculated. VD results from MR PVM experiments were compared with VD results from computational fluid dynamics (CFD) simulations on the image-based computational models. The results showed an overall agreement between MR PVM and CFD. There was a similar ascending tendency in the VD values as the image spatial resolution increased. The most accurate computations of the energy loss were achieved for a CFD grid density that was too high for MR to achieve under current MR system capabilities (in-plane pixel size of less than 0.4 mm). Nevertheless, the agreement between the MR PVM and the CFD VD results under the same resolution settings suggests that the VD method implemented with a clinical imaging modality such as MR has good potential to quantify the energy loss in vascular geometries such as the TCPC

Reliable In-Plane Velocity Measurements With Magnetic Resonance Velocity Imaging

Magnetic resonance (MR) imaging is a well-known diagnostic imaging modality. In addition to its high-quality imaging capabilities, hydrogen-based MR can also provide non-invasively the velocity of water-based fluids in all three spatial directions (through-plane and in-plane) in an image. Many previous studies showed that MR velocity imaging can accurately measure the through-plane velocity. The aim of this study was to evaluate how reliable are the in-plane velocity measurements in an image. The axial velocity of water in horizontal tubes (inner diameter: 14.7–26.2 mm) was measured with segmented (fast) and non-segmented (slow) k-space MR velocity imaging using: (a) an imaging slice placed perpendicular to the tube axis with through-plane velocity-encoding; and (b) an imaging slice placed parallel to the tube axis with in-plane velocity-encoding. The two planes intersected along the vertical tube-centerline. The flow rate was accurately quantified (mean error plane velocity profiles were not significantly different from the through-plane profiles (mean difference =6%, correlation coefficients \u3e0.98). There was no significant difference between the velocity profiles from the segmented and the non-segmented sequences (mean difference 0.95). The results of this study suggest that fast MR velocity imaging can measure the in-plane velocity in an image with reliability

Extent of Left Ventricular Scar Predicts Outcomes in Ischemic Cardiomyopathy Patients With Significantly Reduced Systolic Function A Delayed Hyperenhancement Cardiac Magnetic Resonance Study

ObjectivesThe objective of the study was to determine whether the extent of left ventricular scar, measured with delayed hyperenhancement cardiac magnetic resonance (DHE-CMR), predicts survival in patients with ischemic cardiomyopathy (ICM) and severely reduced left ventricular ejection fraction (LVEF).BackgroundPatients with ICM and reduced LVEF have poor survival. Such patients have a high myocardial scar burden. CMR is highly accurate in delineation of myocardial scar.MethodsWe studied 349 patients (76% men) with severe ICM (≥70% disease in ≥1 epicardial coronary, and mean LVEF of 24%) that underwent DHE-CMR (Siemens 1.5-T scanner, Erlangen, Germany), between 2003 and 2006. Scar (quantified as percentage of myocardium) was defined on DHE-MR images as an intensity >2 standard deviations above the viable myocardium. Transmurality score was semiquantitatively recorded in a 17-segment model as: 0 = no scar, 1 = 1% to 25% scar, 2 = 26% to 50%, 3 = 51% to 75%, and 4 = >75%. The LVEF, demographic data, risk factors, need for cardiac transplantation (CTx), and all-cause mortality were recorded.ResultsThe mean age and follow-up were 65 ± 11 years and 2.6 ± 1.2 years (median 2.4 years [1.1, 3.5]), respectively. There were 56 events (51 deaths and 5 CTx). Mean scar percentage and transmurality score were higher in patients with events versus those without (39 ± 22 vs. 30 ± 20, p = 0.003, and 9.7 ± 5 vs. 7.8 ± 5, p = 0.004). On Cox proportional hazard survival analysis, quantified scar was greater than the median (30% of total myocardium), and female gender predicted events (relative risk 1.75 [95% Confidence Interval: 1.02 to 3.03] and relative risk 1.83 [95% Confidence Interval: 1.06 to 3.16], respectively, both p = 0.03).ConclusionsIn patients with ICM and severely reduced LVEF, a greater extent of myocardial scar, delineated by DHE-CMR is associated with increased mortality or the need for cardiac transplantation, potentially aiding further risk-stratification