In vitro, time-resolved PIV comparison of the effect of stent design on wall shear stress

Smartphone-based particle tracking velocimetry for the in vitro assessment of coronary flows

The present study adopts a smartphone-based approach for the experimental characterization of coronary flows. Technically, Particle Tracking Velocimetry (PTV) measurements were performed using a smartphone camera and a low-power continuous wave laser in realistic healthy and stenosed phantoms of left anterior descending artery with inflow Reynolds numbers approximately ranging from 20 to 200. A Lagrangian-Eulerian mapping was performed to convert Lagrangian PTV velocity data to a Eulerian grid. Eulerian velocity and vorticity data obtained from smartphone-based PTV measurements were compared with Particle Image Velocimetry (PIV) measurements performed with a smartphone-based setup and with a conventional setup based on a high-power double-pulsed laser and a CMOS camera. Smartphone-based PTV and PIV velocity flow fields substantially agreed with conventional PIV measurements, with the former characterized by lower average percentage differences than the latter. Discrepancies emerged at high flow regimes, especially at the stenosis throat, due to particle image blur generated by smartphone camera shutter speed and image acquisition frequency. In conclusion, the present findings demonstrate the feasibility of PTV measurements using a smartphone camera and a low-power light source for the in vitro characterization of cardiovascular flows for research, industrial and educational purposes, with advantages in terms of costs, safety and usability.

Smartphone-based particle image velocimetry for cardiovascular flows applications: A focus on coronary arteries

An experimental set-up is presented for the in vitro characterization of the fluid dynamics in personalized phantoms of healthy and stenosed coronary arteries. The proposed set-up was fine-tuned with the aim of obtaining a compact, flexible, low-cost test-bench for biomedical applications. Technically, velocity vector fields were measured adopting a so-called smart-PIV approach, consisting of a smartphone camera and a low-power continuous laser (30 mW). Experiments were conducted in realistic healthy and stenosed 3D-printed phantoms of left anterior descending coronary artery reconstructed from angiographic images. Time resolved image acquisition was made possible by the combination of the image acquisition frame rate of last generation commercial smartphones and the flow regimes characterizing coronary hemodynamics (velocities in the order of 10 cm/s). Different flow regimes (Reynolds numbers ranging from 20 to 200) were analyzed. The smart-PIV approach was able to provide both qualitative flow visualizations and quantitative results. A comparison between smart-PIV and conventional PIV (i.e., the gold-standard experimental technique for bioflows characterization) measurements showed a good agreement in the measured velocity vector fields for both the healthy and the stenosed coronary phantoms. Displacement errors and uncertainties, estimated by applying the particle disparity method, confirmed the soundness of the proposed smart-PIV approach, as their values fell within the same range for both smart and conventional PIV measured data (≈5% for the normalized estimated displacement error and below 1.2 pixels for displacement uncertainty). In conclusion, smart-PIV represents an easy-to-implement, low-cost methodology for obtaining an adequately robust experimental characterization of cardiovascular flows. The proposed approach, to be intended as a proof of concept, candidates to become an easy-to-handle test bench suitable for use also outside of research labs, e.g., for educational or industrial purposes, or as first-line investigation to direct and guide subsequent conventional PIV measurements.

Numerical and Experimental Analysis of Shear Stress Influence on Cellular Viability in Serpentine Vascular Channels

3D bioprinting has emerged as a tool for developing in vitro tissue models for studying disease progression and drug development. The objective of the current study was to evaluate the influence of flow driven shear stress on the viability of cultured cells inside the luminal wall of a serpentine network. Fluid-structure interaction was modeled using COMSOL Multiphysics for representing the elasticity of the serpentine wall. Experimental analysis of the serpentine model was performed on the basis of a desirable inlet flow boundary condition for which the most homogeneously distributed wall shear stress had been obtained from numerical study. A blend of Gelatin-methacryloyl (GelMA) and PEGDA200 PhotoInk was used as a bioink for printing the serpentine network, while facilitating cell growth within the pores of the gelatin substrate. Human umbilical vein endothelial cells were seeded into the channels of the network to simulate the blood vessels. A Live-Dead assay was performed over a period of 14 days to observe the cellular viability in the printed vascular channels. It was observed that cell viability increases when the seeded cells were exposed to the evenly distributed shear stresses at an input flow rate of 4.62 mm/min of the culture media, similar to that predicted in the numerical model with the same inlet boundary condition. It leads to recruitment of a large number of focal adhesion point nodes on cellular membrane, emphasizing the influence of such phenomena on promoting cellular morphologies.

An Oscillatory Shear Index-Based Model to Describe Progressive Carotid Artery Stenosis

