Prahlad G. Menon

Variability of Computational Fluid Dynamics Solutions for Pressure and Flow in a Giant Aneurysm: The ASME 2012 Summer Bioengineering Conference CFD Challenge

Stimulated by a recent controversy regarding pressure drops predicted in a giant aneurysm with a proximal stenosis, the present study sought to assess variability in the prediction of pressures and flow by a wide variety of research groups. In phase I, lumen geometry, flow rates, and fluid properties were specified, leaving each research group to choose their solver, discretization, and solution strategies. Variability was assessed by having each group interpolate their results onto a standardized mesh and centerline. For phase II, a physical model of the geometry was constructed, from which pressure and flow rates were measured. Groups repeated their simulations using a geometry reconstructed from a micro-computed tomography (CT) scan of the physical model with the measured flow rates and fluid properties. Phase I results from 25 groups demonstrated remarkable consistency in the pressure patterns, with the majority predicting peak systolic pressure drops within 8% of each other. Aneurysm sac flow patterns were more variable with only a few groups reporting peak systolic flow instabilities owing to their use of high temporal resolutions. Variability for phase II was comparable, and the median predicted pressure drops were within a few millimeters of mercury of the measured values but only after accounting for submillimeter errors in the reconstruction of the life-sized flow model from micro-CT. In summary, pressure can be predicted with consistency by CFD across a wide range of solvers and solution strategies, but this may not hold true for specific flow patterns or derived quantities. Future challenges are needed and should focus on hemodynamic quantities thought to be of clinical interest.

Alk1 controls arterial endothelial cell migration in lumenized vessels.

Heterozygous loss of the arterial-specific TGFβ type I receptor, activin receptor-like kinase 1 (ALK1; ACVRL1), causes hereditary hemorrhagic telangiectasia (HHT). HHT is characterized by development of fragile, direct connections between arteries and veins, or arteriovenous malformations (AVMs). However, how decreased ALK1 signaling leads to AVMs is unknown. To understand the cellular mis-steps that cause AVMs, we assessed endothelial cell behavior in alk1-deficient zebrafish embryos, which develop cranial AVMs. Our data demonstrate that alk1 loss has no effect on arterial endothelial cell proliferation but alters arterial endothelial cell migration within lumenized vessels. In wild-type embryos, alk1-positive cranial arterial endothelial cells generally migrate towards the heart, against the direction of blood flow, with some cells incorporating into endocardium. In alk1-deficient embryos, migration against flow is dampened and migration in the direction of flow is enhanced. Altered migration results in decreased endothelial cell number in arterial segments proximal to the heart and increased endothelial cell number in arterial segments distal to the heart. We speculate that the consequent increase in distal arterial caliber and hemodynamic load precipitates the flow-dependent development of downstream AVMs. Summary: Cell tracking in a zebrafish model of hereditary hemorrhagic telangiectasia reveals that the primary defect leading to arteriovenous malformations is misdirected EC migration within lumenized arteries.

Characterization of neonatal aortic cannula jet flow regimes for improved cardiopulmonary bypass.

During pediatric and neonatal cardiopulmonary bypass (CPB), tiny aortic outflow cannulae (2-3 mm inner diameter), with micro-scale blood-wetting features transport relatively large blood volumes (0.3 to 1.0 L/min) resulting in high blood flow velocities (2 to 5 m/s). These severe flow conditions are likely to complement platelet activation, release pro-inflammatory cytokines, and further result in vascular and blood damage. Hemodynamically efficient aortic outflow cannulae are required to provide high blood volume flow rates at low exit force. In addition, optimal aortic insertion strategies are necessary in order to alleviate hemolytic risk, post-surgical neurological complications and developmental defects, by improving cerebral perfusion in the young patient. The methodology and results presented in this study serve as a baseline for design of superior aortic outflow cannulae. In this study, direct numerical simulation (DNS) computational fluid dynamics (CFD) was employed to delineate baseline hemodynamic performance of jet wakes emanating from microCT scanned state-of-the-art pediatric cannula tips in a cuboidal test rig operating at physiologically relevant laminar and turbulent Reynolds numbers (Re: 650-2150 , steady inflow). Qualitative and quantitative validation of CFD simulated device-specific jet wakes was established using time-resolved flow visualization and particle image velocimetry (PIV). For the standard end-hole cannula tip design, blood damage indices were further numerically assessed in a subject-specific cross-clamped neonatal aorta model for different cannula insertion configurations. Based on these results, a novel diffuser type cannula tip is proposed for improved jet flow-control, decreased blood damage and exit force and increased permissible flow rates. This study also suggests that surgically relevant cannula orientation parameters such as outflow angle and insertion depth may be important for improved hemodynamic performance. The jet flow design paradigm demonstrated in this study represents a philosophical shift towards cannula flow control enabling favorable pressure-drop versus outflow rate characteristics.

