Coronary fractional flow reserve measurements of a stenosed side branch: a computational study investigating the influence of the bifurcation angle

Background Coronary hemodynamics and physiology specific for bifurcation lesions was not well understood. To investigate the influence of the bifurcation angle on the intracoronary hemodynamics of side branch (SB) lesions computational fluid dynamics simulations were performed. Methods A parametric model representing a left anterior descending—first diagonal coronary bifurcation lesion was created according to the literature. Diameters obeyed fractal branching laws. Proximal and distal main branch (DMB) stenoses were both set at 60 %. We varied the distal bifurcation angles (40°, 55°, and 70°), the flow splits to the DMB and SB (55 %:45 %, 65 %:35 %, and 75 %:25 %), and the SB stenoses (40, 60, and 80 %), resulting in 27 simulations. Fractional flow reserve, defined as the ratio between the mean distal stenosis and mean aortic pressure during maximal hyperemia, was calculated for the DMB and SB (FFRSB) for all simulations. Results The largest differences in FFRSB comparing the largest and smallest bifurcation angles were 0.02 (in cases with 40 % SB stenosis, irrespective of the assumed flow split) and 0.05 (in cases with 60 % SB stenosis, flow split 55 %:45 %). When the SB stenosis was 80 %, the difference in FFRSB between the largest and smallest bifurcation angle was 0.33 (flow split 55 %:45 %). By describing the ΔPSB−QSB relationship using a quadratic curve for cases with 80 % SB stenosis, we found that the curve was steeper (i.e. higher flow resistance) when bifurcation angle increases (ΔP = 0.451*Q + 0.010*Q 2 and ΔP = 0.687*Q + 0.017*Q 2 for 40° and 70° bifurcation angle, respectively). Our analyses revealed complex hemodynamics in all cases with evident counter-rotating helical flow structures. Larger bifurcation angles resulted in more pronounced helical flow structures (i.e. higher helicity intensity), when 60 or 80 % SB stenoses were present. A good correlation (R2 = 0.80) between the SB pressure drop and helicity intensity was also found. Conclusions Our analyses showed that, in bifurcation lesions with 60 % MB stenosis and 80 % SB stenosis, SB pressure drop is higher for larger bifurcation angles suggesting higher flow resistance (i.e. curves describing the ΔPSB−QSB relationship being steeper). When the SB stenosis is mild (40 %) or moderate (60 %), SB resistance is minimally influenced by the bifurcation angle, with differences not being clinically meaningful. Our findings also highlighted the complex interplay between anatomy, pressure drops, and blood flow helicity in bifurcations

Assessing the Near-Wall Hemodynamics in the Left Coronary Artery Using CFD

The objective of this thesis is to computationally investigate the flow mechanics and the near-wall hemodynamics associated with the different take-off angles in the left coronary artery of the human heart. From this study, we will be able to evaluate if the increase in the take-off angles of the left coronary artery will significantly increases or decrease the likelihood of plaque (atherosclerosis) buildup in the left coronary artery bifurcations. This study quantifies the effects of the varying take-off angles on the branches along the left anterior descending (LAD) of the left coronary artery using computational fluid dynamics (CFD) simulations. The study aims to compare five test cases of the different take off-angles of the left coronary artery (LCA) and four different branch angles between the LAD and the left circumflex (LCx). It also considered the branch angles of the coronary artery downstream the LAD. The idealized geometries used for this study were constructed in SolidWorks 2015 and imported as surface meshes into Star-CCM+, a commercially available CFD solver. In this study, the LCA inlet boundary conditions was set as a pulsatile mass flow inlet and flow split ratios were set for the outlets boundary conditions that are representations of a middle age man at rest. The nature of blood pulsatile flow characteristic was accounted for and the properties of blood which include the density (1,050Kg/m3) and dynamic viscosity (0.0046Pa) were obtained from previous research. The results from the simulations are compared using established scales for the parameters evaluated. The parameters evaluated were: (i) Oscillatory Shear Index (OSI); which quantifies the extent in which the blood flow changes direction as it flows (ii) Time Average Wall Shear Stress (TAWSS); which quantifies the average shear stress experienced by the wall of the artery and (ii) Relative Residence Time (RRT); which defined how long blood spends in a location during blood flow. These parameters are used to predict the likelihood of blood clots, atherosclerosis, endothelial damage, plaque formation, and aneurysm in the blood vessels. The data from the simulations were analyzed using functional macros to quantify and generate threshold values for the parameters. Computational Fluid Dynamics has gain more recognition in field of medicine because it has been used to obtain the various mechanic behaviors of most artificial implanted devices used for endovascular and cardiovascular treatments before these devices are used in patients’ treatment. This can be a useful insight in coronary stenting, solid and stress analysis of biodegradable stent and can also provide insight into stenting for more complex arterial networks like brain stent grafts. In addition, it is important to understand the hemodynamics of the LCA before carrying out stent graft or angioplasty procedures. This will help determine the effectiveness of the stent graft in the coronary artery