Background and aims: This study describes and demonstrates the applicability of a novel in silico method for modeling progressive carotid artery stenosis using the oscillatory shear index (OSI) as the basis of stenosis. Methods: Three-dimensional reconstructions of 11 carotid arteries were generated using patient-derived magnetic resonance angiography and duplex ultrasound data. Computational fluid dynamic simulations were sequentially generated following computational stenosis assessment, and corresponding changes in OSI were observed and used as measure of morphological stabilization. Results: 6 carotid models showed progressive stenosis with statistically significant increases in regions of high OSI (OSI >.2, P < .05) with eventual carotid occlusion in 1 of the cases. Three models remained free or nearly free of increased OSI, whereas 1 model showed an overall decrease in high OSI regions (P < .05) and another trended in that direction but did not achieve statistical significance (P = .145). Conclusions: To our knowledge, this is the first computational model describing progressive stenosis in any peripheral artery including the carotid. Taken together, this study provides a novel framework for computational hemodynamic investigations on progressive atherosclerosis in the carotid artery.

Review of the Development of Hemodynamic Modeling Techniques to Capture Flow Behavior in Arteries Affected by Aneurysm, Atherosclerosis, and Stenting

Cardiovascular diseases (CVDs) are the leading cause of death in the developed world. CVD can include atherosclerosis, aneurysm, dissection, or occlusion of the main arteries. Many CVDs are caused by unhealthy hemodynamics. Some CVDs can be treated with the implantation of stents and stent grafts. Investigations have been carried out to understand the effects of stents and stent grafts have on arteries and the hemodynamic changes post-treatment. Numerous studies on stent hemodynamics have been carried out using computational fluid dynamics (CFD) which has yielded significant insight into the effect of stent mesh design on near-wall blood flow and improving hemodynamics. Particle image velocimetry (PIV) has also been used to capture behavior of fluids that mimic physiological hemodynamics. However, PIV studies have largely been restricted to unstented models or intra-aneurysmal flow rather than peri or distal stent flow behaviors. PIV has been used both as a standalone measurement method and as a comparison to validate the CFD studies. This article reviews the successes and limitations of CFD and PIV-based modeling methods used to investigate the hemodynamic effects of stents. The review includes an overview of physiology and relevant mechanics of arteries as well as consideration of boundary conditions and the working fluids used to simulate blood for each modeling method along with the benefits and limitations introduced.

In Vitro Flow Visualization in a Lactating Human Breast Model

Human milk extraction from the breast is affected by the infant's oral activities. Natural suckling by the infant includes both intraoral vacuum and peripheral oral compression during breastfeeding. However, the contribution of each of these motions to milk extraction at the outlet and at the duct bifurcations is unclear. In this work, we investigated the flow field in a lactating breast model considering bifurcated milk ducts and multiphase breast-infant interactions. A bio-inspired breastfeeding simulator device was utilized to mimic an infant's oral feeding mechanism during breastfeeding and extract the human milk-mimicking Fluid from the transparent and elastic lactating breast phantom during experiments. Using a particle image velocimetry system, we found that the oscillatory flow under vacuum pressure provides a higher velocity field at the outlet compared to that when an infant applies both vacuum and oral compression pressures. Additionally, the intraoral vacuum coordinated with the oral peripheral compression causes stronger vorticities and secondary flows at the adjunction of the bifurcated ducts than the vacuum-only case. Vacuum-only extraction yields an increase in flow velocity at the outlet and could be one of the reasons for nipple pain, whereas infant's oral activities on the breast generated more vortices in the milk duct adjunctions and might cause milk duct clogs. This phenomenon is rationalized based on the validation of a previous in vivo clinical study of milk production compared between commercial pumps and infant suckling. The fact that milk consumption of vacuum-only extraction is less than that of vacuum plus oral compression further explains the effectiveness of applying a natural suckling pattern in human lactation.

Wall shear stress mapping for human femoral artery based on ultrafast ultrasound vector Doppler estimations

PURPOSE

Wall shear stress (WSS), a type of friction exerted on the artery wall by flowing blood, is considered a crucial factor in atherosclerotic plaque development. Currently, achieving a reliable WSS mapping of an artery noninvasively by using existing imaging modalities is still challenging. In this study, a WSS mapping based on vector Doppler flow velocity estimation was proposed to measure the dynamic WSS on the human femoral artery.

METHODS

Because ultrafast ultrasound imaging was used here, flow-enhanced imaging was also performed to observe the moving blood flow condition. The performance of WSS mapping was verified using both straight (8 mm in diameter) and stenosis (70% of stenosis) phantoms under a pulsatile flow condition. A human study was conducted from five healthy volunteers.

RESULTS

Experimental results demonstrated that the WSS estimation was close to the standard value that was obtained from maximum velocity estimation in straight phantom experiments. In a stenosis phantom experiment, a low WSS region was observed at a site downstream of an obstruction, which is a high-risk area for plaque formation. Dynamic WSS mapping was accomplished in measurement in the femoral artery bifurcation. In measurements, the time-averaged WSS of the common femoral artery, superficial femoral artery, and deep femoral artery was 0.52± 0.19, 0.44 ± 0.21, and 0.29 ± 0.16 Pa, respectively, for the anterior wall and 0.29 ± 0.11, 0.54 ± 0.24, and 0.23 ± 0.10 Pa, respectively, for the posterior wall.

