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.

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 &#34;Type IV&#34; 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 inclu...

<ol>
	<li>Dean, W. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Note+on+the+motion+of+a+fluid+in+a+curved+pipe.">Note on the motion of a fluid in a curved pipe.</a> Phil Mag. 7, 208-223 (1927).</li><li>Dean, W. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+streamline+motion+of+a+fluid+in+a+curved+pipe.">The streamline motion of a fluid in a curved pipe.</a> Phil Mag. 7, 673-695 (1928).</li><li>Lyne, W. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Unsteady+viscous+flow+in+a+curved+pipe.">Unsteady viscous flow in a curved pipe.</a> J. Fluid. Mech. 45, 13-31 (1970).</li><li>Glenn, A. L., Bulusu, K. V., Shu, F., Plesniak, M. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Secondary+flow+structures+under+stent-induced+perturbations+for+cardiovascular+flow+in+a+curved+artery+model.">Secondary flow structures under stent-induced perturbations for cardiovascular flow in a curved artery model.</a> Int. J. Heat Fluid Fl. 35, 76-83 (2012).</li><li>Bulusu, K. V., Plesniak, M. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Secondary+flow+morphologies+due+to+model+stent-induced+perturbations+in+a+180&#176;+curved+tube+during+systolic+deceleration.">Secondary flow morphologies due to model stent-induced perturbations in a 180&#176; curved tube during systolic deceleration.</a> Exp. Fluids. 54, 1493 (2013).</li><li>Bulusu, K. V., Hussain, S., Plesniak, M. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Determination+of+secondary+flow+morphologies+by+wavelet+analysis+in+a+curved+artery+model+with+physiological+inflow.">Determination of secondary flow morphologies by wavelet analysis in a curved artery model with physiological inflow.</a> Exp. Fluids. 55, 1832 (2014).</li><li>Womersley, J. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Method+for+the+calculation+of+velocity,+rate+of+flow+and+viscous+drag+in+arteries+when+the+pressure+gradient+is+known.">Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known.</a> J. Physiol. 127, 553-563 (1955).</li><li>Sheriff, J., Bluestein, D., Girdhar, G., Jesty, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=High-shear+stress+sensitizes+platelets+to+subsequent+low-shear+conditions.">High-shear stress sensitizes platelets to subsequent low-shear conditions.</a> Ann. Biomed. Eng. 38 (4), 1442-1450 (2010).</li><li>Popma, J. J., Tiroch, K., Almonacid, A., Cohen, S., Kandzari, D. E., Leon, M. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+qualitative+and+quantitative+angiographic+analysis+of+stent+fracture+late+following+sirolimus-eluting+stent+implantation.">A qualitative and quantitative angiographic analysis of stent fracture late following sirolimus-eluting stent implantation.</a> Am. J. Cardiol. 103 (7), 923-929 (2009).</li><li>Kim, S. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+fractured+sirolimus-eluting+stent+with+a+coronary+aneurysm.">A fractured sirolimus-eluting stent with a coronary aneurysm.</a> Ann. Thorac. Surg. 88, 664-665 (2009).</li><li>Adlakha, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Stent+fracture+in+the+coronary+and+peripheral+arteries.">Stent fracture in the coronary and peripheral arteries.</a> J. Interv. Cardiol. 23 (4), 411-419 (2010).</li><li>Alexopoulos, D., Xanthopoulou, I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Coronary+stent+fracture:+How+frequent+it+is?+Does+it+matter.">Coronary stent fracture: How frequent it is? Does it matter.</a> Hellenic J. Cardiol. 52, 1-5 (2011).</li><li>Nair, R. N., Quadros, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Coronary+stent+fracture:+A+review+of+the+literature.">Coronary stent fracture: A review of the literature.</a> Cardiac. Cath. Lab Director. 1, 32-38 (2011).</li><li>Jaff, M., Dake, M., Popma, J., Ansel, G., Yoder, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Standardized+evaluation+and+reporting+of+stent+fractures+in+clinical+trials+of+noncoronary+devices.">Standardized evaluation and reporting of stent fractures in clinical trials of noncoronary devices.</a> Catheter Cardiovasc. Interv. 70, 460-462 (2007).</li><li>Holdsworth, D., Norley, C. J., Frayne, R., Steinman, D. A., Rutt, B. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Characterization+of+common+carotid+artery+blood-flow+waveforms+in+normal+human+subjects.">Characterization of common carotid artery blood-flow waveforms in normal human subjects.</a> Physiol. Meas. 20 (3), 219-240 (1999).</li><li>Deutsch, S., Tarbell, J. M., Manning, K. B., Rosenberg, G., Fontaine, A. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Experimental+fluid+mechanics+of+pulsatile+artificial+blood+pumps.">Experimental fluid mechanics of pulsatile artificial blood pumps.</a> Annu. Rev. Fluid Mech. 38, 65-86 (2006).</li><li>Yousif, M. Y., Holdsworth, D. W., Poepping, T. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+blood-mimicking+fluid+for+particle+image+velocimetry+with+silicone+vascular+models.">A blood-mimicking fluid for particle image velocimetry with silicone vascular models.</a> Exp. Fluids. 50, 769-774 (2011).</li><li>Budwig, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Refractive+index+matching+methods+for+liquid+flow+investigations.">Refractive index matching methods for liquid flow investigations.</a> Exp. Fluids. 17, 350-355 (1994).</li><li>Hunt, J. C. R., Wray, A. A., Moin, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Eddies,+stream,+and+convergence+zones+in+turbulent+flows.">Eddies, stream, and convergence zones in turbulent flows.</a> Center for Turbulence Research. , (1988).</li><li>Adrian, R. J., Christensen, K. T., Liu, Z. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Analysis+and+interpretation+of+instantaneous+turbulent+velocity+fields.">Analysis and interpretation of instantaneous turbulent velocity fields.</a> Exp. Fluids. 29, 275-290 (2000).</li><li>Chong, M., Perry, A. E., Cantwell, B. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+general+classification+of+three-dimensional+flow+fields.">A general classification of three-dimensional flow fields.</a> Phys. Fluids A. 2 (5), 765-777 (1990).</li><li>Zhou, J., Adrian, R. J., Balachandar, S., Kendall, T. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Mechanisms+for+generating+coherent+packets+of+hairpin+vortices+in+channel+flow.">Mechanisms for generating coherent packets of hairpin vortices in channel flow.</a> J. Fluid Mech. 387, 353-396 (1999).</li><li>Haller, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=An+objective+definition+of+a+vortex.">An objective definition of a vortex.</a> J. Fluid Mech. 525, 1-26 (2005).</li><li>Chakraborty, P., Balachander, S., Adrian, R. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=On+the+relationships+between+local+vortex+identification+schemes.">On the relationships between local vortex identification schemes.</a> J. Fluid Mech. 535, 189-214 (2005).</li><li>Wallace, J. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Twenty+years+of+experimental+and+direct+numerical+simulation+access+to+the+velocity+gradient+tensor:+What+have+we+learned+about+turbulence.">Twenty years of experimental and direct numerical simulation access to the velocity gradient tensor: What have we learned about turbulence.</a> Phys. Fluids. 21, 021301 (2009).</li><li>Farge, M., Guezennec, Y., Ho, C. M., Meneveau, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Continuous+wavelet+analysis+of+coherent+structures.">Continuous wavelet analysis of coherent structures.</a> Center for Turbulence Research, Proceedings of the Summer Program. , 331-348 (1990).</li><li>Himburg, H. A., Friedman, M. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Correspondence+of+Low+Mean+Shear+and+High+Harmonic+Content+in+the+Porcine+Iliac+Arteries.">Correspondence of Low Mean Shear and High Harmonic Content in the Porcine Iliac Arteries.</a> ASME J. Biomedical Eng. 128, 852-856 (2006).</li><li>Dai, G., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Distinct+endothelial+phenotypes+evoked+by+arterial+waveforms+derived+from+atherosclerosis-susceptible+and+-resistant+regions+of+human+vasculature.">Distinct endothelial phenotypes evoked by arterial waveforms derived from atherosclerosis-susceptible and -resistant regions of human vasculature.</a> PNAS. 101 (41), 14871-14876 (2004).</li><li>Hanus, J., Zahora, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Measurement+and+comparison+of+mechanical+properties+of+nitinol+stents.">Measurement and comparison of mechanical properties of nitinol stents.</a> Physica Scripta. 118, 264-267 (2005).</li><li>Segur, J. B., Oberstar, H. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Viscosity+of+glycerol+and+its+aqueous+solutions.">Viscosity of glycerol and its aqueous solutions.</a> Ind. Eng. Chem. 43, 2117-2120 (1951).</li><li>Adrian, R. J., Westerweel, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Particle image velocimetry. , (2011).</li><li>Raffel, M., Willert, C. E., Wereley, S. T., Kompenhans, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Particle image velocimetry - A practical guide, 2nd ed. , (2007).</li><li>Moisy, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> PIVmat 3.01 software. , (2013).</li><li>Ruppert-Felsot, J. E., Praud, O., Sharon, E., Swinney, H. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Extraction+of+coherent+structures+in+a+rotating+turbulent+flow+experiment.">Extraction of coherent structures in a rotating turbulent flow experiment.</a> Physical Review E. 72, 016311 (2005).</li><li>Bulusu, K. V., Plesniak, M. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Shannon+entropy-based+wavelet+transform+methods+for+autonomous+coherent+structure+identification+in+fluid+flow+field+data.">Shannon entropy-based wavelet transform methods for autonomous coherent structure identification in fluid flow field data.</a> Entropy. 17 (10), 6617-6642 (2015).</li></ol>

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, <img alt="Equation 1" src="/files/ftp_upload/51288/51288eq1.jpg" />, suited for identification of vortex and shear-dominated flows. Analysis of experimental data from physiological inflows in the presence of an idealized &#34;Type IV&#34; fracture reveals that secondary flow structures with complica...

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 &#34;Types I-to-IV&#34;, wherein &#34;Type IV&#34; 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 &#34;Type IV&#34; 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 &#34;Type IV&#34; 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 (&#916;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, <img alt="Equation 1" src="/files/ftp_upload/51288/51288eq1.jpg" /> 5, 6, 19 - 24, 26.
The authors note that the velocity gradient tensor is essentially, a 3 x 3 matrix, 
<img alt="Equation 2" src="/files/ftp_upload/51288/51288eq2.jpg" />.
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 <img alt="Equation 3" src="/files/ftp_upload/51288/51288eq3.jpg" /> of the PIV image <img alt="Equation 4" src="/files/ftp_upload/51288/51288eq4.jpg" /> should be a 2 x 2 matrix, <img alt="Equation 5" src="/files/ftp_upload/51288/51288eq5.jpg" />. The z-component vorticity <img alt="Equation 6" src="/files/ftp_upload/51288/51288eq6.jpg" /> for each pixel <img alt="Equation 7" src="/files/ftp_upload/51288/51288eq7.jpg" /> is computed using the anti-symmetric part of the velocity gradient tensor <img alt="Equation 8" src="/files/ftp_upload/51288/51288eq8.jpg" />. The result will be a 2D array of vorticity <img alt="Equation 9" src="/files/ftp_upload/51288/51288eq9.jpg" /> 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., <img alt="Equation 10" src="/files/ftp_upload/51288/51288eq10.jpg" /> and wavelet transformed vorticity <img alt="Equation 11" src="/files/ftp_upload/51288/51288eq11.jpg" /> are applied to velocity field data toward detection of multi-scale, multi-strength occurrences of secondary flow structures downstream of the idealized &#34;Type IV&#34; stent fracture.
The <img alt="Equation 12" src="/files/ftp_upload/51288/51288eq12.jpg" />, 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; <img alt="Equation 13" src="/files/ftp_upload/51288/51288eq13.jpg" />. Norm of the strain rate is then calculated; <img alt="Equation 14" src="/files/ftp_upload/51288/51288eq14.jpg" />. Vorticity is computed from the anti- symmetric part. Enstrophy or square of z-component vorticity, <img alt="Equation 15" src="/files/ftp_upload/51288/51288eq15.jpg" />) is then computed. The <img alt="Equation 16" src="/files/ftp_upload/51288/51288eq16.jpg" /> is finally computed; <img alt="Equation 17" src="/files/ftp_upload/51288/51288eq17.jpg" />. A contour plot of the entire set of <img alt="Equation 18" src="/files/ftp_upload/51288/51288eq18.jpg" /> with iso-regions of <img alt="Equation 19" src="/files/ftp_upload/51288/51288eq19.jpg" />, will indicate secondary flow structures 19. 
The <img alt="Equation 20" src="/files/ftp_upload/51288/51288eq20.jpg" />, also known as &#39;swirling strength&#39; 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&#160;<img alt="Equation 21" src="/files/ftp_upload/51288/51288eq21.jpg" /> are computed. The eigenvalues should be of the form, <img alt="Equation 22" src="/files/ftp_upload/51288/51288eq22.jpg" />. A contour plot of <img alt="Equation 23" src="/files/ftp_upload/51288/51288eq23.jpg" /> with iso-regions of <img alt="Equation 24" src="/files/ftp_upload/51288/51288eq24.jpg" /> 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 <img alt="Equation 25" src="/files/ftp_upload/51288/51288eq25.jpg" />5, 6, 26. By convolving a dilated or contracted wavelet with a 2D vorticity field, wavelet transformed vorticity&#160;<img alt="Equation 26" src="/files/ftp_upload/51288/51288eq26.jpg" /> 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 <img alt="Equation 27" src="/files/ftp_upload/51288/51288eq27.jpg" /> for each pixel <img alt="Equation 21" src="/files/ftp_upload/51288/51288eq21.jpg" /> such that&#160; <img alt="Equation 28" src="/files/ftp_upload/51288/51288eq28.jpg" />, 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.

