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

Haemodynamics analysis of carotid artery stenosis and carotid artery stenting

Carotid stenosis is a local narrowing of the carotid artery, and is usually found in the internal carotid artery. The presence of a high-degree stenosis in a carotid artery may provoke transition from laminar to turbulent flow during part of the cardiac cycle. Turbulence in blood flow can influence haemodynamic parameters such as velocity profiles, shear stress and pressure, which are important in wall remodelling. Patients with severe stenosis could be treated with a minimally invasive clinical procedure, carotid artery stenting (CAS). Although CAS has been widely adopted in clinical practice, the complication of in-stent restenosis (ISR) has been reported after CAS. The incidence of ISR is influenced by stent characteristics and vessel geometry, and correlates strongly with regions of neointimal hyperplasia (NH). Therefore, the main purpose of this study is to provide more insights into the haemodynamics in stenosed carotid artery and in post-CAS geometries via computational simulation. The first part of the thesis presents a computational study on flow features in a stenotic carotid artery bifurcation using two computational approaches, large eddy simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) incorporating the Shear Stress Transport model with the γ-Reθ transition (SST-Tran) models. The computed flow patterns are compared with those measured with particle image velocimetry (PIV). The results show that both SST-Tran and LES can predict the PIV results reasonably well, but LES is more accurate especially at locations distal to the stenosis where flow is highly disturbed. The second part of the thesis is to determine how stent strut design may influence the development of ISR at the carotid artery bifurcation following CAS. Key parameters that can be indicative of ISR are obtained for different stent designs and compared; these include low and oscillating wall shear stress (WSS), high residence time, and wall stress. A computationally efficient methodology is employed to reproduce stent strut geometry. This method facilitates the accurate reconstruction of actual stent geometry and details of strut configuration and its inclusion in the fluid domain. Computational simulations for flow patterns and low-density lipoprotein (LDL) transport are carried out in order to investigate spatial and temporal variations of WSS and LDL accumulation in the stented carotid geometries. Furthermore, finite element (FE) analysis is performed to evaluate the wall stress distribution with different stent designs. The results reveal that the closed-cell stent design is more likely to create atheroprone and procoagulant flow conditions, causing larger area to be exposed to low wall shear stress (WSS), elevated oscillatory shear index, as well as to induce higher wall stress compared to the open-cell stent design. This study also demonstrates the suitability of SST-Tran and LES models in capturing the presence of complex flow patterns in post-stenotic region.Open Acces

Computed poststenotic flow instabilities correlate phenotypically with vibrations measured using laser Doppler vibrometry : perspectives for a promising in vivo device for early detection of moderate and severe carotid stenosis

Early detection of asymptomatic carotid stenosis is crucial for treatment planning in the prevention of ischemic stroke. Auscultation, the current first-line screening methodology, comes with severe limitations that create urge for novel and robust techniques. Laser Doppler vibrometer (LDV) is a promising tool for inferring carotid stenosis by measuring stenosis-induced vibrations. The goal of the current study was to evaluate the feasibility of LDV for carotid stenosis detection. LDV measurements on a carotid phantom were used to validate our previously verified high-resolution computational fluid dynamics methodology, which was used to evaluate the impact of flowrate, flow split, and stenosis severity on the poststenotic intensity of flow instabilities (IFI). We evaluated sensitivity, specificity, and accuracy of using IFI for stenoses detection. Linear regression analyses showed that computationally derived pressure fluctuations correlated (R2 = 0.98) with LDV measurements of stenosis-induced vibrations. The flowrate of stenosed vessels correlated (R2 = 0.90) with the presence of poststenotic instabilities. Receiver operating characteristic analyses of power spectra revealed that the most relevant frequency bands for the detection of moderate (56–76%) and severe (86–96%) stenoses were 80–200 Hz and 0–40 Hz, respectively. Moderate stenosis was identified with sensitivity and specificity of 90%; values decreased to 70% for severe stenosis. The use of LDV as screening tool for asymptomatic stenosis can potentially provide improved accuracy of current screening methodologies for early detection. The applicability of this promising device for mass screening is currently being evaluated clinically

Predictive haemodynamics in a one-dimensional human carotid artery bifurcation. Part II: application to graft design

A Bayesian surrogate modelling technique is proposed that may be able to predict an optimal bypass graft configuration for patients suffering with stenosis in the internal carotid artery (ICA). At the outset, this statistical technique is considered as a means for identifying key geometric parameters influencing haemodynamics in the human carotid bifurcation. This methodology uses a design of experiments (DoE) technique to generate candidate geometries for flow analysis. A pulsatile one dimensional Navier-Stokes solver incorporating fluid-wall interactions for a Newtonian fluid which predicts pressure and flow in the carotid bifurcation (comprising a stenosed segment in the internal carotid artery) is used for the numerical simulations. Two metrics, pressure variation factor (PVF) and maximum pressure (pm) are employed to directly compare the global and local effects, respectively, of variations in the geometry. The values of PVF and pm are then used to construct two Bayesian surrogate models. These models are statistically analysed to visualise how each geometric parameter influences PVF and pm. Percentage of stenosis is found to influence these pressure based metrics more than any other geometric parameter. Later, we identify bypass grafts with optimal geometric and material properties which have low values of PVF on five test cases with 70%, 75%, 80%, 85% and 90% stenosis in the ICA, respectively

