Computational approach to estimating the effects of blood properties on changes in intra-stent flow

In-silico analysis of hemodynamic indicators in idealized stented coronary arteries for varying stent indentation

In this work, we investigate the effects of stent indentation on hemodynamic indicators in stented coronary arteries. Our aim is to assess in-silico risk factors for in-stent restenosis (ISR) and thrombosis after stent implantation. The proposed model is applied to an idealized artery with Xience V stent for four indentation percentages and three mesh refinements. We analyze the patterns of hemodynamic indicators arising from different stent indentations and propose an analysis of time-averaged WSS (TAWSS), topological shear variation index (TSVI), oscillatory shear index (OSI), and relative residence time (RRT). We observe that higher indentations display higher frequency of critically low TAWSS, high TSVI, and non-physiological OSI and RRT. Furthermore, an appropriate mesh refinement is needed for accurate representation of hemodynamics in the stent vicinity. The results suggest that disturbed hemodynamics could play a role in the correlation between high indentation and ISR.

Predicting postinterventional rupture of intracranial aneurysms using arteriography-derived radiomic features after pipeline embolization

Background and purpose

Postinterventional rupture of intracranial aneurysms (IAs) remains a severe complication after flow diverter treatment. However, potential hemodynamic mechanisms underlying independent predictors for postinterventional rupture of IAs remain unclear. In this study, we employed arteriography-derived radiomic features to predict this complication.

Methods

We included 64 patients who underwent pipeline flow diversion for intracranial aneurysms, distinguishing between 16 patients who experienced postinterventional rupture and 48 who did not. We performed propensity score matching based on clinical and morphological factors to match these patients with 48 patients with postinterventional unruptured IAs at a 1:3 ratio. Postinterventional digital subtraction angiography were used to create five arteriography-derived perfusion parameter maps and then radiomics features were obtained from each map. Informative features were selected through the least absolute shrinkage and selection operator method with five-fold cross-validation. Subsequently, radiomics scores were formulated to predict the occurrence of postinterventional IA ruptures. Prediction performance was evaluated with the training and test datasets using area under the curve (AUC) and confusion matrix-derived metrics.

Results

Overall, 1,459 radiomics features were obtained, and six were selected. The resulting radiomics scores had high efficacy in distinguishing the postinterventional rupture group. The AUC and Youden index were 0.912 (95% confidence interval: 0.767-1.000) and 0.847 for the training dataset, respectively, and 0.938 (95% confidence interval, 0.806-1.000) and 0.800 for the testing dataset, respectively.

Conclusion

Radiomics scores generated using arteriography-derived radiomic features effectively predicted postinterventional IA ruptures and may aid in differentiating IAs at high risk of postinterventional rupture.

An Overview of In Vitro Drug Release Methods for Drug-Eluting Stents

The drug release profile of drug-eluting stents (DESs) is affected by a number of factors, including the formulation, design, and physicochemical properties of the utilized material. DES has been around for twenty years and despite its widespread clinical use, and efficacy in lowering the rate of target lesion restenosis, it still requires additional development to reduce side effects and provide long-term clinical stability. Unfortunately, for analyzing these implants, there is still no globally accepted in vitro test method. This is owing to the stent's complexity as well as the dynamic arterial compartments of the blood and vascular wall. The former is the source of numerous biological, chemical, and physical mechanisms that are more commonly observed in tissue, lumen, and DES. As a result, universalizing bio-relevant apparatus, suitable for liberation testing of such complex implants is difficult. This article aims to provide a comprehensive review of the methods used for in vitro release testing of DESs. Aspects related to the correlation of the release profiles in the cases of in vitro and in vivo are also addressed.