Effects of Intraluminal Thrombus on Patient-Specific Abdominal Aortic Aneurysm Hemodynamics via Stereoscopic Particle Image Velocity and Computational Fluid Dynamics Modeling

The pathology of the human abdominal aortic aneurysm (AAA) and its relationship to the later complication of intraluminal thrombus (ILT) formation remains unclear. The hemodynamics in the diseased abdominal aorta are hypothesized to be a key contributor to the formation and growth of ILT. The objective of this investigation is to establish a reliable 3D flow visualization method with corresponding validation tests with high confidence in order to provide insight into the basic hemodynamic features for a better understanding of hemodynamics in AAA pathology and seek potential treatment for AAA diseases. A stereoscopic particle image velocity (PIV) experiment was conducted using transparent patient-specific experimental AAA models (with and without ILT) at three axial planes. Results show that before ILT formation, a 3D vortex was generated in the AAA phantom. This geometry-related vortex was not observed after the formation of ILT, indicating its possible role in the subsequent appearance of ILT in this patient. It may indicate that a longer residence time of recirculated blood flow in the aortic lumen due to this vortex caused sufficient shear-induced platelet activation to develop ILT and maintain uniform flow conditions. Additionally, two computational fluid dynamics (CFD) modeling codes (Fluent and an in-house cardiovascular CFD code) were compared with the two-dimensional, three-component velocity stereoscopic PIV data. Results showed that correlation coefficients of the out-of-plane velocity data between PIV and both CFD methods are greater than 0.85, demonstrating good quantitative agreement. The stereoscopic PIV study can be utilized as test case templates for ongoing efforts in cardiovascular CFD solver development. Likewise, it is envisaged that the patient-specific data may provide a benchmark for further studying hemodynamics of actual AAA, ILT, and their convolution effects under physiological conditions for clinical applications.

Left atrial ligation alters intracardiac flow patterns and the biomechanical landscape in the chick embryo.

Background: Hypoplastic left heart syndrome (HLHS) is a major human congenital heart defect that results in single ventricle physiology and high mortality. Clinical data indicate that intracardiac blood flow patterns during cardiac morphogenesis are a significant etiology. We used the left atrial ligation (LAL) model in the chick embryo to test the hypothesis that LAL immediately alters intracardiac flow streams and the biomechanical environment, preceding morphologic and structural defects observed in HLHS. Results: Using fluorescent dye injections, we found that intracardiac flow patterns from the right common cardinal vein, right vitelline vein, and left vitelline vein were altered immediately following LAL. Furthermore, we quantified a significant ventral shift of the right common cardinal and right vitelline vein flow streams. We developed an in silico model of LAL, which revealed that wall shear stress was reduced at the left atrioventricular canal and left side of the common ventricle. Conclusions: Our results demonstrate that intracardiac flow patterns change immediately following LAL, supporting the role of hemodynamics in the progression of HLHS. Sites of reduced WSS revealed by computational modeling are commonly affected in HLHS, suggesting that changes in the biomechanical environment may lead to abnormal growth and remodeling of left heart structures. Developmental Dynamics 243:652–662, 2014.

Presurgical evaluation of Fontan connection options for patients with apicocaval juxtaposition using computational fluid dynamics.