CONCLUSIONS

All results indicated that WSS mapping has the potential to be a useful tool for vessel duplex scanning in the future.

3D Printing for Cardiovascular Applications: From End-to-End Processes to Emerging Developments

3D printing as a means of fabrication has seen increasing applications in medicine in the last decade, becoming invaluable for cardiovascular applications. This rapidly developing technology has had a significant impact on cardiovascular research, its clinical translation and education. It has expanded our understanding of the cardiovascular system resulting in better devices, tools and consequently improved patient outcomes. This review discusses the latest developments and future directions of generating medical replicas ('phantoms') for use in the cardiovascular field, detailing the end-to-end process from medical imaging to capture structures of interest, to production and use of 3D printed models. We provide comparisons of available imaging modalities and overview of segmentation and post-processing techniques to process images for printing, detailed exploration of latest 3D printing methods and materials, and a comprehensive, up-to-date review of milestone applications and their impact within the cardiovascular domain across research, clinical use and education. We then provide an in-depth exploration of future technologies and innovations around these methods, capturing opportunities and emerging directions across increasingly realistic representations, bioprinting and tissue engineering, and complementary virtual and mixed reality solutions. The next generation of 3D printing techniques allow patient-specific models that are increasingly realistic, replicating properties, anatomy and function.

Integrating multi-fidelity blood flow data with reduced-order data assimilation

High-fidelity patient-specific modeling of cardiovascular flows and hemodynamics is challenging. Direct blood flow measurement inside the body with in-vivo measurement modalities such as 4D flow magnetic resonance imaging (4D flow MRI) suffer from low resolution and acquisition noise. In-vitro experimental modeling and patient-specific computational fluid dynamics (CFD) models are subject to uncertainty in patient-specific boundary conditions and model parameters. Furthermore, collecting blood flow data in the near-wall region (e.g., wall shear stress) with experimental measurement modalities poses additional challenges. In this study, a computationally efficient data assimilation method called reduced-order modeling Kalman filter (ROM-KF) was proposed, which combined a sequential Kalman filter with reduced-order modeling using a linear model provided by dynamic mode decomposition (DMD). The goal of ROM-KF was to overcome low resolution and noise in experimental and uncertainty in CFD modeling of cardiovascular flows. The accuracy of the method was assessed with 1D Womersley flow, 2D idealized aneurysm, and 3D patient-specific cerebral aneurysm models. Synthetic experimental data were used to enable direct quantification of errors using benchmark datasets. The accuracy of ROM-KF in reconstructing near-wall hemodynamics was assessed by applying the method to problems where near-wall blood flow data were missing in the experimental dataset. The ROM-KF method provided blood flow data that were more accurate than the computational and synthetic experimental datasets and improved near-wall hemodynamics quantification.

Parametric Optimization of 3D Printed Hydrogel-Based Cardiovascular Stent

PURPOSE

This study aimed to develop personalized biodegradable stent (BDS) for the treatment of coronary heart disease. Three-dimensional (3D) printing technique has offered easy and fast fabrication of BDS with enhanced reproducibility and efficacy.

METHODS

A variety of BDS were printed with 3 types of hydrogel (~5 ml) resources (10%w/v sodium alginate (SA), 10%w/v cysteine-sodium alginate (SA-CYS), and 10%w/v cysteine-sodium alginate with 0.4%w/v PLA-nanofibers (SA-CYS-NF)) dispersed from an 22G print head nozzle attached to the BD-syringe. The printability of hydrogels into 3D structures was examined based on such variables as hydrogel's viscosity, printing distance, printing speed and the nozzle size.

RESULTS

It was demonstrated that alginate composition (10%w/v) offered BDS with sufficient viscosity that defined the thickness and swelling ratio of the stent struts. The thickness of the strut was found to be 338.7 ± 29.3 μm, 262.5 ± 14.7 μm and 237.1 ± 14.7 μm for stents made of SA, SA-CYS and SA-CYS-NF, respectively. SA-CYS-NF stent displayed the highest swelling ratio of 38.8 ± 2.9% at the initial 30 min, whereas stents made of SA and SA-CYS had 23.1 ± 2.4% and 22.0 ± 2.4%, respectively.

CONCLUSION

The printed stents had sufficient mechanical strength and were stable against pseudo-physiological wall shear stress. An addition of nanofibers to alginate hydrogel significantly enhanced the biodegradation rates of the stents. In vitro cell culture studies revealed that stents had no cytotoxic effects on human umbilical vein endothelial cells (HUVECs) and Raw 264.7 cells (i.e., Monocyte/macrophage-like cells), supporting that stents are biocompatible and can be explored for future clinical applications.

PIV Analysis of Haemodynamics Distal to the Frozen Elephant Trunk Stent Surrogate

