Comparison of CFD and MRI Flow and Velocities in an In Vitro Large Artery Bypass Graft Model

3D MR flow analysis in realistic rapid-prototyping model systems of the thoracic aorta: comparison with in vivo data and computational fluid dynamics in identical vessel geometries

The knowledge of local vascular anatomy and function in the human body is of high interest for the diagnosis and treatment of cardiovascular disease. A comprehensive analysis of the hemodynamics in the thoracic aorta is presented based on the integration of flow-sensitive 4D MRI with state-of-the-art rapid prototyping technology and computational fluid dynamics (CFD). Rapid prototyping was used to transform aortic geometries as measured by contrast-enhanced MR angiography into realistic vascular models with large anatomical coverage. Integration into a flow circuit with patient-specific pulsatile in-flow conditions and application of flow-sensitive 4D MRI permitted detailed analysis of local and global 3D flow dynamics in a realistic vascular geometry. Visualization of characteristic 3D flow patterns and quantitative comparisons of the in vitro experiments with in vivo data and CFD simulations in identical vascular geometries were performed to evaluate the accuracy of vascular model systems. The results indicate the potential of such patient-specific model systems for detailed experimental simulation of realistic vascular hemodynamics. Further studies are warranted to examine the influence of refined boundary conditions of the human circulatory system such as fluid-wall interaction and their effect on normal and pathological blood flow characteristics associated with vascular geometry.

Feasibility of patient specific aortic blood flow CFD simulation

Patient specific modelling of the blood flow through the human aorta is performed using computational fluid dynamics (CFD) and magnetic resonance imaging (MRI). Velocity patterns are compared between computer simulations and measurements. The workflow includes several steps: MRI measurement to obtain both geometry and velocity, an automatic levelset segmentation followed by meshing of the geometrical model and CFD setup to perform the simulations follwed by the actual simulations. The computational results agree well with the measured data.

Effects of Exercise and Respiration on Hemodynamic Efficiency in CFD Simulations of the Total Cavopulmonary Connection

Congenital heart defects with a single functional ventricle, such as hypoplastic left heart syndrome and tricuspid atresia, require a staged surgical approach to separate the systemic and pulmonary circulations. Ultimately, the venous or pulmonary side of the heart is bypassed by directly connecting the vena cava to the pulmonary arteries with a modified t-shaped junction. The Fontan procedure (total cavopulmonary connection, TCPC) completes this process of separation. To date, computational fluid dynamics (CFD) simulations in this low pressure, passive flow, intrathoracic system have neglected the presumed important effects of respiration on physiology and higher "stress" states such as with exercise have never been considered. We hypothesize that incorporating effects of respiration and exercise would provide more realistic estimates of TCPC performance. Time-dependent, 3D blood flow simulations are performed by a custom finite element solver for two patient-specific Fontan models with a novel respiration model, developed to generate physiologic time-varying flow conditions. Blood flow features, pressure, and energy efficiency are analyzed at rest and with increasing flow rates to simulate exercise conditions. The simulations produce realistic pressure and flow data, comparable to that measured by catheterization and echocardiography, and demonstrate substantial increases in energy dissipation (i.e. decreased performance) with exercise and respiration due to increasing intensity of small scale vortices in the flow. As would be expected, these changes are highly dependent on patient-specific anatomy and Fontan geometry. We propose that respiration and exercise should be incorporated into TCPC CFD simulations to provide increasingly realistic evaluations of TCPC performance.

Aortoiliac hemodynamic and morphologic adaptation to chronic spinal cord injury

BACKGROUND

Reduced lower limb blood flow and resistive hemodynamic conditions potentially promote aortic inflammation and aneurysmal degeneration. We used abdominal ultrasonography, magnetic resonance imaging, and computational flow modeling to determine the relationship between reduced infrarenal aortic blood flow in chronic spinal cord injury (SCI) subjects and risk for abdominal aortic aneurysm (AAA) disease.

METHODS

Aortic diameter in consecutive SCI subjects (n = 123) was determined via transabdominal ultrasonography. Aortic anatomic and physiologic data were acquired via magnetic resonance angiography (MRA; n = 5) and cine phase-contrast magnetic resonance flow imaging (n = 4) from SCI subjects whose aortic diameter was less than 3.0 cm by ultrasonography. Computational flow models were constructed from magnetic resonance data sets. Results were compared with those obtained from ambulatory control subjects (ultrasonography, n = 129; MRA/phase-contrast magnetic resonance flow imaging, n = 6) who were recruited at random from a larger pool of risk factor-matched individuals without known AAA disease.

RESULTS

Age, sex distribution, and smoking histories were comparable between the SCI and control groups. In the SCI group, time since injury averaged 26 +/- 13 years (mean +/- SD). Aortic diameter was larger (P < .01), and the prevalence of large (> or = 2.5 cm; P < .01) or aneurysmal (> or = 3.0 cm; P < .05) aortas was greater in SCI subjects. Paradoxically, common iliac artery diameters were reduced in SCI subjects (< 1.0 cm; 48% SCI vs 26% control; P < .0001). Focal preaneurysmal enlargement was noted in four of five SCI subjects by MRA. Flow modeling revealed normal flow volume, biphasic and reduced oscillatory flow, slower pressure decay, and reduced wall shear stress in the SCI infrarenal aorta.

CONCLUSIONS

Characteristic aortoiliac hemodynamic and morphologic adaptations occur in response to chronic SCI. Slower aortic pressure decay and reduced wall shear stress after SCI may contribute to mural degeneration, enlargement, and an increased prevalence of AAA disease.

Validation of CFD simulations of cerebral aneurysms with implication of geometric variations

BACKGROUND