Apicocaval juxtaposition (ACJ) is a rare congenital heart defect associated with single ventricle physiology where optimal positioning of the Fontan conduit for completion of total cavopulmonary connection (TCPC) is still controversial. In ACJ, the cardiac apex is ipsilateral with the inferior vena cava (IVC), risking kinking and collapse of the Fontan conduit at the apex of the heart. The purpose of this study is to evaluate two viable routes for Fontan conduit connection in patients with ACJ, using computational fluid dynamics. Internal energy loss evaluations were used to determine contribution of conduit curvature to the energy efficiency of each cavopulmonary anastomosis configuration. This percentage of energy loss contribution was found to be greater in the case of a curved extracardiac conduit connection (44%, 4.1 mW) traveling behind the ventricular apex, connecting the IVC to the left pulmonary artery, than the straighter lateral tunnel conduit (6%, 1.4 mW) installed through the ventricular apex. In contrast, net energy loss across the anastomosis was significantly lower with extracardiac TCPC (9.3 mW) in comparison with lateral tunnel TCPC (23.2 mW), highlighting that a curved Fontan conduit is favorable provided that it is traded off for a superior cavopulmonary connection efficiency. Therefore, a relatively longer and curved Fontan conduit has been demonstrated to be a suitable connection option independent of anatomical situations.

Rapid Quantification of Mean Myocardial Wall Velocity in Ischemic Cardiomyopathy by Cardiac Magnetic Resonance: An Index of Cardiac Functional Abnormalities during the Cardiac Cycle

Background: Characterization of regional Left Ventricular (LV) function in Heart Failure (HF) patients may have application in assessing response and choice of HF therapy. Methods: 4D LV endocardial function from cine cMR data in 12 NYHA Class III /IV Heart Failure (HF) patients with antero-apical aneurysms was compared against normal controls. The endocardial surface contours were extracted at each cardiac phase using Medviso Segment. An in-house plugin was used to compute a signed Hausdorff Distance (HD), establishing point-correspondences between endocardial segmentations at consecutive cardiac phases and tracking their surface motion throughout the cardiac cycle. The LV averaged myocardial velocity, was computed based on phase-to-phase displacement at several uniformly spaced endocardial surface points and the displacement histories of each point were recorded. The average and standard deviation in the endocardiumaveraged displacement history characteristic curves for the HF patients and normal controls were compared. Results: The endocardium-averaged phase-to-phase displacement history curves of HF cohort revealed diminished displacement magnitudes; and indistinct individual end-systolic and diastole instants with large standard deviations across the LV. When the cumulative phase to phase displacement plots for controls and HF patients were compared, there was a statistically significant difference between the two characteristic curves during the ejection phase, early filling phase of diastole, and the end diastolic phase. Conclusions: Shape-based Mean Myocardial Velocity (MMV) characterization constitutes a promising paradigm which may have application to prognosticate response to therapy.

Total cavopulmonary connection in patients with apicocaval juxtaposition: optimal conduit route using preoperative angiogram and flow simulation †

OBJECTIVES Single ventricle with apicocaval juxtaposition (ACJ) is a rare, complex anomaly, in which the optimal position of the conduit for completion of total cavopulmonary connection (TCPC) is still controversial. The purpose of this study was to identify a preoperative method for optimal conduit position using the IVC anatomy and computational fluid dynamics (CFD). METHODS Twenty-four patients with ACJ (5.3 ± 5.7 years) who underwent TCPC were enrolled. A conduit was placed ipsilateral to the cardiac apex in each of 11 patients, of which 9 were intra-atrial and 2 extracardiac (group A) and, in a further 13 patients, extracardiac on the contralateral side (group B). As control, 10 patients with tricuspid atresia were also enrolled (group C). The location of the IVC in relation to the spine was evaluated from the frontal view of preoperative angiogram, using the following index: IVC-index = IVC width overlapping the vertebra/width of the vertebra × 100%. Energy loss was calculated by CFD simulation. RESULTS IVC-index of group B was larger than groups A and C (45 ± 26 vs. 20 ± 21 and 28 ± 19%, P = 0.03). Postoperative catheterizations showed that, due to its curvature, conduit length in group B was significantly longer than the others (65 ± 12 vs. 36 ± 14 and 44 ± 10 mm, P < 0.001), although there was no statistical difference in central venous pressure or cardiac output. CFD studies revealed less energy loss in group A conduits compared with group B (1.6 ± 0.3 vs. 3.6 ± 0.6 mW, P = 0.05), although this did not appear to be clinically significant. Moreover, CFD simulation showed significant energy loss within the Fontan circulation when the conduit was either compressed or kinked: 4.9 and 18.2 mW respectively. CONCLUSIONS In patients with ACJ, placement of a straighter and shorter conduit on the ventricular apical side provides better laminar blood flow with less energy loss. However, conduit compression and kinking are far more detrimental to the Fontan circulation. A preoperative IVC-index is pivotal for avoiding these factors and deciding the optimal conduit route.