A Rapid and Computationally Inexpensive Method to Virtually Implant Current and Next-Generation Stents into Subject-Specific Computational Fluid Dynamics Models

Computational modeling is often used to quantify hemodynamic alterations induced by stenting, but frequently uses simplified device or vascular representations. Based on a series of Boolean operations, we developed an efficient and robust method for assessing the influence of current and next-generation stents on local hemodynamics and vascular biomechanics quantified by computational fluid dynamics. Stent designs were parameterized to allow easy control over design features including the number, width and circumferential or longitudinal spacing of struts, as well as the implantation diameter and overall length. The approach allowed stents to be automatically regenerated for rapid analysis of the contribution of design features to resulting hemodynamic alterations. The applicability of the method was demonstrated with patient-specific models of a stented coronary artery bifurcation and basilar trunk aneurysm constructed from medical imaging data. In the coronary bifurcation, we analyzed the hemodynamic difference between closed-cell and open-cell stent geometries. We investigated the impact of decreased strut size in stents with a constant porosity for increasing flow stasis within the stented basilar aneurysm model. These examples demonstrate the current method can be used to investigate differences in stent performance in complex vascular beds for a variety of stenting procedures and clinical scenarios

Engineering design of artificial vascular junctions for 3D printing

Vascular vessels, including arteries, veins and capillaries, are being printed using additive manufacturing technologies, also known as 3D printing. This paper demonstrates that it is important to follow the vascular design by nature as close as possible when 3D printing artificial vascular branches. In previous work, the authors developed an algorithm of computational geometry for constructing smooth junctions for 3D printing. In this work, computational fluid dynamics (CFDs) is used to compare the wall shear stress and blood velocity field for the junctions of different designs. The CFD model can reproduce the expected wall shear stress at locations remote from the junction. For large vessels such as veins, it is shown that ensuring the smoothness of the junction and using smaller joining angles as observed in nature is very important to avoid high wall shear stress and recirculation. The issue is however less significant for capillaries. Large joining angles make no difference to the hemodynamic behavior, which is also consistent with the fact that most capillary junctions have large joining angles. The combination of the CFD analysis and the junction construction method form a complete design method for artificial vascular vessels that can be 3D printed using additive manufacturing technologies

Coronary fractional flow reserve measurements of a stenosed side branch: A computational study investigating the influence of the bifurcation angle

Background: Coronary hemodynamics and physiology specific for bifurcation lesions was not well understood. To investigate the influence of the bifurcation angle on the intracoronary hemodynamics of side branch (SB) lesions computational fluid dynamics simulations were performed. Methods: A parametric model representing a left anterior descending-first diagonal coronary bifurcation lesion was created according to the literature. Diameters obeyed fractal branching laws. Proximal and distal main branch (DMB) stenoses were both set at 60%. We varied the distal bifurcation angles (40°, 55°, and 70°), the flow splits to the DMB and SB (55%:45%, 65%:35%, and 75%:25%), and the SB stenoses (40, 60, and 80%), resulting in 27 simulations. Fractional flow reserve, defined as the ratio between the mean distal stenosis and mean aortic pressure during maximal hyperemia, was calculated for the DMB and SB (FFRSB) for all simulations. Results: The largest differences in FFRSB comparing the largest and smallest bifurcation angles were 0.02 (in cases with 40% SB stenosis, irrespective of the assumed flow split) and 0.05 (in cases with 60% SB stenosis, flow split 55%:45%). When the SB stenosis was 80%, the difference in FFRSB between the largest and smallest bifurcation angle was 0.33 (flow split 55%:45%). By describing the PSB-QSB relationship using a quadratic curve for cases with 80% SB stenosis, we found that the curve was steeper (i.e. higher flow resistance) when bifurcation angle increases (P=0.451*Q+0.010*Q 2 and P=0.687*Q+0.017*Q 2 for 40° and 70° bifurcation angle, respectively). Our analyses revealed complex hemodynamics in all cases with evident counter-rotating helical flow structures. Larger bifurcation angles resulted in more pronounced helical flow structures (i.e. higher helicity intensity), when 60 or 80% SB stenoses were present. A good correlation (R2=0.80) between the SB pressure drop and helicity intensity was also found. Conclusions: Our analyses showed that, in bifurcation lesions with 60% MB stenosis and 80% SB stenosis, SB pressure drop is higher for larger bifurcation angles suggesting higher flow resistance (i.e. curves describing the PSB-QSB relationship being steeper). When the SB stenosis is mild (40%) or moderate (60%), SB resistance is minimally influenced by the bifurcation angle, with differences not being clinically meaningful. Our findings also highlighted the complex interplay between anatomy, pressure drops, and blood flow helicity in bifurcations