Computational fluid dynamics (CFD) simulations using medical-image-based anatomical vascular geometry are now gaining clinical relevance. This study aimed at validating the CFD methodology for studying cerebral aneurysms by using particle image velocimetry (PIV) measurements, with a focus on the effects of small geometric variations in aneurysm models on the flow dynamics obtained with CFD.

METHOD OF APPROACH

An experimental phantom was fabricated out of silicone elastomer to best mimic a spherical aneurysm model. PIV measurements were obtained from the phantom and compared with the CFD results from an ideal spherical aneurysm model (S1). These measurements were also compared with CFD results, based on the geometry reconstructed from three-dimensional images of the experimental phantom. We further performed CFD analysis on two geometric variations, S2 and S3, of the phantom to investigate the effects of small geometric variations on the aneurysmal flow field. Results. We found poor agreement between the CFD results from the ideal spherical aneurysm model and the PIV measurements from the phantom, including inconsistent secondary flow patterns. The CFD results based on the actual phantom geometry, however, matched well with the PIV measurements. CFD of models S2 and S3 produced qualitatively similar flow fields to that of the phantom but quantitatively significant changes in key hemodynamic parameters such as vorticity, positive circulation, and wall shear stress.

CONCLUSION

CFD simulation results can closely match experimental measurements as long as both are performed on the same model geometry. Small geometric variations on the aneurysm model can significantly alter the flow-field and key hemodynamic parameters. Since medical images are subjected to geometric uncertainties, image-based patient-specific CFD results must be carefully scrutinized before providing clinical feedback.

AAA Disease: Mechanism, Stratification, and Treatment

Abdominal aortic aneurysm (AAA) is a common and frequently lethal disease of older Americans. No medical therapy has been proven effective in retarding progression of small AAAs prior to surgical repair. With the emerging ability of magnetic resonance (MR) flow imaging and MR-based computational analysis to define aortic hemodynamic conditions, and bio-imaging strategies to monitor aortic inflammation real time in vivo, the opportunity now exists to confirm the potential value of medical interventions such as supervised exercise training as first line therapy for small AAA disease.

Allometric scaling of wall shear stress from mice to humans: quantification using cine phase-contrast MRI and computational fluid dynamics

Allometric scaling laws relate structure or function between species of vastly different sizes. They have rarely been derived for hemodynamic parameters known to affect the cardiovascular system, e.g., wall shear stress (WSS). This work describes noninvasive methods to quantify and determine a scaling law for WSS. Geometry and blood flow velocities in the infrarenal aorta of mice and rats under isoflurane anesthesia were quantified using two-dimensional magnetic resonance angiography and phase-contrast magnetic resonance imaging at 4.7 tesla. Three-dimensional models constructed from anatomic data were discretized and used for computational fluid dynamic simulations using phase-contrast velocity imaging data as inlet boundary conditions. WSS was calculated along the infrarenal aorta and compared between species to formulate an allometric equation for WSS. Mean WSS along the infrarenal aorta was significantly greater in mice and rats compared with humans (87.6, 70.5, and 4.8 dyn/cm2, P < 0.01), and a scaling exponent of −0.38 ( R2 = 0.92) was determined. Manipulation of the murine genome has made small animal models standard surrogates for better understanding the healthy and diseased human cardiovascular system. It has therefore become increasingly important to understand how results scale from mouse to human. This noninvasive methodology provides the opportunity to serially quantify changes in WSS during disease progression and/or therapeutic intervention.

Predicting changes in blood flow in patient-specific operative plans for treating aortoiliac occlusive disease

Traditionally, a surgeon will select a procedure for a particular patient on the basis of past experience with patients with a similar state of disease. The experience gained from this patient will be selectively used when treating the next patient with similar symptoms. This article describes a surgical planning system that was developed to enable a vascular surgeon to create and test alternative operative plans prior to surgery for a given patient. One-dimensional and three-dimensional hemodynamic (i.e., blood flow) simulations were performed for rest and exercise for operative plans for two aorto-femoral bypass patients and compared with actual postoperative data. The information obtained from one-dimensional (volume flow distribution and pressure losses) and three-dimensional (flow, pressure, and wall shear stress) hemodynamic simulations may be clinically relevant to vascular surgeons planning interventions.

Helical flow as fluid dynamic signature for atherogenesis risk in aortocoronary bypass. A numeric study

The main purpose of the study was to verify if helical flow, widely observed in several vessels, might be a signature of the blood dynamics of vein graft anastomosis. We investigated the existence of a relationship between helical flow structures and vascular wall indexes of atherogenesis in aortocoronary bypass models with different geometric features. In particular, we checked for the existence of a relationship between the degree of helical motion and the magnitude of oscillating shear stress in conventional hand-sewn proximal anastomosis. The study is based on the numerical evaluation of four bypass geometries that are attached to a simplified computer representation of the ascending aorta with different angulations relative to aortic outflow. The finite volume technique was used to simulate realistic graft fluid dynamics, including aortic compliance and proper aortic and graft flow rates. A quantitative method was applied to evaluate the level of helicity in the flow field associated with the four bypass models under investigation. A linear inverse relationship (R = -0.97) was found between the oscillating shear index and the helical flow index for the models under investigation. The results obtained support the hypothesis that an arrangement of the flow field in helical patterns may elicit damping in wall shear stress temporal gradients at the proximal graft. Accordingly, helical flow might play a significant role in preventing plaque deposition or in tuning the mechanotransduction pathways of cells. Therefore, results confirm that helical flow constitutes an important flow signature in vessels, and its strength as a fluid dynamic index (for instance in combination with magnetic resonance imaging flow visualization techniques) for risk stratification, in the activation of both mechanical and biological pathways leading to fibrointimal hyperplasia.