PURPOSE

The Frozen Elephant Trunk (FET) stent is a hybrid endovascular device that may be implemented in the event of an aneurysm or aortic dissection of the aortic arch or superior descending aorta. A Type 1B endoleak can lead to intrasaccular flow during systole and has been identified as a known failure of the FET stent graft. The purpose was to develop in-vitro modelling techniques to enable the investigation of the known failure.

METHODS

A silicone aortic phantom and 3D printed surrogate stent graft were manufactured to investigate the haemodynamics of a Type 1B endoleak. Physiological pulsatile flow dynamics distal to the surrogate stent graft were investigated in-vitro using Particle Image Velocimetry (PIV).

RESULTS

PIV captured recirculation zones and an endoleak distal to the surrogate stent graft. The endoleak was developed at the peak of systole and sustained until the onset of diastole. The endoleak was asymmetric, indicating a potential variation in the phantom artery wall thickness or stent alignment. Recirculation was identified on the posterior dorsal line during late systole.

CONCLUSIONS

The identification of the Type 1B endoleak proved that in-vitro modelling can be used to investigate complex compliance changes and wall motions. The recirculation may indicate the potential for long term intimal layer inflammatory issues such as atherosclerosis. These results may aid future remediation techniques or stent design. Further development of the methods used in this experiment may assist with the future testing of stents prior to animal or human trial.

Influence of near-wall PIV data on recirculation hemodynamics in a patient-specific moderate stenosis: Experimental-numerical comparison

BACKGROUND

Recirculation zones within the blood vessels are known to influence the initiation and progression of atherosclerotic lesions. Quantification of recirculation parameters with accuracy remains subjective due to uncertainties in measurement of velocity and derived wall shear stress (WSS).

OBJECTIVE

The primary aim is to determine recirculation height and length from PIV experiments while validating with two different numerical methods: finite-element (FE) and -volume (FV). Secondary aim is to analyze how FE and FV compare within themselves.

METHODS

PIV measurements were performed to obtain velocity profiles at eight cross sections downstream of stenosis at flow rate of 200 ml/min. WSS was obtained by linear/quadratic interpolation of experimental velocity measurements close to wall.

RESULTS

Recirculation length obtained from PIV technique was 1.47 cm and was within 2.2% of previously reported in-vitro measurements. Derived recirculation length from PIV agreed within 6.8% and 8.2% of the FE and FV calculations, respectively. For lower shear rate, linear interpolation with five data points results in least error. For higher shear rate either higher order (quadratic) interpolation with five data points or lower order (linear) with lesser (three) data points leads to better results.

CONCLUSION

Accuracy of the recirculation parameters is dependent on number of near wall PIV data points and the type of interpolation algorithm used.

Determining Haemodynamic Wall Shear Stress in the Rabbit Aorta In Vivo Using Contrast-Enhanced Ultrasound Image Velocimetry

Abnormal blood flow and wall shear stress (WSS) can cause and be caused by cardiovascular disease. To date, however, no standard method has been established for mapping WSS in vivo. Here we demonstrate wide-field assessment of WSS in the rabbit abdominal aorta using contrast-enhanced ultrasound image velocimetry (UIV). Flow and WSS measurements were made independent of beam angle, curvature or branching. Measurements were validated in an in silico model of the rabbit thoracic aorta with moving walls and pulsatile flow. Mean errors over a cardiac cycle for velocity and WSS were 0.34 and 1.69%, respectively. In vivo time average WSS in a straight segment of the suprarenal aorta correlated highly with simulations (PC = 0.99) with a mean deviation of 0.29 Pa or 5.16%. To assess fundamental plausibility of the measurement, UIV WSS was compared to an analytic approximation derived from the Poiseuille equation; the discrepancy was 17%. Mapping of WSS was also demonstrated in regions of arterial branching. High time average WSS (TAWSS_xz = 3.4 Pa) and oscillatory flow (OSI_xz = 0.3) were observed near the origin of conduit arteries. In conclusion, we have demonstrated that contrast-enhanced UIV is capable of measuring spatiotemporal variation in flow velocity, arterial wall location and hence WSS in vivo with high accuracy over a large field of view.

PIV Analysis of Stented Haemodynamics in the Descending Aorta

Cardiovascular diseases (CVD) are the leading cause of death in the developed world and aortic aneurysm is a key contributor. Aortic aneurysms typically occur in the thoracic aorta and can extend into the descending aorta. The Frozen Elephant Trunk stent (FET) is one of the leading treatments for the aneurysms extending into the descending aorta. This study focuses on the in-vitro experimentation of a stented descending aorta, investigating the haemodynamics in a compliant phantom. A silicone phantom of the descending aorta was manufactured using a lost core casting method. A PVC stent was manufactured using the same mould core. Particle Image Velocimetry (PIV) was used for pulsatile studies, focusing specifically on the passive fixation at the distal end of the FET. The results showed an apparent expansion in the diastolic period that was identified to be a collapse in the lateral plane. Flow recirculation regions were identified during the collapse. The collapse was attributed to low upstream and high downstream pressures causing a vacuum effect. The findings may imply a potential risk introduced by the FET stent that requires further investigation.