An in silico trials platform for the evaluation of effect of the arterial anatomy configuration on stent implantation. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:4213-4217. [PMID: 34892153 DOI: 10.1109/embc46164.2021.9629950] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The introduction of Bioresorbable Vascular Scaffolds (BVS) has revolutionized the treatment of atherosclerosis. InSilc is an in silico clinical trial (ISCT) platform in a Cloud-based environment used for the design, development and evaluation of BVS. Advanced multi-disciplinary and multiscale models are integrated in the platform towards predicting the short/acute and medium/long term scaffold performance. In this study, InSilc platform is employed in a use case scenario and demonstrates how the whole in silico pipeline allows the interpretation of the effect of the arterial anatomy configuration on stent implantation.

The effects of stenting on hemorheological parameters: An in vitro investigation under various blood flow conditions

 Despite their wide clinical usage, stent functionality may be compromised by complications at the site of implantation, including early/late stent thrombosis and occlusion. Although several studies have described the effect of fluid-structure interaction on local haemodynamics, there is yet limited information on the effect of the stent presence on specific hemorheological parameters. The current work investigates the red blood cell (RBC) mechanical behavior and physiological changes as a result of flow through stented vessels. Blood samples from healthy volunteers were prepared as RBC suspensions in plasma and in phosphate buffer saline at 45% haematocrit. Self-expanding nitinol stents were inserted in clear perfluoroalkoxy alkane tubing which was connected to a syringe, and integrated in a syringe pump. The samples were tested at flow rates of 17.5, 35 and 70 ml/min, and control tests were performed in non-stented vessels. For each flow rate, the sample viscosity, RBC aggregation and deformability, and RBC lysis were estimated. The results indicate that the presence of a stent in a vessel has an influence on the hemorheological characteristics of blood. The viscosity of all samples increases slightly with the increase of the flow rate and exposure. RBC aggregation and elongation index (EI) decrease as the flow rate and exposure increases. RBC lysis for the extreme cases is evident. The results indicate that the stresses developed in the stent area for the extreme conditions could be sufficiently high to influence the integrity of the RBC membrane.

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.

Numerical simulation on the effects of drug-eluting stents with different bending angles on hemodynamics and drug distribution

Implantation of drug-eluting stents in curved blood vessels may cause changes in hemodynamics and drug distribution, and even provoke in-stent restenosis. Due to the complexity of building three-dimensional (3-D) curved stent model, few studies have gone through such numerical simulations. In this study, three virtual stent models with different bending angles (0°, 30° and 90°) were developed to numerically investigate the distribution of wall shear stress (WSS) and drug concentration. The results showed that (1) the low-WSS regions in the inner bend of the stent models increased with the bending angles; (2) the drug concentration differed between the inner and outer bends of the stents but irrespective to the changes of bending angle; (3) the pattern of drug concentration for the curved stents found similar to that of the straight stents, and the phenomenon, 'proximal part low and distal part high' in the drug concentration showed in both the straight and curved stents. The increase in bending angles from 30° to 90° had little effect on the WSS and drug concentration; however, the largest effect of the curved stents was the remarkable difference of drug concentration between the inner and outer bends of the stents-about 20 %. Hence, it is feasible that quick analysis focused on the straight stents instead of the curved stents.

The quantification of hemodynamic parameters downstream of a Gianturco Zenith stent wire using newtonian and non-newtonian analog fluids in a pulsatile flow environment