Hemodynamics of the hepatic venous three-vessel confluences using particle image velocimetry.

Despite rapid advancements in the patient-specific hemodynamic analysis of systemic arterial anatomies, limited attention has been given to the characterization of major venous flow components, such as the hepatic venous confluence. A detailed investigation of hepatic flow structures is essential to better understand the origin of characteristic abnormal venous flow patterns observed in patients with cardiovascular venous disease. The present study incorporates transparent rapid-prototype replicas of two pediatric hepatic venous confluence anatomies and two-component particle image velocimetry to investigate the primary flow structures influencing the inferior vena cava outflow. Novel jet flow regimes are reported at physiologically relevant mean venous conditions. The sensitivity of fluid unsteadiness and hydraulic resistance to multiple-inlet flow regimes is documented. Pressure drop measurements, jet flow characterization, and blood damage assessments are also performed. Results indicate that the orientation of the inlets significantly influences the major unsteady flow structures and power loss characteristics of this complex venous flow junction. Compared to out-of-plane arranged inlet vessel configuration, the internal flow field observed in planar inlet configurations was less sensitive to the venous inlet flow split. Under pathological flow conditions, the effective pressure drop increased as much as 77% compared to the healthy flow state. Experimental flow field results presented here can serve as a benchmark case for the surgical optimization of complex anatomical confluences including visceral hemodynamics as well as for the experimental validation of high-resolution computational fluid dynamics solvers applied to anatomical confluences with multiple inlets and outlets.

cMRI-BED: A novel informatics framework for cardiac MRI biomarker extraction and discovery applied to pediatric cardiomyopathy classification

BackgroundPediatric cardiomyopathies are a rare, yet heterogeneous group of pathologies of the myocardium that are routinely examined clinically using Cardiovascular Magnetic Resonance Imaging (cMRI). This gold standard powerful non-invasive tool yields high resolution temporal images that characterize myocardial tissue. The complexities associated with the annotation of images and extraction of markers, necessitate the development of efficient workflows to acquire, manage and transform this data into actionable knowledge for patient care to reduce mortality and morbidity.MethodsWe develop and test a novel informatics framework called cMRI-BED for biomarker extraction and discovery from such complex pediatric cMRI data that includes the use of a suite of tools for image processing, marker extraction and predictive modeling. We applied our workflow to obtain and analyze a dataset of 83 de-identified cases and controls containing cMRI-derived biomarkers for classifying positive versus negative findings of cardiomyopathy in children. Bayesian rule learning (BRL) methods were applied to derive understandable models in the form of propositional rules with posterior probabilities pertaining to their validity. Popular machine learning methods in the WEKA data mining toolkit were applied using default parameters to assess cross-validation performance of this dataset using accuracy and percentage area under ROC curve (AUC) measures.ResultsThe best 10-fold cross validation predictive performance obtained on this cMRI-derived biomarker dataset was 80.72% accuracy and 79.6% AUC by a BRL decision tree model, which is promising from this type of rare data. Moreover, we were able to verify that mycocardial delayed enhancement (MDE) status, which is known to be an important qualitative factor in the classification of cardiomyopathies, is picked up by our rule models as an important variable for prediction.ConclusionsPreliminary results show the feasibility of our framework for processing such data while also yielding actionable predictive classification rules that can augment knowledge conveyed in cardiac radiology outcome reports. Interactions between MDE status and other cMRI parameters that are depicted in our rules warrant further investigation and validation. Predictive rules learned from cMRI data to classify positive and negative findings of cardiomyopathy can enhance scientific understanding of the underlying interactions among imaging-derived parameters.