Blood flow dynamics at the pulmonary artery bifurcation

Knowledge of physiologic hemodynamics is a fundamental requirement to establish pathological findings. However, little is known about the normal flow fields in the pulmonary arteries, especially for children. The purpose of this study is to characterize flow patterns in the pulmonary artery bifurcation of healthy pediatric subjects using direct numerical simulations. A realistic geometry is obtained via statistical shape modeling, by averaging five subject-specific digital models extracted from cardiovascular magnetic resonance datasets of healthy volunteers. Boundary conditions are assigned to mimic physiological conditions at rest, corresponding to a peak Reynolds number equal to 3400 and a Womersley number equal to 15. Results show that the normal bifurcation is highly hemodynamically efficient, as measured by an energy dissipation index. The curvature of the pulmonary arteries is sufficiently small to prevent flow separation along the inner walls, and no signs of a turbulent-like state are found. In line with previous imaging studies, a helical structure protruding into the right pulmonary artery is detected, and its formation mechanism is elucidated in the paper. These findings might help to identify abnormal flow features in patients with altered anatomic and physiologic states, particularly those with repaired congenital heart disease.Peer ReviewedPostprint (published version

Numerical investigation of hemodynamic performance of a stent in the main branch of a coronary artery bifurcation

Introduction: The effect of a bare-metal stent on the hemodynamics in the main branch of a coronary artery bifurcation with a particular type of stenosis was numerically investigated by the computational fluid dynamics (CFD). Methods: Three-dimensional idealized geometry of bifurcation was constructed in Catia modelling commercial software package. The Newtonian blood flow was assumed to be incompressible and laminar. CFD was utilized to calculate the shear stress and blood pressure distributions on the wall of main branch. In order to do the numerical simulations, a commercial software package named as COMSOL Multiphysics 5.3 was employed. Two types of stent, namely, one-part stent and two-part stent were applied to prevent the build-up and progression of the atherosclerotic plaques in the main branch. Results: A particular type of stenosis in the main branch was considered in this research. It occurred before and after the side branch. Moreover, it was found that the main branch with an inserted one-part stent had the smallest region with the wall shear stress (WSS) below 0.5 Pa which was the minimum WSS in the main branch without the stenosis. Conclusion: The use of a one-part stent in the main branch of a coronary artery bifurcation for the aforementioned type of stenosis is recommended

Improving Cardiovascular Stent Design Using Patient-Specific Models and Shape Optimization

Stent geometry influences local hemodynamic alterations (i.e. the forces moving blood through the cardiovascular system) associated with adverse clinical outcomes. Computational fluid dynamics (CFD) is frequently used to quantify stent-induced hemodynamic disturbances, but previous CFD studies have relied on simplified device or vascular representations. Additionally, efforts to minimize stent-induced hemodynamic disturbances using CFD models often only compare a small number of possible stent geometries. This thesis describes methods for modeling commercial stents in patient-specific vessels along with computational techniques for determining optimal stent geometries that address the limitations of previous studies. An efficient and robust method was developed for virtually implanting stent models into patient-specific vascular geometries derived from medical imaging data. Models of commercial stent designs were parameterized to allow easy control over design features. Stent models were then virtually implanted into vessel geometries using a series of Boolean operations. This approach allowed stented vessel models to be automatically regenerated for rapid analysis of the contribution of design features to resulting hemodynamic alterations. The applicability of the method was demonstrated with patient-specific models of a stented coronary artery bifurcation and basilar trunk aneurysm to reveal how it can be used to investigate differences in hemodynamic performance in complex vascular beds for a variety of clinical scenarios. To identify hemodynamically optimal stents designs, a computational framework was constructed to couple CFD with a derivative-free optimization algorithm. The optimization algorithm was fully-automated such that solid model construction, mesh generation, CFD simulation and time-averaged wall shear stress (TAWSS) quantification did not require user intervention. The method was applied to determine the optimal number of circumferentially repeating stent cells (NC) for a slotted-tube stents and various commercial stents. Optimal stent designs were defined as those minimizing the area of low TAWSS. It was determined the optimal value of NC is dependent on the intrastrut angle with respect to the primary flow direction. Additionally, the geometries of current commercial stents were found to generally incorporate a greater NC than is hemodynamically optimal. The application of the virtual stent implantation and optimization methods may lead to stents with superior hemodynamic performance and the potential for improved clinical outcomes. Future in vivo studies are needed to validate the findings of the computational results obtained from the methods developed in this thesis