Although deployed in the vasculature to expand vessel diameter and improve blood flow, protruding stent struts can create complex flow environments associated with flow separation and oscillating shear gradients. Given the association between magnitude and direction of wall shear stress (WSS) and endothelial phenotype expression, accurate representation of stent-induced flow patterns is critical if we are to predict sites susceptible to intimal hyperplasia. Despite the number of stents approved for clinical use, quantification on the alteration of hemodynamic flow parameters associated with the Gianturco Z-stent is limited in the literature. In using experimental and computational models to quantify strut-induced flow, the majority of past work has assumed blood or representative analogs to behave as Newtonian fluids. However, recent studies have challenged the validity of this assumption. We present here the experimental quantification of flow through a Gianturco Z-stent wire in representative Newtonian and non-Newtonian blood analog environments using particle image velocimetry (PIV). Fluid analogs were circulated through a closed flow loop at physiologically appropriate flow rates whereupon PIV snapshots were acquired downstream of the wire housed in an acrylic tube with a diameter characteristic of the carotid artery. Hemodynamic parameters including WSS, oscillatory shear index (OSI), and Reynolds shear stresses (RSS) were measured. Our findings show that the introduction of the stent wire altered downstream hemodynamic parameters through a reduction in WSS and increases in OSI and RSS from nonstented flow. The Newtonian analog solution of glycerol and water underestimated WSS while increasing the spatial coverage of flow reversal and oscillatory shear compared to a non-Newtonian fluid of glycerol, water, and xanthan gum. Peak RSS were increased with the Newtonian fluid, although peak values were similar upon a doubling of flow rate. The introduction of the stent wire promoted the development of flow patterns that are susceptible to intimal hyperplasia using both Newtonian and non-Newtonian analogs, although the magnitude of sites affected downstream was appreciably related to the rheological behavior of the analog. While the assumption of linear viscous behavior is often appropriate in quantifying flow in the largest arteries of the vasculature, the results presented here suggest this assumption overestimates sites susceptible to hyperplasia and restenosis in flow characterized by low and oscillatory shear.

Clinical benefits of integrating cardiac and vascular models

INTRODUCTION

Recent advances in computational methods and medical imaging techniques have enabled non-invasive exploration of cardiovascular pathologies, from cardiac level to complex arterial networks. The potential of cardiac and vascular modeling in guiding and monitoring therapies could be further extended through the integration of the two systems.

AREAS COVERED

This review includes advances in methods for cardiac electromechanics and vascular flow simulations. The results of a literature search depicting the state of the art in cardiac and vascular modeling are reviewed. The paper goes on to address the benefits and challenges of combined cardiovascular modeling, highlighting the relevance of specific cardiovascular features and implementation. Various alternative approaches and insights on future directions are presented and analyzed with respect to their applicability to clinical practice.

EXPERT OPINION

The article has emerged from the exploration of currently available cardiac and vascular mathematical tools and their corresponding clinical application. The summarized analysis suggests that future efforts should be aimed at developing more accurate and patient-specific mathematical models integrating cardiac and vascular functions to enhance the knowledge of cardiovascular pathologies.

Accurate Prediction of Wall Shear Stress in a Stented Artery: Newtonian Versus Non-Newtonian Models

A significant amount of evidence linking wall shear stress to neointimal hyperplasia has been reported in the literature. As a result, numerical and experimental models have been created to study the influence of stent design on wall shear stress. Traditionally, blood has been assumed to behave as a Newtonian fluid, but recently that assumption has been challenged. The use of a linear model; however, can reduce computational cost, and allow the use of Newtonian fluids (e.g., glycerine and water) instead of a blood analog fluid in an experimental setup. Therefore, it is of interest whether a linear model can be used to accurately predict the wall shear stress caused by a non-Newtonian fluid such as blood within a stented arterial segment. The present work compares the resulting wall shear stress obtained using two linear and one nonlinear model under the same flow waveform. All numerical models are fully three-dimensional, transient, and incorporate a realistic stent geometry. It is shown that traditional linear models (based on blood’s lowest viscosity limit, 3.5 Pa s) underestimate the wall shear stress within a stented arterial segment, which can lead to an overestimation of the risk of restenosis. The second linear model, which uses a characteristic viscosity (based on an average strain rate, 4.7 Pa s), results in higher wall shear stress levels, but which are still substantially below those of the nonlinear model. It is therefore shown that nonlinear models result in more accurate predictions of wall shear stress within a stented arterial segment.

Estimation of distributed arterial mechanical properties using a wave propagation model in a reverse way