Evaluating medical device and material thrombosis under flow: current and emerging technologies

This review highlights the importance of flow in medical device thrombosis and explores current and emerging technologies to evaluate dynamic biomaterial Thrombosis in vitro.

Evaluation of a Desktop 3D Printed Rigid Refractive-Indexed-Matched Flow Phantom for PIV Measurements on Cerebral Aneurysms

PURPOSE

Fabrication of a suitable flow model or phantom is critical to the study of biomedical fluid dynamics using optical flow visualization and measurement methods. The main difficulties arise from the optical properties of the model material, accuracy of the geometry and ease of fabrication.

METHODS

Conventionally an investment casting method has been used, but recently advancements in additive manufacturing techniques such as 3D printing have allowed the flow model to be printed directly with minimal post-processing steps. This study presents results of an investigation into the feasibility of fabrication of such models suitable for particle image velocimetry (PIV) using a common 3D printing Stereolithography process and photopolymer resin.

RESULTS

An idealised geometry of a cerebral aneurysm was printed to demonstrate its applicability for PIV experimentation. The material was shown to have a refractive index of 1.51, which can be refractive matched with a mixture of de-ionised water with ammonium thiocyanate (NH₄SCN). The images were of a quality that after applying common PIV pre-processing techniques and a PIV cross-correlation algorithm, the results produced were consistent within the aneurysm when compared to previous studies.

CONCLUSIONS

This study presents an alternative low-cost option for 3D printing of a flow phantom suitable for flow visualization simulations. The use of 3D printed flow phantoms reduces the complexity, time and effort required compared to conventional investment casting methods by removing the necessity of a multi-part process required with investment casting techniques.

Benchmark for Numerical Models of Stented Coronary Bifurcation Flow

In-stent restenosis ails many patients who have undergone stenting. When the stented artery is a bifurcation, the intervention is particularly critical because of the complex stent geometry involved in these structures. Computational fluid dynamics (CFD) has been shown to be an effective approach when modeling blood flow behavior and understanding the mechanisms that underlie in-stent restenosis. However, these CFD models require validation through experimental data in order to be reliable. It is with this purpose in mind that we performed particle image velocimetry (PIV) measurements of velocity fields within flows through a simplified coronary bifurcation. Although the flow in this simplified bifurcation differs from the actual blood flow, it emulates the main fluid dynamic mechanisms found in hemodynamic flow. Experimental measurements were performed for several stenting techniques in both steady and unsteady flow conditions. The test conditions were strictly controlled, and uncertainty was accurately predicted. The results obtained in this research represent readily accessible, easy to emulate, detailed velocity fields and geometry, and they have been successfully used to validate our numerical model. These data can be used as a benchmark for further development of numerical CFD modeling in terms of comparison of the main flow pattern characteristics.

Effect of pannus formation on the prosthetic heart valve: In vitro demonstration using particle image velocimetry

Although hemodynamic influence of the subprosthetic tissue, termed as pannus, may contribute to prosthetic aortic valve dysfunction, the relationship between pannus extent and hemodynamics in the prosthetic valve has rarely been reported. We investigated the fluid dynamics of pannus formation using in vitro experiments with particle image velocimetry. Subvalvular pannus formation caused substantial changes in prosthetic valve transvalvular peak velocity, transvalvular pressure gradient (TPG) and opening angle. Maximum flow velocity and corresponding TPG were mostly affected by pannus width. When the pannus width was 25% of the valve diameter, pannus formation elevated TPG to >2.5 times higher than that without pannus formation. Opening dysfunction was observed only for a pannus involvement angle of 360°. Although circumferential pannus with an involvement angle of 360° decreased the opening angle of the valve from approximately 82° to 58°, eccentric pannus with an involvement angle of 180° did not induce valve opening dysfunction. The pannus involvement angle largely influenced the velocity flow field at the aortic sinus and corresponding hemodynamic indices, including wall shear stress, principal shear stress and viscous energy loss distributions. Substantial discrepancy between the velocity-based TPG estimation and direct pressure measurements was observed for prosthetic valve flow with pannus formation.

Spatio-Temporal Flow and Wall Shear Stress Mapping Based on Incoherent Ensemble-Correlation of Ultrafast Contrast Enhanced Ultrasound Images