<table border="0" cellpadding="0" cellspacing="0" width="698">
	<tbody>
		<tr height="63">
			<td height="63" style="height:63px;width:216px;">Acrylic tubes and sheet</td>
			<td style="width:155px;">McMaster-Carr Supply Company</td>
			<td style="width:119px;"></td>
			<td style="width:208px;">Inlet and outlet pipes and, material of the curved artery test section</td>
		</tr>
		<tr height="63">
			<td height="63" style="height:63px;width:216px;">Object24 Desktop 3D printer</td>
			<td style="width:155px;">Stratasys</td>
			<td style="width:119px;"></td>
			<td style="width:208px;">Desktop rapid prototyping machine. http://www.stratasys.com/.&#160;</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:216px;">VeroWhitePlus Opaque material</td>
			<td style="width:155px;">Stratasys</td>
			<td style="width:119px;"></td>
			<td style="width:208px;">Building material for Object24 Desktop 3D printer</td>
		</tr>
		<tr height="84">
			<td height="84" style="height:84px;width:216px;">Fullcure 705</td>
			<td style="width:155px;">Stratasys</td>
			<td style="width:119px;"></td>
			<td style="width:208px;">&#160;Non-toxic gel-like photopolymer Support material for Object24 Desktop 3D printer</td>
		</tr>
		<tr height="63">
			<td height="63" style="height:63px;width:216px;">Ubbelhode viscometer</td>
			<td style="width:155px;">Cole Parmer</td>
			<td>YO-98934-12</td>
			<td style="width:208px;">Toward measurement of kinematic viscosity of the blood-analog fluid</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:216px;">VELP scientifica - ESP stirrer&#160;</td>
			<td style="width:155px;">VELP Scientifica</td>
			<td>F206A0179</td>
			<td style="width:208px;">Magnetic stirrer</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:216px;">Ohaus Scout Pro SP 601&#160;</td>
			<td style="width:155px;">The Lab Depot</td>
			<td>SP4001</td>
			<td style="width:208px;">Weigh scale</td>
		</tr>
		<tr height="63">
			<td height="63" style="height:63px;width:216px;">Refractometer</td>
			<td>Atago</td>
			<td style="width:119px;">PAL-RI</td>
			<td style="width:208px;">Toward measurement of refractive index of blood-analog fluid</td>
		</tr>
		<tr height="105">
			<td height="105" style="height:105px;width:216px;">Beakers, pipettes, syringes and spatula</td>
			<td style="width:155px;">Sigma-Aldrich&#160;</td>
			<td style="width:119px;">CLS710110,&#160; CLS10031L, CLS71015, CLS71011 Z193216</td>
			<td style="width:208px;">Toward handling materials required for blood-analog solution preparation</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:216px;">Sodium Iodide</td>
			<td style="width:155px;">Sigma-Aldrich</td>
			<td>383112-2.5KG&#160;</td>
			<td style="width:208px;">Crystalline</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:216px;">Glycerol</td>
			<td style="width:155px;">Sigma-Aldrich</td>
			<td>G5516-1L</td>
			<td style="width:208px;">Liquid</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:216px;">Deionized Water</td>
			<td style="width:155px;">-</td>
			<td style="width:119px;">-</td>
			<td style="width:208px;">Liquid</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:216px;">Sodium thiosulfate anhydrous</td>
			<td style="width:155px;">Sigma-Aldrich</td>
			<td>72049-250G</td>
			<td style="width:208px;">Powder</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:216px;">PIV Recording medium</td>
			<td style="width:155px;">&#160;LaVision&#160;</td>
			<td style="width:119px;">Imager Intense 10Hz</td>
			<td style="width:208px;">PIV Image acquisition CCD camera</td>
		</tr>
		<tr height="63">
			<td height="63" style="height:63px;width:216px;">PIV Illumination source</td>
			<td style="width:155px;">New Wave Research</td>
			<td style="width:119px;">Solo III-15</td>
			<td style="width:208px;">PIV Laser source, Nd:YAG laser, 532 nm, dual pulse 70 mJ/pulse</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:216px;">PIV Imaging software</td>
			<td style="width:155px;">&#160;LaVision&#160;</td>
			<td style="width:119px;">DaVis 7.2&#160;</td>
			<td style="width:208px;">PIV data acquisition and instrument control&#160;</td>
		</tr>
		<tr height="211">
			<td height="211" style="height:211px;width:216px;">PIV Seeding material</td>
			<td style="width:155px;">Thermo-scientific&#160;&#160;</td>
			<td style="width:119px;">Flouro-Max</td>
			<td style="width:208px;">Red fluorescent polymer microspheres (&#8776; 7 &#181;m); Dry dyed polystyrene (DVB) fluorescent microspheres emit bright and distinct colors when illuminated by the light of shorter&#160; wavelengths than the emission wavelength.&#160;</td>
		</tr>
	</tbody>
</table>

experimental investigation of secondary flow structures downstream of a model type iv stent failure in a 180&#176; curved artery test section

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-&#224;-vis those under unperturbed flows 5. These stent failures have been classified as &#34;Types I-to-IV&#34; 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 (&#34;Type IV&#34;) 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.

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 &#34;Type IV&#34; stent fracture was the carotid artery waveform shown in Figure 4B. Our previous studies have demonstrated the sensitivity of secondary flow structures to de...

Watch this Scientific Journal Video about Experimental Investigation of Secondary Flow Structures Downstream of a Model Type IV Stent Failure in a 180Â° Curved Artery Test Section at JoVE.com

Experimental Investigation of Secondary Flow Structures Downstream of a Model Type IV Stent Failure in a 180&amp;#176; Curved Artery Test Section

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.

Experimental Investigation of Secondary Flow Structures Downstream of a Model Type IV Stent Failure in a 180&#176; Curved Artery Test Section

The arterial network in the human vasculature comprises of ubiquitously present blood vessels with complex geometries (branches, curvatures and ...

This video shows how to use Particle Image Velocimetry, and Coherent Structure Detection, to detect secondary flow or wortical structures in a curved artery model, to better understand the clinical relevance of Type IV stent failures. This method can help answer key questions in cardiovascular hemodynamics. Such as, how stents, which can deteriorate over time and ultimately fracture, can influence blood flow patterns and affect disease progression in vascular remodeling.The main advantage of this technique is that it facilitates investigation of in-vivo and in-vitro measurements in complex physiological arterial flows involving stent and stent fractures. Demonstrating procedure will be Mohammad Najjari, a graduate student, and Jessica Hinke, an undergraduate student. Prepare a blood-analog solution comprised of 79%saturated sodium iodide solution, 20%glycerol, and 1%deionized water.Start with saturated sodium iodide solution, then, add the glycerol in small increments using a syringe until the entire volume of glycerol has been added. Between additions, wait for the solution to become visibly homogenized, and take note of the volume of each glycerol addition. Next, add the required volume of water, continue stirring until the blood-analog solution is visibly homogenized.Finally, add very small quantities of sodium thiosulfate, using a spatula until the solution is optically clear. To begin, create high resolution STL files from the straight and curved stent Cad models. In the software, select Export and Model from the file menu.Choose the STL option, set the Chord Height to zero, and set the Angle Control to one. The angle control value regulates the amount of tessellation along the surface with small radii, and the setting can be adjusted between zero and one. Then, click OK to create the STL file.Next, fabricate the stent models on a Rapid Prototyping machine. To do this, start the 3D printing software and click on Insert to locate the STL file on the 3D printer computer. Then, select the desired file and drive the mouse to place the STL file on a virtual platform on the screen.Next, using the 3D-printed parts, install the stent into the curved artery test section to recreate an Idealized Type IV fracture scenario. Entailing a complete transverse fracture of stents and linear displacement of fragmented parts. Now, assemble the setup on an optical table.Connect the straight acrylic pipes to the inlet and outlet of a 180 degree curved artery test section. To begin, grossly adjust the laser's position by illuminating a small piece of paper. Make the laser sheet approximately 2mm thick by adjusting the focus.Next, target the laser sheet along the 90 degree portion of the measurement region with the sheet perpendicular to the table. Then, place a camera near the zero degree or 180 degree location to view the cross sectional region illuminated by the laser sheet. Incrementally adjust the alignment of the laser and camera so the camera views the circular cross section of the curved artery with minimal particle distortion, which is assessed using the Grab function in the software.Once the camera and laser are in position, acquire images of the secondary flow fields. First, turn on the programmable pump. In the pump instrument control program, set the Amplitude to one volt.Set the DC Offset to zero volts. Set the number of timed steps to 1000, and set the time period to four seconds. Additional settings are discussed in the Text Protocol.Then, load the text file that has the values for the voltage time wave form, and run the program to supply the blood-analog fluid to the experiment. In the PIV Data Acquisition program, after clicking the New Recordings tab, select the device under the Settings section, and confirm that the laser is set to On with the appropriate power settings. Navigate to laser control to confirm the values.Navigate to Timing and set it to External Cyclic Trigger. Click on Acquisit, select Table Scan, Edit Table Scan, Append Scan, and populate scan settings by entering start time representative of systolic deceleration. Time increment and end time of PIV data acquisition.Reference Time, DT1, or the time between laser pulses should be adjusted. Then, click on Close. Under Acquisit, click on image acquisition and enter the number of images as 200.Now, the PIV system is ready to acquire data. Select Start Recording to acquire phase-wise measurements, using the trigger signal from the pump instrument control. The set number of plane or velocity fields at each time instance, and at the prearranged 90 degree test section location will then be acquired.When the recording is done, power down the laser. Turn off the pump, turn off the camera, and replace the lens cover. Using the supplemental code file, compute the Phase Averaged and RMS Secondary Flow Velocities, Vorticity, and Swirling Strength fields.In the software, initialize setting pertaining to defined scale and send to Mask Definition. Review the Supplemental Code file for details. In the software, right click on any PIV data, and initiate Batch Processing.Select the operation Vector Statistics, Vector Field Results from the group Statistics, and click on Parameter in the Dialog box. Toggle the average V and RMSV options under the Vector Field section. Then, select the operation Vorticity from the group Extract Scalar Fields Rotation and Shear.This finds the two dimensional vorticity in the planar cross section. Now, start the processing by making a right click on any PIV data under the Project window, and selecting Hyperloop All Sets. To compute the Swirling Strength fields for the detection of Secondary Flow Structures, right click on any PIV data, select Velocity Vector Field in the tree, and initiate batch processing.Then, select the operation Swirling Strength from the group Extract Scalar Fields Rotation and Shear, and click on Close. Now, right click on the same Velocity Vector Field, select Hyperloop All Sets, and adjust the parameters for Batch Processing as before. Finally, execute the computation of the Swirling Strength fields.Secondary Flow Velocity data at the 90 degree cross-sectional location, was acquired from the 2C 2D PIV system. The inflow condition supplied to the curved artery test section with an Idealized Type IV stent fracture was the carotid artery wave form. The Secondary Flow Velocity field data was generated using the PIV technique, via synchronization of Trigger produced by the pump instrument control computer.Post-processing of the data was applied to the pixelated data, to determine the Q-criterion and Swirling Strength or Landis ci-criterion during the systolic deceleration. The Swirling Strength and Q-criterion-based flow fields were compared. The wortical patterns marked with DLW characters, representative of deformed Dean-Lynn-and Wall-type DLW wortical patterns, and strain-dominated regions, were examined at four distinct cross-sectional locations during systolic deceleration.The continuous Wavelet Transformed-based method was applied to the secondary flow vorticity field, that enabled vortex detection with greater detail in shape, size, and strength. Finally, the Q-criterion, Swirling Strength, and Wavelet Transformed-based flow fields were compared. The Wavelet Transformed-based method revealed the presence of complex wortical patterns with great detail in size, scale, and strengths.Once mastered, this technique can be done in approximately eight hours, including PIV data acquisition and post-processing, after the model stents are printed on the 3D printer. While attempting this procedure, it's important to remember to make sure the blood-analog fluid refractive index is adequately matched using minute quantities of sodium thiosulfate. In future, to improve paratheory result accuracy, and minimize the optical distortion, we will be using a PIV calibration map like this.Following this procedure, flows in patient-specific blood vessels can be evaluated in order to answer additional questions pertaining to pathological conditions, allowing unprecedented access to hemodynamics of prosthetics, stent implants, and fractured stents.

experimental-investigation-secondary-flow-structures-downstream-model

Research

JoVE Journal

Bioengineering

11.5K vistas.  The George Washington University. This video shows how to use Particle Image Velocimetry, and Coherent Structure Detection, to detect secondary flow or wortical structures in a curved artery model, to better understand the clinical relevance of Type IV stent failures. This method can help answer key questions in cardiovascular hemodynamics. Such as, how stents, which can deteriorate over time and ultimately fracture, can influence blood flow patterns and affect disease progression in vascular remodeling.The main advant...