Turbulence in blood damage modeling.

PURPOSE: To account for the impact of turbulence in blood damage modeling, a novel approach based on the generation of instantaneous flow fields from RANS simulations is proposed. METHODS: Turbulent flow in a bileaflet mechanical heart valve was simulated using RANS-based (SST k-ω) flow solver using FLUENT 14.5. The calculated Reynolds shear stress (RSS) field is transformed into a set of divergence-free random vector fields representing turbulent velocity fluctuations using procedural noise functions. To consider the random path of the blood cells, instantaneous flow fields were computed for each time step by summation of RSS-based divergence-free random and mean velocity fields. Using those instantaneous flow fields, instantaneous pathlines and corresponding point-wise instantaneous shear stresses were calculated. For a comparison, averaged pathlines based on mean velocity field and respective viscous shear stresses together with RSS values were calculated. Finally, the blood damage index (hemolysis) was integrated along the averaged and instantaneous pathlines using a power law approach and then compared. RESULTS: Using RSS in blood damage modeling without a correction factor overestimates damaging stress and thus the blood damage (hemolysis). Blood damage histograms based on both presented approaches differ. CONCLUSIONS: A novel approach to calculate blood damage without using RSS as a damaging parameter is established. The results of our numerical experiment support the hypothesis that the use of RSS as a damaging parameter should be avoided.

The Computational Fluid Dynamics Rupture Challenge 2013--Phase II: Variability of Hemodynamic Simulations in Two Intracranial Aneurysms.

With the increased availability of computational resources, the past decade has seen a rise in the use of computational fluid dynamics (CFD) for medical applications. There has been an increase in the application of CFD to attempt to predict the rupture of intracranial aneurysms, however, while many hemodynamic parameters can be obtained from these computations, to date, no consistent methodology for the prediction of the rupture has been identified. One particular challenge to CFD is that many factors contribute to its accuracy; the mesh resolution and spatial/temporal discretization can alone contribute to a variation in accuracy. This failure to identify the importance of these factors and identify a methodology for the prediction of ruptures has limited the acceptance of CFD among physicians for rupture prediction. The International CFD Rupture Challenge 2013 seeks to comment on the sensitivity of these various CFD assumptions to predict the rupture by undertaking a comparison of the rupture and blood-flow predictions from a wide range of independent participants utilizing a range of CFD approaches. Twenty-six groups from 15 countries took part in the challenge. Participants were provided with surface models of two intracranial aneurysms and asked to carry out the corresponding hemodynamics simulations, free to choose their own mesh, solver, and temporal discretization. They were requested to submit velocity and pressure predictions along the centerline and on specified planes. The first phase of the challenge, described in a separate paper, was aimed at predicting which of the two aneurysms had previously ruptured and where the rupture site was located. The second phase, described in this paper, aims to assess the variability of the solutions and the sensitivity to the modeling assumptions. Participants were free to choose boundary conditions in the first phase, whereas they were prescribed in the second phase but all other CFD modeling parameters were not prescribed. In order to compare the computational results of one representative group with experimental results, steady-flow measurements using particle image velocimetry (PIV) were carried out in a silicone model of one of the provided aneurysms. Approximately 80% of the participating groups generated similar results. Both velocity and pressure computations were in good agreement with each other for cycle-averaged and peak-systolic predictions. Most apparent "outliers" (results that stand out of the collective) were observed to have underestimated velocity levels compared to the majority of solutions, but nevertheless identified comparable flow structures. In only two cases, the results deviate by over 35% from the mean solution of all the participants. Results of steady CFD simulations of the representative group and PIV experiments were in good agreement. The study demonstrated that while a range of numerical schemes, mesh resolution, and solvers was used, similar flow predictions were observed in the majority of cases. To further validate the computational results, it is suggested that time-dependent measurements should be conducted in the future. However, it is recognized that this study does not include the biological aspects of the aneurysm, which needs to be considered to be able to more precisely identify the specific rupture risk of an intracranial aneurysm.

Non-dimensional modeling in flow simulation studies of coronary arteries including side-branches: a novel diagnostic tool in coronary artery disease.

AIMS: Blood flow, vascular shape and size and local remodeling of the vascular wall are linked through wall shear stress (WSS) signaling. Inter-individual comparison of shape and WSS is hampered by large differences in size of flow and shape. We performed non-dimensional modeling to discriminate different types of coronary artery remodeling based on WSS patterns and vessel morphology. METHODS AND RESULTS: Blood flow was simulated in three-dimensional reconstructed right coronary artery trees from seven controls, five patients with coronary artery disease (CAD) and five patients with aneurysmatic CAD (AnCAD) classified by expert visual diagnosis. A discriminant model using low WSS area, a remodeling index, and cross-correlation of WSS in main trunks and complete trees (K) as non-dimensional parameters classified CAD and AnCAD correctly and identified three patients with high risk profile and functional disease in controls. The new model was compared with discriminant analysis of identical cases simulated without side-branches. The inclusion of K (information from side-branches) and replacement of the mean diameter by a non-dimensional remodeling index improved the model. We found significant (p<0.005) gender differences in the remodeling index. CONCLUSION: The combination of non-dimensional modeling and WSS profiling should be further investigated as a novel diagnostic tool in CAD beyond local stenosis.

Flow simulation studies in coronary arteries--impact of side-branches.

AIMS: Wall shear stress (WSS) may induce local remodeling of the vascular wall and the WSS pattern in turn depends on vascular geometry. We aimed to elucidate the impact of side-branches on local WSS. METHODS AND RESULTS: Steady numerical flow simulation studies were performed in three-dimensional reconstructed right coronary artery (RCA) trees. RCA from seven controls, five patients with coronary artery disease (CAD) and five patients with aneurysmatic CAD (AnCAD) classified by expert visual diagnosis were studied. Then three transient flow simulations were performed with cases representative for each group in order to evaluate the impact of pulsatile flow simulation. As vascular size and flow rates vary considerably between patients, non-dimensional approaches were applied for group comparison. A point-to-point comparison of the WSS in the same tree with and without side-branches revealed local differences in WSS of up to 12.0 Pa. This was caused by a reduction of volume flow of up to 78.7% in the trunk. Differences are not only limited to bifurcation sites but also affect local narrowings and strongly curved segments. The point-to-point comparison of steady and transient simulations found an average increase of WSS of below 7% in transient simulations. No significant differences were found between histograms of pulsatile and steady simulations, showing a high cross-correlation of >0.97. CONCLUSION: Side-branches must not be neglected in numerical flow simulation (steady and transient) studies. Steady simulations are valid for an assessment of time-averaged WSS distributions.