In this study, a technique for high-frame-rate ultrasound imaging velocimetry (UIV) is extended first to provide more robust quantitative flow velocity mapping using ensemble correlation of images without coherent compounding, and second to generate spatio-temporal wall shear stress (WSS) distribution. A simulation model, which couples the ultrasound simulator with analytical flow solution, was implemented to evaluate its accuracy. It is shown that the proposed approach can reduce errors in velocity estimation by up to 10-fold in comparison with the coherent correlation approach. Mean errors (ME) of 3.2% and 8.6% were estimated under a steady flow condition, while 3.0% and 10.6% were found under a pulsatile condition for the velocity and wall shear rate (WSR) measurement, respectively. Appropriate filter parameters were selected to constrain the velocity profiles before WSR estimations and the effects of incorrect wall tracking were quantified under a controlled environment. Although accurate wall tracking is found to be critical in WSR measurement (as a 200 µm deviation from the wall may yield up to a 60% error), this can be mitigated by HFR imaging (of up to 10 kHz) with contrast agents, which allow for improved differentiation of the wall-fluid boundaries. In vitro investigations on two carotid bifurcation phantoms, normal and diseased, were conducted, and their relative differences in terms of the flow patterns and WSR distribution were demonstrated. It is shown that high-frame-rate UIV technique can be a non-invasive tool to measure quantitatively the spatio-temporal velocity and WSS distribution.

The stentable in vitro artery: an instrumented platform for endovascular device development and optimization

Although vascular disease is a leading cause of mortality, in vitro tools for controlled, quantitative studies of vascular biological processes in an environment that reflects physiological complexity remain limited. We developed a novel in vitro artery that exhibits a number of unique features distinguishing it from tissue-engineered or organ-on-a-chip constructs, most notably that it allows deployment of endovascular devices including stents, quantitative real-time tracking of cellular responses and detailed measurement of flow velocity and lumenal shear stress using particle image velocimetry. The wall of the stentable in vitro artery consists of an annular collagen hydrogel containing smooth muscle cells (SMCs) and whose lumenal surface is lined with a monolayer of endothelial cells (ECs). The system has in vivo dimensions and physiological flow conditions and allows automated high-resolution live imaging of both SMCs and ECs. To demonstrate proof-of-concept, we imaged and quantified EC wound healing, SMC motility and altered shear stresses on the endothelium after deployment of a coronary stent. The stentable in vitro artery provides a unique platform suited for a broad array of research applications. Wide-scale adoption of this system promises to enhance our understanding of important biological events affecting endovascular device performance and to reduce dependence on animal studies.

Endoleak Assessment Using Computational Fluid Dynamics and Image Processing Methods in Stented Abdominal Aortic Aneurysm Models

Endovascular aortic aneurysm repair (EVAR) is a predominant surgical procedure to reduce the risk of aneurysm rupture in abdominal aortic aneurysm (AAA) patients. Endoleak formation, which eventually requires additional surgical reoperation, is a major EVAR complication. Understanding the etiology and evolution of endoleak from the hemodynamic perspective is crucial to advancing the current posttreatments for AAA patients who underwent EVAR. Therefore, a comprehensive flow assessment was performed to investigate the relationship between endoleak and its surrounding pathological flow fields through computational fluid dynamics and image processing. Six patient-specific models were reconstructed, and the associated hemodynamics in these models was quantified three-dimensionally to calculate wall stress. To provide a high degree of clinical relevance, the mechanical stress distribution calculated from the models was compared with the endoleak positions identified from the computed tomography images of patients through a series of imaging processing methods. An endoleak possibly forms in a location with high local wall stress. An improved stent graft (SG) structure is conceived accordingly by increasing the mechanical strength of the SG at peak wall stress locations. The presented analytical paradigm, as well as numerical analysis using patient-specific models, may be extended to other common human cardiovascular surgeries.

Hemodynamics of Stent Implantation Procedures in Coronary Bifurcations: An In Vitro Study

Stent implantation in coronary bifurcations presents unique challenges and currently there is no universally accepted stent deployment approach. Despite clinical and computational studies, the effect of each stent implantation method on the coronary artery hemodynamics is not well understood. In this study the hemodynamics of stented coronary bifurcations under pulsatile flow conditions were investigated experimentally. Three implantation methods, provisional side branch (PSB), culotte (CUL), and crush (CRU), were investigated using time-resolved particle image velocimetry to measure the velocity fields. Subsequently, hemodynamic parameters including wall shear stress, oscillatory shear index (OSI), and relative residence time (RRT) were calculated. The pressure field through the vessel was non-invasively quantified and pressure wave speeds were computed. The effects of each stented case were evaluated and compared against an un-stented case. CRU provided the lowest compliance mismatch, but demonstrated detrimental stent interactions. PSB, the clinically preferred method, and CUL maintained many normal flow conditions. However, PSB provided about a 300% increase in both OSI and RRT. CUL yielded a 10 and 85% increase in OSI and RRT, respectively. The results of this study support the concept that different bifurcation stenting techniques result in hemodynamic environments that deviate from that of un-stented bifurcations, to varying degrees.

Local blood flow patterns in stented coronary bifurcations: an experimental and numerical study

PURPOSE