To estimate arterial stiffness, different methods based either on distensibility, pulse wave velocity or a pressure-velocity loop, have been proposed. These methods can be employed to determine the arterial mechanical properties either locally or globally, e.g. averaged over an entire arterial segment. The aim of this study was to investigate the feasibility of a new method that estimates distributed arterial mechanical properties non-invasively. This new method is based on a wave propagation model and several independent ultrasound and pressure measurements. Model parameters (including arterial mechanical properties) are obtained from a reverse method in which differences between modeling results and measurements are minimized using a fitting procedure based on local sensitivity indices. This study evaluates the differences between in vivo measured and simulated blood pressure and volume flow waveforms at the brachial, radial and ulnar arteries of 6 volunteers. The estimated arterial Young's modulus range from 1.0 to 6.0MPa with an average of (3.8±1.7)MPa at the brachial artery and from 1.2 to 7.8MPa with an average of (4.8±2.2)MPa at the radial artery. A good match between measured and simulated waveforms and the realistic stiffness parameters indicate a good in vivo suitability.

Assessment of blood volume flow in slightly curved arteries from a single velocity profile

Non-invasive estimation of arterial blood volume flow (BVF) has become a central issue in assessment of cardiovascular risk. Poiseuille and Womersley approaches are commonly used to assess the BVF from centerline velocity, but both methods neglect the influence of curvature. Based on the assumption that the velocity in curved tubes as function of the circumferential position for a given radial position can be approximated by a cosine, the BVF can also be estimated by averaging velocities at opposite radial positions, referred to as the cosine theta model (CTM). This study investigates the accuracy of BVF estimation in slightly curved arteries for BVF waveforms obtained in the brachial artery of 6 volunteers. Computational fluid dynamics simulations were used to compute the influence of curvature on velocity profiles. The BVF was then estimated from the simulation results with the CTM and methods based on Poiseuille, Womersley and using the center stream velocity and the velocity waveform at the position where the maximum velocity is observed, and compared to the prescribed BVF. The simulations show that the influence of curvature is strongest when the flow decelerates. For Poiseuille and Womersley, the time average BVF was underestimated by maximally 10.4% and 7.8% for a radius of curvature of 50 and 100 mm, respectively. The estimation error is lower for the CTM and equals 4.2% and 1.2% for a radius of curvature of 50 and 100mm, respectively. From this study, we can conclude that the velocity waveform at the position of the maximum rather than the center stream velocity waveform combined with the Womersley method should be chosen. The CTM improves current estimation techniques if in-vivo velocity distributions are available.

Evaluation of the effect of stent strut profile on shear stress distribution using statistical moments

BACKGROUND

In-stent restenosis rates have been closely linked to the wall shear stress distribution within a stented arterial segment, which in turn is a function of stent design. Unfortunately, evaluation of hemodynamic performance can only be evaluated with long term clinical trials. In this work we introduce a set of metrics, based on statistical moments, that can be used to evaluate the hemodynamic performance of a stent in a standardized way. They are presented in the context of a 2D flow study, which analyzes the impact of different strut profiles on the wall shear stress distribution for stented coronary arteries.

RESULTS

It was shown that the proposed metrics have the ability to evaluate hemodynamic performance quantitatively and compare it to a common standard. In the context of the simulations presented here, they show that stent's strut profile significantly affect the shear stress distribution along the arterial wall. They also demonstrates that more streamlined profiles exhibit better hemodynamic performance than the standard square and circular profiles. The proposed metrics can be used to compare results from different research groups, and provide an improved method of quantifying hemodynamic performance in comparison to traditional techniques.

CONCLUSION

The strut shape found in the latest generations of stents are commonly dictated by manufacturing limitations. This research shows, however, that strut design can play a fundamental role in the improvement of the hemodynamic performance of stents. Present results show that up to 96% of the area between struts is exposed to wall shear stress levels above the critical value for the onset of restenosis when a tear-drop strut profile is used, while the analogous value for a square profile is 19.4%. The conclusions drawn from the non-dimensional metrics introduced in this work show good agreement with an ordinary analysis of the wall shear stress distribution based on the overall area exposed to critically low wall shear stress levels. The proposed metrics are able to predict, as expected, that more streamlined profiles perform better hemodynamically. These metrics integrate the entire morphology of the shear stress distribution and as a result are more robust than the traditional approach, which only compares the relative value of the local wall shear stress with a critical value of 0.5 Pa. In the future, these metrics could be employed to compare, in a standardized way, the hemodynamic performance of different stent designs.

