Stent implants in stenosed arterial curvatures are prone to “Type IV” failures involving the complete transverse fracture of stents and linear displacement of the fractured parts. We present a protocol for detection of secondary flow (vortical) structures in a curved artery model, downstream of clinically relevant “Type IV” stent failures.
The arterial network in the human vasculature comprises of ubiquitously present blood vessels with complex geometries (branches, curvatures and tortuosity). Secondary flow structures are vortical flow patterns that occur in curved arteries due to the combined action of centrifugal forces, adverse pressure gradients and inflow characteristics. Such flow morphologies are greatly affected by pulsatility and multiple harmonics of physiological inflow conditions and vary greatly in size-strength-shape characteristics compared to non-physiological (steady and oscillatory) flows 1 – 7.
Secondary flow structures may ultimately influence the wall shear stress and exposure time of blood-borne particles toward progression of atherosclerosis, restenosis, sensitization of platelets and thrombosis 4 – 6, 8 – 13. Therefore, the ability to detect and characterize these structures under laboratory-controlled conditions is precursor to further clinical investigations.
A common surgical treatment to atherosclerosis is stent implantation, to open up stenosed arteries for unobstructed blood flow. But the concomitant flow perturbations due to stent installations result in multi-scale secondary flow morphologies 4 – 6. Progressively higher order complexities such as asymmetry and loss in coherence can be induced by ensuing stent failures vis-à-vis those under unperturbed flows 5. These stent failures have been classified as "Types I-to-IV" based on failure considerations and clinical severity 14.
This study presents a protocol for the experimental investigation of the complex secondary flow structures due to complete transverse stent fracture and linear displacement of fractured parts ("Type IV") in a curved artery model. The experimental method involves the implementation of particle image velocimetry (2C-2D PIV) techniques with an archetypal carotid artery inflow waveform, a refractive index matched blood-analog working fluid for phase-averaged measurements 15 – 18. Quantitative identification of secondary flow structures was achieved using concepts of flow physics, critical point theory and a novel wavelet transform algorithm applied to experimental PIV data 5, 6, 19 – 26.
Secondary flow structures are vortical flow patterns that occur in internal flow geometries with curvatures such as curved pipes and channels. These vortical structures arise due to the combined action of centrifugal forces, adverse pressure gradients and inflow characteristics. In general, secondary flow structures appear in planar cross-sections of curved pipes as symmetric Dean-type vortices under steady inflow and, symmetric Dean- and Lyne-type vortices under oscillatory inflow conditions 1 – 3. Secondary flow morphologies are greatly affected by pulsatility and multiple harmonics of pulsatile, physiological inflow conditions. These structures acquire markedly different size-strength-shape characteristics compared to non-physiological (steady and oscillatory) flows 1 – 6. Atherosclerotic lesion development in arteries is affected by the existence of high frequency shear oscillations in regions experiencing low mean shear 27, 28. Secondary flow structures may influence the progress of diseases such as atherosclerosis and possibly, mediate the endothelial response due to pulsatile blood flow by altering wall shear stresses and exposure times of blood-borne particles.
A common treatment to atherosclerosis, a complication resulting in narrowing of arteries by obstructive lesions, is the implantation of stents. Stent fractures are structural failures of implanted stents that lead to further medical complications such as in-stent restenosis (ISR), stent thrombosis and aneurysm formation 9 – 13. Stent fractures have been categorized into various failure "Types I-to-IV", wherein "Type IV" characterizes the highest clinical severity and is defined as the complete transverse fracture of stent struts along with linear displacements of the stent fragments 14. The protocol presented in this study describes an experimental method of visualization of secondary flow structures downstream of an idealized "Type IV" stent fracture in a curved artery model.
The suggested protocol has the following four essential features:
Design and fabrication of laboratory-scale stent models: Geometric description of stents can be associated with a set of self-expandable spirals (springs or helices) intertwined using Nitinol (an alloy of nickel and titanium) wires 29. The length of the stent and its strut diameter depend on the length scale of arterial lesions encountered during clinical implantation 5. Parametric variation of strut diameter and the rising of the winding (or pitch) leads to stents of various geometric configurations. A summary of stent design parameters chosen for 3D printing are presented in Table 1.
Preparation of a blood analog working fluid matched with kinematic viscosity of blood and refractive index of the test section: Optical access to the curved artery test section is required in order to make non-invasive velocity measurements. Accordingly, a Newtonian blood-mimicking working fluid with the refractive index of the vascular model and ideally, a dynamic viscosity, matching human blood is used to obtain accurate blood flow measurements 16 – 18, 30. The working fluid used in this study was reported by Deutsch et al. (2006), that comprised of 79% saturated aqueous sodium iodide (NaI), 20% pure glycerol, and 1% water (by volume)16.
Experimental arrangement for the detection of coherent secondary flow structures using a two-component, two-dimensional particle image velocimetry (2C-2D PIV): Experiments were designed to acquire phase-averaged secondary flow velocity data at various planar cross-sectional locations downstream of a combination of straight and curved stent sections embodying an idealized "Type IV" stent fracture 5, 6, 9, 14. The protocol-steps pertaining to the acquisition of secondary flow velocity fields using particle image velocimetry (PIV) technique involves a PIV system that comprises of a laser (light sheet) source, optics to focus and illuminate the regions of flow, a special cross-correlation charge coupled device (CCD-sensor or camera) and tracer particles to be illuminated by the light sheet within a short time interval (Δt; see Table 4) 31, 32.
The steps in the protocol assume the following: First, a calibrated, experimental set-up of a two-component, two-dimensional (2C-2D) PIV system that evaluates images by double-frame, single-exposure recordings. Second, the 2C-2D PIV system calculates the mean displacements of tracer particles by performing cross-correlation between two image frames acquired during each recording. A brief summary of PIV specifications and image acquisition software is presented in the materials and equipment table. Third, all safety precautions needed to operate the laser are followed by trained laboratory personnel according to the guidelines provided by the host institution. The authors suggest Refs. 31 and 32 for a holistic understanding of the implementation, functionality and application of PIV technique in aero-, hydro- and microfluid dynamics, correlation peak detection and displacement estimation, material and density of tracer particles and, measurement noise and accuracy. Also note that the laser and camera can be controlled by the PIV data acquisition computer (Figure 3A) and data processing software.
Data acquisition and post-processing for coherent structure detection: Phase-averaged secondary flow velocity measurements using a 2C-2D PIV were generated using the protocol description that follows. Post-processing of the data involved coherent secondary flow structure detection using the following three methods: continuous wavelet transforms, 5, 6, 19 – 24, 26.
The authors note that the velocity gradient tensor is essentially, a 3 x 3 matrix,
.
The protocol presents a method of acquiring two-dimensional experimental measurements (from 2C-2D PIV technique). Therefore, full experimental access to the velocity gradient tensor will not be attainable using this method. The velocity gradient tensor for each pixel of the PIV image should be a 2 x 2 matrix, . The z-component vorticity for each pixel is computed using the anti-symmetric part of the velocity gradient tensor . The result will be a 2D array of vorticity that can be visualized in a contour plot. The authors strongly suggest Ref. 25 for an eloquent discussion experimental access to the velocity gradient tensor toward enhancing the knowledge of vorticity dissipation, strain rates and coherent structure detection. Furthermore, the authors do not attempt to explore the inter-relationships between the aforementioned coherent structure detection methods and suggest Ref. 23, 24 for a comprehensive discussion on that subject.
The focus of the steps in the protocol is the quantitative identification of secondary flow (vortical) structures (also known as coherent structures). Three methods of coherent structure detection viz., and wavelet transformed vorticity are applied to velocity field data toward detection of multi-scale, multi-strength occurrences of secondary flow structures downstream of the idealized "Type IV" stent fracture.
The , defines a vortex as a spatial region where the Euclidean norm of the vorticity tensor dominates that of the rate of strain 19, 23, 24.The velocity gradient matrix is decomposed into symmetric (strain rate) and anti-symmetric (rotation) parts. Eigenvalues of strain rate matrix are computed; . Norm of the strain rate is then calculated; . Vorticity is computed from the anti- symmetric part. Enstrophy or square of z-component vorticity, ) is then computed. The is finally computed; . A contour plot of the entire set of with iso-regions of , will indicate secondary flow structures 19.
The , also known as 'swirling strength' is a vortex identification method performed by critical-point analysis of the local velocity gradient tensor and its corresponding eigenvalues 20 – 24. Eigenvalues of the velocity gradient tensor at each pixel are computed. The eigenvalues should be of the form, . A contour plot of with iso-regions of will indicate secondary flow structures 20 – 22.
Wavelet transform method utilizes an analyzing function (or wavelet) that has smoothness in physical and spectral spaces, is admissible (or has zero mean) and has a finite 5, 6, 26. By convolving a dilated or contracted wavelet with a 2D vorticity field, wavelet transformed vorticity field is generated comprising of coherent structures with a wide range of scales and strengths 5, 6, 26. Shannon entropy of the 2D wavelet-transformed vorticity field is computed to estimate the optimal wavelet scale at which all the coherent structures are adequately resolved. This entropy estimation involves a set of probabilities for each pixel such that , the normalized square modulus of the vorticity associated with the pixel at location m, n 5, 6. The procedural steps are presented graphically in Figure 6. The restrictions placed on the choice of the wavelet are presented in detail in Ref. 26. This protocol step describes the procedure for coherent structure detection using a 2D Ricker wavelet. The justification for the use of this wavelet for vortical pattern matching is presented in Ref. 5, 6 and the pertinent references cited therein.
1. Design and Fabrication of Stent Models
Note: The following steps have been followed to create laboratory-scale models of straight and curved stents. The installation of the two stent models will embody a "Type IV" fracture (fragmentation and linear displacement of fractured stent parts).
Note: The authors used Pro/Engineer software at the time of the research for creating CAD models of the stent geometry. The procedure below is generalized and may not include terms generic to the CAD software used. Other CAD packages available can also be used. The steps that follow are applicable for the CAD software that the authors used at the time of the research and have been adapted from the manufacturer's website. For further description of the rapid prototyping machine used by the authors see the Materials List. The parametric equations and initialized values for stent design are presented in Table 1 and Figure 1D and 1E are examples of the straight and curved stent models after rapid prototyping.
2. Prepare Kinematic Viscosity– and Refractive Index–matched Blood-analog Fluid
Note: The following procedure will yield approximately 600 ml of blood-analog solution. A summary of the chemical reagents and solvents with relevant properties used in the solution preparation are presented in the Materials List. Relevant material properties, suggested laboratory equipment and the guidelines for volumetric calculations are presented in Tables 2, 3 and 4, respectively.
3. Arrange Experiment for Measurement of Secondary Flow Velocity Fields Downstream of a "Type IV" Stent Failure
Note: The 180° curved artery test section comprises of two acrylic blocks glued together, 180° curved channel machined on each block and provision for inlet and outlet pipes as shown in Figures 1F, 3A and 3B. The materials selected for the test section are ensured to have optical accessibility. The authors report the refractive index of material used in the test section as ≈1.4914 5, 6 (See Table 2).
4. Acquire of Secondary Flow Velocity Fields
Note: The following description in the protocol pertains to the acquisition of secondary flow velocity fields using particle image velocimetry (PIV) technique. Figure 3B (schematic drawing) shows that there four locations (45°, 90°, 135° and 180°) with angular notches to facilitate laser sheet projection and making planar cross-sectional secondary flow velocity. The protocol steps pertain to measurements acquired for the 90° location. If the laser sheet is placed at 45° location, camera is placed at the 135° location to gain optical access for secondary flow measurements at the 45° location.
Note: The procedure below is generalized and may not include terms generic to the image acquisition and post processing software and the instrument control software used (see Materials List). Other image and data acquisition packages available can also be used in the protocol.
5. Detect Coherent Secondary Flow Structures
Note: Use the image acquisition and post processing software and a set of command-line functions (MATLAB-based toolbox, PIVMat 3.01) to import, post-process and analyze 2- component vector fields from the PIV system 5, 6, 33.
Results presented in Figure 7A-D were generated after post processing secondary flow velocity data (see Figures 5, 6) acquired from 2C-2D PIV system shown in Figure 3A. The inflow condition supplied to the curved artery test section with an idealized "Type IV" stent fracture was the carotid artery waveform shown in Figure 4B. Our previous studies have demonstrated the sensitivity of secondary flow structures to decelerating conditions in a variety of pulsatile inflow waveforms within a certain range of Womersley numbers 4 – 6. Accordingly, the time instances of results presented in Figure 7A-D, were chosen to correspond to the systolic deceleration phase of the carotid artery inflow waveform. Coherent secondary flow structures of varying size-strength-morphological characteristics are presented at various planar cross-sections as shown in Figure 7A-D. Large-scale coherent secondary flow structures that emerged in the curved artery test section have been classified as deformed Dean-, Lyne- and Wall-type (D-L-W) vortices. Typically, D-L-W vortices evolve during the systolic acceleration phase. During systolic deceleration phase, D-L-W structures experience an atypical loss in coherence, asymmetry and, changes in vortical positions, sizes, strengths and morphologies. The following is a description of results presented in Figure 7A-D:
At location (Figure 7A): A single pair of symmetric, coherent, deformed Dean vortices (D) are observed in the fields t/T= 0.23 and 0.27. These D-type vortices appear to translate toward the outer wall during deceleration. Q-fields present strain and shear-dominated patterns at t/T= 0.23, in addition to D-type vortices. As a possible effect of deceleration at t/T= 0.27, a reduction in strength of D-type vortices and near-wall strain dominated flow fields is observed. Multi-scale secondary flow morphologies in addition to the D-type are detected in the wavelet transformed vorticity fields indicating presence of several strain-dominated vortical patterns.
At location (Figure 7B): A transition from one pair of D-vortices at the location, to three pairs of symmetric, coherent, deformed Dean-, Lyne- and Wall-type vortices (D-L-W) are observed in the fields. As evidenced by the magnitudes of swirling strengths, the L- and W- type vortices have higher circulation than D-type vortices. Flow perturbations that emanate from the fractured-stent at the location have likely contributed to the formation of D-L-W vortices. The effect of deceleration is observed as the reduction in strength in L- and W-type vortices. There is good agreement in the location of large-scale coherent D-L-W structures between the and fields. Additional smaller scale secondary flow morphologies are detected in the field. Q-fields mostly comprise of shear-dominated near-wall regions allude to the presence of symmetric large-scale D-L-W vortices.
At location (Figure 7C): field at t/T= 0.23 indicates the loss of L-type vortices and the presence of elongated D- and W-type vortices. At t/T= 0.27 there is a loss of swirling strength in both D- and W-type vortices. The effect of deceleration is indicated by the asymmetry of vortical structures observed in the field. Along with the presence of elongated D-type vortices a multitude of small scale W-type vortices are observed. Q-fields indicate the presence of near-wall shear-dominated regions that are suggestive of increased unsteadiness from stent fracture-induced flow perturbations.
At location (Figure 7D): field at t/T= 0.23 comprises of weak, D-L-W structures. Due to the effect of flow deceleration these D-L-W structures tend to deplete further at t/T = 0.27. Loss in near-wall shearing flow is observed in the Q-fields at both instances of time. At t/T = 0.23, field shows that D-vortices are located closer to the inner wall along with multi-scale W-type vortices and surrounding strain-dominated structures in agreement with the corresponding field. field clearly shows a loss in coherence of D-L-W structures and asymmetry at both instances of time while fields do not capture that phenomenon.
Broad inferences after successful execution of the protocol The detected large-scale secondary flow structures and their changing flow morphologies. detected areas of high strain-rate that are normally encountered in the near-wall regions. Continuous wavelet transform algorithm detected the large-scale secondary flow structures in good agreement with unthresholded . The 2D Ricker wavelet kernel additionally resolved several low-circulation, multi-scale secondary flow morphologies that were undetected with and unthresholded . A combination of these three metrics holistically identified secondary flow vortical and strain-dominated structures.
Figure 1. Design, fabrication and installation of straight and curved stents. (A) CAD model of straight stent configuration using a combination of left and right turning helices. (B) CAD model of curved stent configuration. (C) 3D printer used for the fabrication of the stents. (D) and (E) Straight and curved stents after 3D printing. (F) Stents installed in the 180° curved artery test section. Please click here to view a larger version of this figure.
Figure 2. Schematic drawing of the particle image velocimetry (PIV) system. The following system-components are indicated: 1. Nd-YAG laser with optics to produce a laser sheet and 2. CCD camera that is controlled by the PIV-data acquisition computer, 3. Pump instrument control computer that provides the voltage-time waveform to the pump and synchronization triggers to the PIV-data acquisition computer, 4. Programmable gear pump that produces physiological flow rates, 5. A closed loop, experimental test section having inlet and outlet pipes, 180° curved artery test section and a reservoir for the blood-analog fluid. Inset: Various planar cross-sections where PIV measurements can be made. Please click here to view a larger version of this figure.
Figure 3. Experimental arrangement of PIV-system with location of stent installation. (A) Arrangement of the PIV system on the optical table with various system-components. (B) Schematic drawing of the 180° curved artery test section with the important dimensions, location of the straight and curved stent that embody of the 'Type IV' stent fracture and spacing between the fractured stent parts (dspace). Please click here to view a larger version of this figure.
Figure 4. Physiological waveform produced by the programmable pump possessing the characteristics features such as systolic peak at time t/T = 0.19. (A) Flow rate (ml/sec) measured upstream of the 180° curved artery test section over 20 waveform cycles. (B) Flow rate waveform with the standard deviations at various instances of time measured over 20 waveform cycles. Please click here to view a larger version of this figure.
Figure 5. Sequence of PIV measurements and detection of secondary flow structures in the 180° curved artery test section. (A) Generation of secondary flow velocity field data using the PIV technique via synchronization of trigger produced by the pump instrument control computer. (B) Post-processing sequence using the secondary flow field data by treating the pixelated images (or matrices) for Q- and λci – criteria, and wavelet-transformed vorticity (Ω'). Please click here to view a larger version of this figure.
Figure 6. Algorithmic representation of continuous wavelet transform approach for arterial secondary flow structure detection. Insets: 2D-Ricker wavelet at an arbitrary scale (ℓ), an example of 2D vorticity field, Shannon entropy variation with wavelet scale (ℓ). Please click here to view a larger version of this figure.
Figure 7. Secondary flow structures in the 180° curved artery test section at the 45°, 90°, 135° and 180° planar locations and time instances, t/T= 0.23, 0.27, during systolic deceleration. Insets: Schematic drawing depicting the measurement locations, comparison of Q- and λci – criteria, and wavelet-transformed vorticity (Ω') data fields at each planar locations and instances during systolic deceleration, colorbars indicating the range of values of acquired by of Q- and λci – criteria, and wavelet-transformed vorticity (Ω') data and their interpretation. Please click here to view a larger version of this figure.
Parameter | Equation no. | Initialized value | Stent model category | Description | ||
θ = 360 nturns t | 1 | nturns = 4 | Straight; Curved | Number of turns in the helix (nturns) |
||
2 | pitch = 22.225 mm per turn | Straight; Curved | Pitch of the helix (pitch) |
|||
3 | D = 11.84 mm | Straight | Nominal diameter of the stent (D) |
|||
4 | D = 11.84 mm | Straight | Nominal diameter of the stent (D) |
|||
5 | pitch = 22.225 mm per turn | Straight | Pitch of the helix (pitch) Length of straight stent model (z) |
|||
6 | Curved | Radius of the 180° curved artery model Rarc |
||||
β = 180 t | 7 | β = 45 | Curved | Angle subtended by the curved stent at the center of curvature | ||
dwire | – | dwire = 0.85 mm | Straight; Curved | Diameter of stent strut | ||
Lstraight = z | – | Lstraight = 88.9 mm | Straight | Length of straight stent model |
Table 1. Parametric equations of left and right helices and, initialized parameter values.
Chemical Solvents and Reagents | Chemical formula | Density at 20 °C (g/cm3) |
Refractive Index | Kinematic Viscosity (m2/sec) x 10-6 |
Form | CAS Registry number |
Sodium Iodide | NaI | 3.67 | 1.7745 | – | Crystalline | 7681-82-5 |
Glycerol | C3H8O3 | 1.262 | 1.4746 | ≈1115a | Liquid | 56-81-5 |
Deionized Water | H2O | 1 | 1.333 | 1.002 | Liquid | – |
Sodium thiosulfate anhydrous | Na2O3S2 | 1.01 | – | – | Powder | 7772-98-7 |
a Measurements reported by Segur and Oberstar 16 |
Table 2. Description of chemical solvents and reagents used in creating the blood-analog solution.
Parameter | Equation no. | Description | Suggested laboratory equipment |
8 | Density of saturated sodium iodide solution (NaI) is calculated by measuring the mass of the solution and the volume added in small quantities to a 50 ml beaker. | 1. Beaker (50 ml) 2. Weigh Scale 3. Graduated or volumetric pipette |
|
9 | Volume of the entire batch of saturated sodium iodide solution prepared | 1. Beaker with saturated NaI solution (2,000 ml) 2. Weigh scale |
|
10 | Total volume of the blood analog solution expected after volumetric solution preparation | 1. Beaker with saturated NaI solution (2,000 ml) to mix with glycerol and DI water. 2. Weigh scale |
|
11 | Total volume of glycerol to be added to saturated sodium iodide solution | 1. Beaker with saturated NaI solution (2,000 ml) 2. Weigh scale 3. Beaker (100 ml) to transfer glycerol to saturated NaI solution |
|
12 | Total volume of DI water to be added to the solution of saturated NaI and glycerol | 1. Graduated or volumetric pipette to transfer DI water to solution of saturated NaI and glycerol |
Table 3. Table of percentage-by-volume calculations for blood-analog solution: 79% NaI, 20% Glycerol and 1% DI Water.
PIV system specification | Geometry or characteristic value | Description |
Flow geometry | Circular cross-section parallel to light sheet | Curved artery test section |
Maximum in-plane velocity | 0.16 m sec-1 | Secondary flow velocity scale |
Image size | x 1,376 pixels y 1,040 pixels | PIV Camera CCD Array size |
Time interval between laser pulses (∆t) | 600 – 3,200 µsec | Input into PIV Image acquisition software (DaVis 7.2) |
Final number of vectors | x 86, y 65 | Output from PIV data post processing (DaVis 7.2) |
Table 4. Specifications of the two-component, two-dimensional (2C-2D) PIV System.
Supplemental Code File 1. Creating a Mask. Please click here to download this file.
Supplemental Code File 2. Creating a post processing routine. Please click here to download this file.
Supplemental Code File 3: MATLAB codes. Please click here to download this file.
The protocol presented in this paper describes the acquisition of high fidelity experimental data using particle image velocimetry technique (PIV) and coherent structure detection methods, viz., continuous wavelet transforms, , suited for identification of vortex and shear-dominated flows. Analysis of experimental data from physiological inflows in the presence of an idealized "Type IV" fracture reveals that secondary flow structures with complicated hydrodynamic effects such as flow structure asymmetry and variation in spatio-temporal distributions that cannot be predicted from simple fluid dynamic theories.
There are four critical steps in the execution of this protocol viz., (i) Design and fabrication of laboratory-scale stent models, (ii) Preparation of a blood analog working fluid matched with kinematic viscosity of blood and refractive index of the curved artery model, (iii) Non-invasive experimental arrangement (2C-2D PIV) and (iv) Advanced coherent structure detection methods for identifying arterial blood flow patterns.
Womersley number is a dimensionless parameter that relates pulsatile flow frequency to viscous effects 7. Reynolds number relates inertial forces to viscous flow forces. Dean number relates centripetal forces that arise in flow through curved pipes to inertial and viscous forces 1, 2. Details pertaining to the scaling of the physiological waveform with Womersley and Reynolds numbers are presented in 5, 6. The inflow waveform used in this study was reconstructed from archetypal carotid artery flow rate measurements (averaged) of 17-20 healthy patients by Holdworth et al. 15. The pipes leading to the curved artery test section are long enough to allow the flow to be fully developed such that the pulsatile flow conditions at the inlet to the curved artery test section are in-phase with the pump (Figures 3a, 3b and 4a). Repeatability of supplied physiological waveform was ensured by making axial PIV measurements of flow rate and bulk velocity upstream to the model artery using a 2C-2D PIV system (see Fig. 4b).
The hydrodynamic stimuli from arterial hemodynamics toward the aforementioned clinical complications are not well known. Physiological flows involving stent and stent-fractures pose complexities for in vivo and in vitro measurements. The protocol presented herein can be modified to include compliance in the pipes to study the influence of arterial secondary flow structures under non-ideal and more realistic flow scenarios. Such experiments will pose additional challenges in measurement and post-processing of data. The use of stereo- or tomographic-PIV techniques, capable of mapping three-dimensional velocity fields can significantly improve our understanding of the dynamics of secondary flow structures.
The limitations of the experimental arrangement lie in the lack resolution in the near-wall (model-artery lumen) regions and the lack of optical access to the blood flow within the stent implanted regions. These limitations however, pose elegant extensions of the protocol presented. The use of optically clear material for 3D-printing of stents, realistic and patient-specific arterial geometries would allow unprecedented access to hemodynamics of stent-implants and fractured-stents.
An extended outcome of the protocol presented herein relates to the selection of the "best" wavelet scale for coherent structure detection. Steps 5.6.3 – 5.6.7 are a suggested solution to the problem of the "best" wavelet scale (or basis function) in coherent structure detection. The authors found that following steps 5.6.3 – 5.6.7 resolved all the large-scale coherent structures and in addition, detected smaller scale coherent structures that were hitherto undetected in curved artery model experiments. The authors suggest Ref. 34, 35 wherein Shannon entropy is used to evaluate the "best" basis in a discrete wavelet packet transform (DWPT) algorithm toward detecting coherent structures in a turbulent flow experiment. For further information on the approach pertaining to a continuous wavelet transform algorithm, the authors suggest Ref. 5, 6, 35 and the references cited therein.
The incidence of fractures in stent implants and concomitant flow perturbations result in secondary flow structures with complex, multi-scale morphologies and varying size-strength characteristics. Significance of the methodologies such as particle image velocimetry (PIV) in combination with coherent structure detection especially, wavelet transforms allows for resolution of multi-scale, multi-strength secondary flow structures under stent and stent-fracture-induced flow scenarios. The protocol presented herein paves the way for investigating medical complications such as in-stent restenosis (ISR), stent thrombosis and aneurysm formation 8, 11 – 14 due to secondary flows. In addition, secondary flow vortical patterns encountered in the core regions will tend to affect the motion and exposure time of blood-borne particles such as platelets, sensitizing them for activation toward thrombosis. Strain-dominated near-wall (lumen) secondary flow structures will ultimately influence wall shear stress that is closely related to atherogenesis, especially in arterial curvatures.
The analytical procedures to predict secondary flow (vortical) structures are complicated, requiring Navier-Stokes equations in toroidal coordinates and asymptotic theories 1 – 3, 7. A combination of experiments and higher order analytical methods will promote new insights into the hemodynamics of curved arteries prone to several cardiovascular diseases and clinical complication associated with stent implantations and stent fractures.
The authors have nothing to disclose.
The authors acknowledge support from NSF grant CBET-0909678 and funding from the GW Center for Biomimetics and Bioinspired Engineering (COBRE). We thank the students, Mr. Christopher Popma, Ms. Leanne Penna, Ms. Shannon Callahan, Mr. Shadman Hussain, Mr. Mohammed R. Najjari, and Ms. Jessica Hinke for help in the laboratory and Mr. Mathieu Barraja for assisting in CAD drawings.
Acrylic tubes and sheet | McMaster-Carr Supply Company | Inlet and outlet pipes and, material of the curved artery test section | |
Object24 Desktop 3D printer | Stratasys | Desktop rapid prototyping machine. http://www.stratasys.com/. | |
VeroWhitePlus Opaque material | Stratasys | Building material for Object24 Desktop 3D printer | |
Fullcure 705 | Stratasys | Non-toxic gel-like photopolymer Support material for Object24 Desktop 3D printer | |
Ubbelhode viscometer | Cole Parmer | YO-98934-12 | Toward measurement of kinematic viscosity of the blood-analog fluid |
VELP scientifica – ESP stirrer | VELP Scientifica | F206A0179 | Magnetic stirrer |
Ohaus Scout Pro SP 601 | The Lab Depot | SP4001 | Weigh scale |
Refractometer | Atago | PAL-RI | Toward measurement of refractive index of blood-analog fluid |
Beakers, pipettes, syringes and spatula | Sigma-Aldrich | CLS710110, CLS10031L, CLS71015, CLS71011 Z193216 | Toward handling materials required for blood-analog solution preparation |
Sodium Iodide | Sigma-Aldrich | 383112-2.5KG | Crystalline |
Glycerol | Sigma-Aldrich | G5516-1L | Liquid |
Deionized Water | – | – | Liquid |
Sodium thiosulfate anhydrous | Sigma-Aldrich | 72049-250G | Powder |
PIV Recording medium | LaVision | Imager Intense 10Hz | PIV Image acquisition CCD camera |
PIV Illumination source | New Wave Research | Solo III-15 | PIV Laser source, Nd:YAG laser, 532 nm, dual pulse 70 mJ/pulse |
PIV Imaging software | LaVision | DaVis 7.2 | PIV data acquisition and instrument control |
PIV Seeding material | Thermo-scientific | Flouro-Max | Red fluorescent polymer microspheres (≈ 7 µm); Dry dyed polystyrene (DVB) fluorescent microspheres emit bright and distinct colors when illuminated by the light of shorter wavelengths than the emission wavelength. |