Despite the atheroprone environment of blood flow in coronary bifurcations, limited quantitative information is available on the hemodynamics occurring in these geometries, both before and after their treatment with endovascular stents. Previous studies have focused on computational fluid dynamics (CFD) analyses and have bypassed the challenges associated with experimentally representing the flow environment, providing no means for validation. This study analyzed steady flow conditions in 3 bifurcation angles and 4 different single- and double-stenting procedures, which are used clinically in coronary bifurcations.

METHODS

The numerical aspect of this study utilized geometries derived from CAD models (nonstented cases) and finite element simulations (stented cases). Digital particle image velocimetry (DPIV) testing was conducted within compliant bifurcating models for which an uncertainty analysis was performed at each measurement location for CFD validation purposes. Results were analyzed in terms of velocity magnitude contour maps and axial velocity profiles at several locations in the bifurcated vessels.

RESULTS AND CONCLUSIONS

Qualitatively, the 2 approaches showed agreement in the bulk flow patterns. However, the velocity computed with CFD was outside the DPIV uncertainty estimates, which can be attributed to the intrinsic difference and modeling assumptions of the 2 approaches. The findings reveal that wider bifurcation angles and double-stenting procedures are both characterized by increased areas of low flow and recirculation. Additionally, inferior performance in terms of viscous and wall shear stresses was observed in double-stented cases.

Effects of elastic modulus change in helical tubes under the influence of dynamic changes in curvature and torsion

The incidence of stent late restenosis is high (Zwart et al., 2010, "Coronary Stent Thrombosis in the Current Era: Challenges and Opportunities for Treatment," Curr. Treat. Options Cardiovasc. Med., 12(1), pp. 46-57) despite the extensive use of stents, and is most prevalent at the proximal and distal ends of the stent. Elastic modulus change in stented coronary arteries subject to the motion of the myocardium is not studied extensively. It is our objective to understand and reveal the mechanism by which changes in elastic modulus and geometry contribute to the generation of nonphysiological wall shear stress (WSS). Such adverse hemodynamic conditions could have an effect on the onset of restenosis. Three-dimensional (3D), spatiotemporally resolved computational fluid dynamics (CFD) simulations of pulsatile flow with moving wall boundaries and fluid structure interaction (FSI) were carried out for a helical artery with physiologically relevant flow parameters. To study the effect of coronary artery (CA) geometry change on stent elastic modulus mismatch, models where the curvature, torsion and both curvature and torsion change were examined. The elastic modulus is increased by a factor of two, five, and ten in the stented section for all three modes of motion. The changes in elastic modulus and arterial geometry cause critical variations in the local pressure and velocity gradients and secondary flow patterns. The pressure gradient change is 47%, with respect to the unstented baseline when the elastic modulus is increased to 10. The corresponding WSS change is 15.4%. We demonstrate that these changes are attributed to the production of vorticity (vorticity flux) caused by the wall movement and elastic modulus discontinuity. The changes in curvature dominate torsion changes in terms of the effects to local hemodynamics. The elastic modulus discontinuities along with the dynamic change in geometry affected the secondary flow patterns and vorticity flux, which in turn affects the WSS.

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.

Sequential Structural and Fluid Dynamic Numerical Simulations of a Stented Bifurcated Coronary Artery

Despite their success, stenting procedures are still associated to some clinical problems like sub-acute thrombosis and in-stent restenosis. Several clinical studies associate these phenomena to a combination of both structural and hemodynamic alterations caused by stent implantation. Recently, numerical models have been widely used in the literature to investigate stenting procedures but always from either a purely structural or fluid dynamic point of view. The aim of this work is the implementation of sequential structural and fluid dynamic numerical models to provide a better understanding of stenting procedures in coronary bifurcations. In particular, the realistic geometrical configurations obtained with structural simulations were used to create the fluid domains employed within transient fluid dynamic analyses. This sequential approach was applied to investigate the final kissing balloon (FKB) inflation during the provisional side branch technique. Mechanical stresses in the arterial wall and the stent as well as wall shear stresses along the arterial wall were examined before and after the FKB deployment. FKB provoked average mechanical stresses in the arterial wall almost 2.5 times higher with respect to those induced by inflation of the stent in the main branch only. Results also enlightened FKB benefits in terms of improved local blood flow pattern for the side branch access. As a drawback, the FKB generates a larger region of low wall shear stress. In particular, after FKB the percentage of area characterized by wall shear stresses lower than 0.5 Pa was 79.0%, while before the FKB it was 62.3%. For these reasons, a new tapered balloon dedicated to bifurcations was proposed. The inclusion of the modified balloon has reduced the mechanical stresses in the proximal arterial vessel to 40% and the low wall shear stress coverage area to 71.3%. In conclusion, these results show the relevance of the adopted sequential approach to study the wall mechanics and the hemodynamics created by stent deployment.

Drug-eluting stent restenosis: effect of drug type, release kinetics, hemodynamics and coating strategy