The role of biomechanics in the assessment of carotid atherosclerosis severity: a numerical approach

Numerical fluid biomechanics has been proved to be an efficient tool for understanding vascular diseases including atherosclerosis. There are many evidences that atherosclerosis plaque formation and rupture are associated with blood flow behavior. In fact, zones of low wall shear stress are vivid areas of proliferation of atherosclerosis, and in particular, in the carotid artery. In this paper a model is presented for investigating how the presence of the plaque influences the distribution of the wall shear stress. In complement to a first approach with rigid walls, an FSI model is developed as well to simulate the coupling between the blood flow and the carotid artery deformation. The results show that the presence of the plaque causes an attenuation of the WSS in the after-plaque region as well as the emergence of recirculation areas

Computational modeling of low-density lipoprotein accumulation at the carotid artery bifurcation after stenting

Restenosis typically occurs in regions of low and oscillating wall shear stress, which also favor the accumulation of atherogenic macromolecules such as low-density lipoprotein (LDL). This study aims to evaluate LDL transport and accumulation at the carotid artery bifurcation following carotid artery stenting (CAS) by means of computational simulation. The computational model consists of coupled blood flow and LDL transport, with the latter being modeled as a dilute substance dissolved in the blood and transported by the flow through a convection-diffusion transport equation. The endothelial layer was assumed to be permeable to LDL, and the hydraulic conductivity of LDL was shear-dependent. Anatomically realistic geometric models of the carotid bifurcation were built based on pre- and post-stent computed tomography (CT) scans. The influence of stent design was investigated by virtually deploying two different types of stents (open- and closed-cell stents) into the same carotid bifurcation model. Predicted LDL concentrations were compared between the post-stent carotid models and the relatively normal contralateral model reconstructed from patient-specific CT images. Our results show elevated LDL concentration in the distal section of the stent in all post-stent models, where LDL concentration is 20 times higher than that in the contralateral carotid. Compared with the open-cell stents, the closed-cell stents have larger areas exposed to high LDL concentration, suggesting an increased risk of stent restenosis. This computational approach is readily applicable to multiple patient studies and, once fully validated against follow-up data, it can help elucidate the role of stent strut design in the development of in-stent restenosis after CAS

Quantification of Local Hemodynamic Alterations Caused by Virtual Implantation of Three Commercially Available Stents for the Treatment of Aortic Coarctation

Patients with coarctation of the aorta (CoA) are prone to morbidity including atherosclerotic plaque that has been shown to correlate with altered wall shear stress (WSS) in the descending thoracic aorta (dAo). We created the first patient-specific computational fluid dynamics (CFD) model of a CoA patient treated by Palmaz stenting to date, and compared resulting WSS distributions to those from virtual implantation of Genesis XD and modified NuMED CP stents, also commonly used for CoA. CFD models were created from magnetic resonance imaging, fluoroscopy and blood pressure data. Simulations incorporated vessel deformation, downstream vascular resistance and compliance to match measured data and generate blood flow velocity and time-averaged WSS (TAWSS) results. TAWSS was quantified longitudinally and circumferentially in the stented region and dAo. While modest differences were seen in the distal portion of the stented region, marked differences were observed downstream along the posterior dAo and depended on stent type. The Genesis XD model had the least area of TAWSS values exceeding the threshold for platelet aggregation in vitro, followed by the Palmaz and NuMED CP stents. Alterations in local blood flow patterns and WSS imparted on the dAo appear to depend on the type of stent implanted for CoA. Following confirmation in larger studies, these findings may aid pediatric interventional cardiologists in selecting the most appropriate stent for each patient, and ultimately reduce long-term morbidity following treatment for CoA by stenting

Computational Fluid Dynamics in Cardiovascular Disease

Computational fluid dynamics (CFD) is a mechanical engineering field for analyzing fluid flow, heat transfer, and associated phenomena, using computer-based simulation. CFD is a widely adopted methodology for solving complex problems in many modern engineering fields. The merit of CFD is developing new and improved devices and system designs, and optimization is conducted on existing equipment through computational simulations, resulting in enhanced efficiency and lower operating costs. However, in the biomedical field, CFD is still emerging. The main reason why CFD in the biomedical field has lagged behind is the tremendous complexity of human body fluid behavior. Recently, CFD biomedical research is more accessible, because high performance hardware and software are easily available with advances in computer science. All CFD processes contain three main components to provide useful information, such as pre-processing, solving mathematical equations, and post-processing. Initial accurate geometric modeling and boundary conditions are essential to achieve adequate results. Medical imaging, such as ultrasound imaging, computed tomography, and magnetic resonance imaging can be used for modeling, and Doppler ultrasound, pressure wire, and non-invasive pressure measurements are used for flow velocity and pressure as a boundary condition. Many simulations and clinical results have been used to study congenital heart disease, heart failure, ventricle function, aortic disease, and carotid and intra-cranial cerebrovascular diseases. With decreasing hardware costs and rapid computing times, researchers and medical scientists may increasingly use this reliable CFD tool to deliver accurate results. A realistic, multidisciplinary approach is essential to accomplish these tasks. Indefinite collaborations between mechanical engineers and clinical and medical scientists are essential. CFD may be an important methodology to understand the pathophysiology of the development and progression of disease and for establishing and creating treatment modalities in the cardiovascular field