Wall shear stress in intracranial self-expanding stents studied using ultra-high-resolution 3D reconstructions

BACKGROUND AND PURPOSE

Imaging of intracranial stents is constrained by resolution limits of current clinical imaging techniques providing insufficient visualization of deployment details and impeding its use for computational hemodynamic (CHD) simulations. The purpose of our study was to evaluate whether ultra-high-resolution MicroCT scans can illuminate detailed aspects of realistic in vitro stent deployment and serve as a reliable basis for CHD simulations of blood flow through self-expanding intracranial stents.

MATERIALS AND METHODS

A Neuroform Treo (NF) stent and an Enterprise (ENT) stent were deployed in identical straight polytetrafluoroethylene tubes filled with contrast agent. MicroCT scans were obtained at a spatial resolution of 14 mum and used for ultra-high-resolution 3D reconstructions. CHD simulations were performed, with particular emphasis on local flow behavior near the wall and struts. Flow differences between the geometrically different stents were studied.

RESULTS

MicroCT data revealed strut prolapse near the markers for the closed-cell design (ENT) stent and at some of the unconnected vertices of the open-cell design (NF) stent, which also showed some misalignments. CHD simulations showed that reverse wall shear stress occurred near some of the strut vertices and markers for the NF but only near the markers for the ENT.

CONCLUSIONS

This work demonstrates the feasibility of ultra-high-resolution MicroCT imaging in elucidating important details of intracranial stent deployment as a basis for accurate CHD simulations and in enabling a structural and hemodynamic study of realistically deployed stents with different geometry and design.

Luminal flow patterns dictate arterial drug deposition in stent-based delivery

Endovascular stents reside in a dynamic flow environment and yet the impact of flow on arterial drug deposition after stent-based delivery is only now emerging. We employed computational fluid dynamic modeling tools to investigate the influence of luminal flow patterns on arterial drug deposition and distribution. Flow imposes recirculation zones distal and proximal to the stent strut that extend the coverage of tissue absorption of eluted drug and induce asymmetry in tissue drug distribution. Our analysis now explains how the disparity in sizes of the two recirculation zones and the asymmetry in drug distribution are determined by a complex interplay of local flow and strut geometry. When temporal periodicity was introduced as a model of pulsatile flow, the net luminal flow served as an index of flow-mediated spatio-temporal tissue drug uptake. Dynamically changing luminal flow patterns are intrinsic to the coronary arterial tree. Coronary drug-eluting stents should be appropriately considered where luminal flow, strut design and pulsatility have direct effects on tissue drug uptake after local delivery.

Studies on Design Optimization of Coronary Stents

The stent design itself seems to be one of the factors responsible for restenosis. As a remedy, the present work attempts to perform a design optimization of coronary stents from a hemodynamic point of view. For the purpose, we have applied the principles of modern exploration of design space restricting ourselves to two-dimensional considerations. Width, thickness, and spacing of the struts of the stent formed the design variables. The objectives chosen for optimization were the vorticity generated, length of recirculation zone, and the reattachment distance in between the struts. Both semicircular and rectangular cross sections of stents were included. Starting with the range of design variables, sample stent cases were generated using Latin hypercube sampling. Objective functions were calculated for each of these by computing the two-dimensional flow using software FLUENT under the assumption of a steady, Newtonian flow considering a model stent with three struts. This was followed by Kriging to construct a response surface, which gives the relationship between the objectives and the design variables. The procedure gave nondominated fronts, which consist of optimized designs. Stents with minimum vorticity, with minimum recirculation distance, and the ones with maximum reattachment length in between struts were generated. The procedure is capable of producing the optimum set of design variables to achieve the prescribed objectives.