Restenosis following stent implantation diminishes the procedure's efficacy influencing long-term clinical outcomes. Stent-based drug delivery emerged a decade ago as an effective means of reducing neointimal hyperplasia by providing localized pharmacotherapy during the acute phase of the stent-induced injury and the ensuing pathobiological mechanisms. However, drug-eluting stent (DES) restenosis may still occur especially when stents are used in complex anatomical and clinical scenarios. A DES consists of an intravascular metallic frame and carriers which allow controlled release of active pharmaceutical agents; all these components are critical in determining drug distribution locally and thus anti-restenotic efficacy. Furthermore, dynamic flow phenomena characterizing the vascular environment, and shear stress distribution, are greatly influenced by stent implantation and play a significant role in drug deposition and bioavailability within local vascular tissue. In this review, we discuss the performance of DES and the interaction of the different DES components with the hemodynamic milieu emphasizing on the inhibition of clinical restenosis.

Multilaboratory particle image velocimetry analysis of the FDA benchmark nozzle model to support validation of computational fluid dynamics simulations

This study is part of a FDA-sponsored project to evaluate the use and limitations of computational fluid dynamics (CFD) in assessing blood flow parameters related to medical device safety. In an interlaboratory study, fluid velocities and pressures were measured in a nozzle model to provide experimental validation for a companion round-robin CFD study. The simple benchmark nozzle model, which mimicked the flow fields in several medical devices, consisted of a gradual flow constriction, a narrow throat region, and a sudden expansion region where a fluid jet exited the center of the nozzle with recirculation zones near the model walls. Measurements of mean velocity and turbulent flow quantities were made in the benchmark device at three independent laboratories using particle image velocimetry (PIV). Flow measurements were performed over a range of nozzle throat Reynolds numbers (Re(throat)) from 500 to 6500, covering the laminar, transitional, and turbulent flow regimes. A standard operating procedure was developed for performing experiments under controlled temperature and flow conditions and for minimizing systematic errors during PIV image acquisition and processing. For laminar (Re(throat)=500) and turbulent flow conditions (Re(throat)≥3500), the velocities measured by the three laboratories were similar with an interlaboratory uncertainty of ∼10% at most of the locations. However, for the transitional flow case (Re(throat)=2000), the uncertainty in the size and the velocity of the jet at the nozzle exit increased to ∼60% and was very sensitive to the flow conditions. An error analysis showed that by minimizing the variability in the experimental parameters such as flow rate and fluid viscosity to less than 5% and by matching the inlet turbulence level between the laboratories, the uncertainties in the velocities of the transitional flow case could be reduced to ∼15%. The experimental procedure and flow results from this interlaboratory study (available at http://fdacfd.nci.nih.gov) will be useful in validating CFD simulations of the benchmark nozzle model and in performing PIV studies on other medical device models.

Hydrodynamic effects of compliance mismatch in stented arteries

Cardiovascular diseases are the number one cause of death in the world, making the understanding of hemodynamics and development of treatment options imperative. The most common modality for treatment of occlusive coronary artery diseases is the use of stents. Stent design profoundly influences the postprocedural hemodynamic and solid mechanical environment of the stented artery. However, despite their wide acceptance, the incidence of stent late restenosis is still high (Zwart et al., 2010, "Coronary Stent Thrombosis in the Current Era: Challenges and Opportunities for Treatment," Current Treatment Options in Cardiovascular Medicine, 12(1), pp. 46-57), and it is most prevailing at the proximal and distal ends of the stent. In this work, we focus our investigation on the localized hemodynamic effects of compliance mismatch due to the presence of a stent in an artery. The compliance mismatch in a stented artery is maximized at the proximal and distal ends of the stent. Hence, it is our objective to understand and reveal the mechanism by which changes in compliance contribute to the generation of nonphysiological wall shear stress (WSS). Such adverse hemodynamic conditions could have an effect on the onset of restenosis. Three-dimensional, spatiotemporally resolved computational fluid dynamics simulations of pulsatile flow with fluid-structure interaction were carried out for a simplified coronary artery with physiologically relevant flow parameters. A model with uniform elastic modulus is used as the baseline control case. In order to study the effect of compliance variation on local hemodynamics, this baseline model is compared with models where the elastic modulus was increased by two-, five-, and tenfold in the middle of the vessel. The simulations provided detailed information regarding the recirculation zone dynamics formed during flow reversals. The results suggest that discontinuities in compliance cause critical changes in local hemodynamics, namely, altering the local pressure and velocity gradients. The change in pressure gradient at the discontinuity was as high as 90%. The corresponding changes in WSS and oscillatory shear index calculated were 9% and 15%, respectively. We demonstrate that these changes are attributed to the physical mechanism associating the pressure gradient discontinuities to the production of vorticity (vorticity flux) due to the presence of the stent. The pressure gradient discontinuities and augmented vorticity flux are affecting the wall shear stresses. As a result, this work reveals how compliance variations act to modify the near wall hemodynamics of stented arteries.