Wall shear stress variations in basilar tip aneurysms investigated with computational fluid dynamics

Integration of the computational fluid dynamics technique with MRI in aortic dissections

Short-term and long-term prognosis and their determining factors of Type III/Stanford B aortic dissections (TB-AD), which separate the aorta distal at the origin of the subclavian artery into a true lumen and false lumen, have been elusive: One quarter of patients thought to be treated successfully, either by medical or by surgical means, do not survive 3 years. Unfavorable hemodynamic conditions are believed to lead to false lumen pressure increases and complications. A better characterization of TB-AD hemodynamics may therefore impact therapeutic decision making and improve outcome. The large variations in TB-AD morphology and hemodynamics favor a patient-specific approach. Magnetic resonance imaging with its capability to provide high-resolution structural images of the lumen and aortic wall and also to quantify aortic flow and kinetics of an exogenous tracer is a promising clinical modality for developing a deeper understanding of TB-AD hemodynamics in an individual patient. With the information obtained with magnetic resonance imaging, computational fluid dynamics simulations can be performed to augment the image information. Here, an overview of the interplay of magnetic resonance imaging and computational fluid dynamics techniques is given illustrating the synergy of these two approaches toward a comprehensive morphological and hemodynamic characterization of TB-AD.

Temporal variations of wall shear stress parameters in intracranial aneurysms--importance of patient-specific inflow waveforms for CFD calculations

PURPOSE

To assess reliability of wall shear stress (WSS)calculations using computational fluid dynamics (CFD) dependent on inflow in internal carotid artery aneurysms (ICA).

MATERIALS AND METHODS

Six unruptured ICA aneurysms were studied. 3D computational meshes were created from 3D digital subtraction angiographic images (Axiom Artis dBA, Siemens Medical Solutions). Transient CFD simulations(Fluent, ANSYS Inc.) were performed for two inflow conditions: (1) idealized averaged waveform from normal subjects (ID) and (2) patient-specific waveform (PS)measured with 2D phase contrast magnetic resonance imaging. Stability of calculation was assessed by comparing mean WSS (<WSS>), temporal wall shear stress magnitude variation (Delta WSS), and oscillatory shear index(OSI, a measure of variation in the WSS direction) on the aneurysmal wall for both conditions.

RESULTS

For all cases, mean relative difference (PS-ID) of WSS (<WSS>) was -15% (range -32% to 11%). Mean Delta WSS difference was -29.3% ( -100% to 67%). Mean OSI difference was 7.5% (-12% to 40%). Large variations in histograms of these parameters were noted.

CONCLUSION

For accurate calculations of WSS parameters,patient-specific information on physiological flow may be necessary. Results obtained with averaged or idealized flow waveforms may have to be interpreted with caution.

Tetrahedral and polyhedral mesh evaluation for cerebral hemodynamic simulation--a comparison

Computational fluid dynamic (CFD) based on patient-specific medical imaging data has found widespread use for visualizing and quantifying hemodynamics in cerebrovascular disease such as cerebral aneurysms or stenotic vessels. This paper focuses on optimizing mesh parameters for CFD simulation of cerebral aneurysms. Valid blood flow simulations strongly depend on the mesh quality. Meshes with a coarse spatial resolution may lead to an inaccurate flow pattern. Meshes with a large number of elements will result in unnecessarily high computation time which is undesirable should CFD be used for planning in the interventional setting. Most CFD simulations reported for these vascular pathologies have used tetrahedral meshes. We illustrate the use of polyhedral volume elements in comparison to tetrahedral meshing on two different geometries, a sidewall aneurysm of the internal carotid artery and a basilar bifurcation aneurysm. The spatial mesh resolution ranges between 5,119 and 228,118 volume elements. The evaluation of the different meshes was based on the wall shear stress previously identified as a one possible parameter for assessing aneurysm growth. Polyhedral meshes showed better accuracy, lower memory demand, shorter computational speed and faster convergence behavior (on average 369 iterations less).