Name: the diffusion of passive tracers in laminar shear flow
Uploaded: 2018-05-01T13:30:00
Duration: 8 min 1 s
Description: Watch this Scientific Journal Video about The Diffusion of Passive Tracers in Laminar Shear Flow at JoVE.com

We acknowledge funding from the Office of Naval Research (grant DURIP N00014-12-1-0749) and the National Science Foundation (grants RTG DMS-0943851, CMG ARC-1025523, DMS-1009750, and DMS-1517879). Additionally, we acknowledge the work of Sarah C. Burnett who helped develop an early version of the experimental setup and protocol.

1 . Prepare the parts to build the experimental setup
<ol>
	<li>Utilize the 3D CAD drawings attached (.stl format) to 3D-print an injector post, a reservoir, a hexagonal connector, and two plates to be used as mounts for the pipes (two for each geometry). 
	NOTE: Alternatively, certain parts of the setup can be laser-cut. In this report, the square thick pipe has been mounted with laser-cut plates, while the rectangular thin pipe has been mounted with 3D-printed plates.</li>
	<li>Obtain smooth glass capillary pipes of the desired geometry. 
	NOTE: In this report, two pipe geometries are used: 30 cm-long pipe of square cross-section - internal cross-section 1 mm x 1 mm and wall thickness 0.2 mm; 30 cm-long pipe of rectangular cross-section - internal cross-section 1 mm x 10 mm and wall thickness 0.7 mm. The square pipe is henceforth referred to as the thick pipe, whereas the rectangular pipe is referred to as the thin pipe.</li>
</ol>
2 . Assembly of the experimental setup
<ol>
	<li>Tapping of the 3D-Printed parts
	<ol>
		<li>Tap the injector post on both sides with a 1/8&#34; (0.32 cm) NPT tap where the injection needle and dye input will be installed. Tap the reservoir in the back with a 10-32 tap where the draining tube will be installed.</li>
		<li>Tap the four screw holes with a 6-32 tap on the front of the reservoir. Tap the hexagonal connector piece on the top and bottom with a 6-32 tap.</li>
	</ol>
	</li>
	<li>Prepare the tapped 3D-printed parts
	<ol>
		<li>Injector post
		<ol>
			<li>Cover the threads of a barbed hose fitting with PTFE sealing tape. Screw the prepared fitting onto the back hole of the injector post. Cut a 30 cm-long piece of plastic tubing (inner diameter 3.30 mm). Insert the tube on the hose adaptor.</li>
			<li>Cover the threads of the stainless-steel dispensing needle (outer diameter 0.71 mm) with PTFE sealing tape. Screw the stainless-steel dispensing needle on the front (large) hole onto the injector post.</li>
		</ol>
		</li>
		<li>Reservoir
		<ol>
			<li>Cover the threads of a small barbed hose fitting with PTFE sealing tape. Screw the prepared fitting onto the back hole of the reservoir (smaller hole).</li>
			<li>Cut a 30 cm-long piece of plastic tubing (inner diameter 3.30 mm). Insert the tube onto the hose adaptor. Close the other end of the tube with a small cap. 
			NOTE: This will be the draining system for the reservoir.</li>
			<li>Place a rubber O-Ring (Oil-Resistant Buna-N O-Ring, 1/16&#34; (0.16 cm) Fractional Width, Dash Number 016) in the circular recession on the pipe side of the reservoir.</li>
		</ol>
		</li>
		<li>Hexagonal connector
		<ol>
			<li>Cover the threads of a small barbed hose fitting with PTFE sealing tape. Screw the prepared fitting onto the bottom hole of the hexagonal connector.</li>
			<li>Cut a 30 cm-long piece of plastic tubing (inner diameter 3.30 mm). Insert the tube on the hose adaptor.</li>
			<li>Cover a hose adaptor with PTFE sealing tape. Make sure to cover the hose adaptor going against the threads.</li>
			<li>Cut a 4 cm-long piece of plastic tubing (inner diameter 3.30 mm). Insert the tube on the hose adaptor.</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>Prepare the pipe
	<ol>
		<li>Distribute a thin layer of RTV rubber sealant 2 mm away from each end of the pipe. Spread the sealant evenly around the outside of the pipe and make sure not to obstruct the pipe access with the sealant.</li>
		<li>Mount the pipe onto the 3D-printed plates by inserting it carefully into the pre-cut holes on the 3D-printed pipe adapters. Make sure to push the pipe in at least 2 mm so that the sealant along each side contacts with the plates.</li>
		<li>Carefully spread the sealant onto the edge of the plate so that the pipe gets sealed into the cutout. Wait at least 12 h for the sealant to fully vulcanize thus sealing the pipe onto the plates.</li>
	</ol>
	</li>
	<li>Measure 0.40 g of fluorescein powder to prepare the dye solution. Dilute powder into 0.50 L of distilled water to obtain the desired dye concentration (0.80 g/L concentration). 
	NOTE: The diffusivity of fluorescein in water is estimated by performing a least-squares fit of the analytical expression for the second moment of the cross-sectionally averaged tracer distribution in the circular pipe geometry14 to the experimental measurement of the same quantity. The molecular diffusion coefficient is estimated to be &#954; = 5.7 x 10-6 cm2/s, consistent with previously published values of diffusivity of fluorescein in pure water.</li>
	<li>Assembly
	<ol>
		<li>Syringe pump A setup
		<ol>
			<li>Fill a 12 mL plastic syringe with a rubber plunger with distilled water. Insert a plastic dispensing tip onto the syringe. Mount the syringe onto syringe pump A. Connect the syringe to the 30 cm-long tube inserted at the bottom of the hexagonal connector.</li>
			<li>Fill a 1 mL plastic syringe with a rubber plunger with distilled water. Mount the syringe onto syringe pump A. Cut a 30 cm-long piece of plastic tubing (inner diameter 3.30 mm). Attach it to the 1 mL plastic syringe. 
			NOTE: Both syringes filled with distilled water are mounted on syringe pump A. As the pump is activated, water will be ejected from both syringes. The first one to be used is the 12 mL syringe, so the 1 mL syringe needs to be connected to a draining tube to avoid water spills. This step is not necessary for the thin rectangular pipe.</li>
		</ol>
		</li>
		<li>Injector Setup
		<ol>
			<li>Fill a 3 mL plastic syringe with a rubber plunger with the fluorescein solution. Insert a plastic dispensing tip onto the syringe.</li>
			<li>Attach the tube connected to the back of the injector to the dye syringe.</li>
			<li>Fill the injector post with the dye solution by manually injecting dye through the syringe while holding the injector post horizontally (i.e. with the needle oriented upwards and above the syringe). Keep pushing on the syringe until the injector is completely full of dye and no air is trapped inside.</li>
			<li>Mount the syringe onto syringe pump B. Clamp the injector post to the edge of the lab bench in a way that it is reachable by the tube connected to the syringe pump.</li>
			<li>Insert small washers on four long screws (Stainless Pan Head Phillips Machine Screws 6-32 Thread, 2-1/4&#34; (5.76 cm) length). Insert the four screws in the four holes surrounding the needle. 
			NOTE: Make sure the head of the screw is on the back of the injector post (on the same side as the tube connected to the dye syringe).</li>
		</ol>
		</li>
		<li>Hexagonal connector
		<ol>
			<li>Place two O-Rings (Oil-Resistant Buna-N O-Ring, 1/16&#34; (0.16 cm) Fractional Width, Dash Number 016) in the circular cutouts on each side of the hexagonal connector.</li>
			<li>Attach the hexagonal connector to the injector post by aligning its holes to the four screws and inserting it on them. Make sure to have the side with the larger hole facing the injector post. Check and ensure that the O-Ring does not move out of place when clamped between the two parts.</li>
		</ol>
		</li>
		<li>Pipe
		<ol>
			<li>Attach one of the end-plates connected to the pipe to the hexagonal connector by aligning its holes to the four screws and inserting it onto them. Pay close attention to the needle which needs to enter the pipe as it is being mounted.</li>
			<li>Secure the four long screws to compress together the injector, the hexagonal connector, and the pipe-adapter plate by attaching four 6-32 stainless steel nuts to the end of the long bolts. Ensure that the O-Rings do not move out of place when clamped between the parts.</li>
			<li>Attach the opposite end of the pipe to the reservoir by using four short screws and washers (Stainless Pan Head Phillips Machine Screws 6-32 Thread, 1/2&#34; (1.27 cm) length). Check that the O-Ring does not move out of place when compressed between the two parts.</li>
		</ol>
		</li>
		<li>Clamp the reservoir to the table. Make sure the reservoir is aligned with the injector post to not bend the pipe.</li>
		<li>Air extraction system: Insert a plastic dispensing tip into the tube connected to the top of the hexagonal connector. Attach a 3 mL syringe to the plastic tip. 
		NOTE: This syringe will be used to extract any air bubbles trapped in the system.</li>
		<li>Lights and Camera
		<ol>
			<li>Place two 61 cm-long UV-A tube lights on each side of the experimental setup. 
			NOTE: There is a specifically designed track on each side of both the injector and reservoir. The experiment should be run in the dark with the UV-A tube lights turned on.</li>
			<li>Place a camera with memory card above the experimental setup facing down. 
			NOTE: The camera should be positioned at least 1 m above the pipe. In this way, the frame will include the entire pipe length. A DSLR camera was used with a lens of adjustable focal length, 24 - 120 mm.</li>
			<li>Program the camera using a remote trigger to take pictures every 1 s with aperture 5.6f, shutter speed 5, and ISO 200.</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
3 . Experimental run
<ol>
	<li>Setup
	<ol>
		<li>Fill the reservoir with distilled water to a level slightly above the pipe. Fill pipe with distilled water by pushing on the syringe pump. Turn on the UV-A tube lights and pull the blackout curtains.</li>
		<li>Run the programmable syringe pump A to flush the pipe of any residual dye.</li>
		<li>Take a single reference image of the pipe filled with pure distilled water. 
		NOTE: This is the reference&#160;shot that will be used in the data processing steps later. This picture needs to be taken in the dark in conditions as similar as possible to the experimental run.</li>
		<li>Switch the tube connecting to the injector post to the 1 mL syringe mounted on syringe pump A. Connect the 12 mL syringe to the draining tube (previously connected to the 1 mL syringe). 
		NOTE: This step is not necessary for the thin rectangular pipe.</li>
	</ol>
	</li>
	<li>Initial Condition
	<ol>
		<li>Inject a 1 mm-thick dollop of dye (3 mm-thick for the thin rectangular tube) in the pipe by running the analog syringe pump B. 
		NOTE: This step creates the dye initial condition. The amount of dye injected depends on the geometry of the pipe used. The thin tube requires a larger amount of dye because its cross-sectional area is larger. Before the experimental run, the dye will have to diffuse across the cross-section and injecting a larger amount of dye ensures that it will be bright enough to be captured in photographs even after it has diffused.</li>
		<li>Program syringe pump A to inject distilled water at the very slow flow rate of 0.193 mL/h for the thick square pipe (the flow rate is 1.93 mL/h for the thin rectangular pipe). Run the syringe pump for 5 min to allow the bolus of dye to be transported down the pipe away from the needle. 
		NOTE: After 5 min, the dye should be approximately 1 cm away from the needle. The increase in the flow rate by one order of magnitude for the thin pipe is because the volume of the thin pipe is 10 times that of the thick pipe.</li>
		<li>Pull the dye syringe backwards manually, making sure the dye does not reach the needle. 
		NOTE: This will ensure that there is distilled water at the end of the needle so that no more dye will be dispersed into the pipe during the experimental run.</li>
		<li>Wait for a time tw &#62; t*d for the dye bolus to diffuse across the cross-section of the pipe. 
		NOTE: The diffusive time t*d = b2/&#954; considers the characteristic length b to be half the long cross-sectional side. This way of computing the wait time is generalizable to any cross-section with an appropriate choice of b. For our representative results, the wait time was 15 min for the thick square pipe and 15 h for the thin rectangular pipe.</li>
	</ol>
	</li>
	<li>Flow
	<ol>
		<li>Program syringe pump A to the desired flow rate of 1.93 mL/h for the thick square pipe and 19.3 mL/h for the thin rectangular pipe.</li>
		<li>Start the syringe pump and the remote trigger on the camera at the same time. Run the experiment for 5 min, with an interval between pictures of 1 s.</li>
		<li>Turn the room lights on and take an image of a&#160;ruler placed at the same height as the pipe and parallel to it. 
		NOTE: This will help determine the length scale (pixels/mm) used in data processing.</li>
	</ol>
	</li>
</ol>
4. Data Processing
<ol>
	<li>Extract the memory card from the camera and download the data to a computer where image processing software will be used to analyze it.</li>
	<li>MATLAB analysis
	<ol>
		<li>First subtract the reference image shot (snapped in step 3.1.3) from the first experimental image.</li>
		<li>Crop the image along the upper and lower edges of the pipe. Make sure to rotate the image if the pipe is not aligned with the frame.</li>
		<li>Sum the intensity reading of the green channel vertically in the resulting image. 
		NOTE: This is proportional to the total cross-sectional dye intensity as a function of the length along the pipe.</li>
		<li>Convert the units of length from pixels to mm by using the physical length scale from the calibration image (see step 3.3.3).</li>
		<li>Repeat for all remaining images. This results in a time sequence of the curves measuring the total dye concentration along the length of the pipe.</li>
	</ol>
	</li>
</ol>

<ol>
	<li>Dutta, D., Leighton, D. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Dispersion+in+Large+Aspect+Ratio+Microchannels+for+Open-Channel+Liquid+Chromatography.">Dispersion in Large Aspect Ratio Microchannels for Open-Channel Liquid Chromatography.</a> Anal. Chem. 75 (1), 57-70 (2003).</li><li>Blom, M. T., Chmela, E., Oosterbroek, R. E., Tijssen, R., van den Berg, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=On-Chip+Hydrodynamic+Chromatography+Separation+and+Detection+on+Nanoparticles+and+Biomolecules.">On-Chip Hydrodynamic Chromatography Separation and Detection on Nanoparticles and Biomolecules.</a> Anal. Chem. 75 (24), 6761-6768 (2003).</li><li>Betteridge, D., Fields, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Construction+of+pH+Gradients+in+Flow-Injection+Analysis+and+Their+Potential+Use+for+Multielement+Analysis+in+a+Single+Sample+Bolus.">Construction of pH Gradients in Flow-Injection Analysis and Their Potential Use for Multielement Analysis in a Single Sample Bolus.</a> Anal. Chem. 50 (4), 654-656 (1978).</li><li>Trojanowicz, M., Ko&#322;acin&#769;ska, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Recent+advances+in+flow+injection+analysis.">Recent advances in flow injection analysis.</a> Analyst. 141, 2085-2139 (2016).</li><li>Ajdari, A., Bontoux, N., Stone, H. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Hydrodynamic+Dispersion+in+Shallow+Microchannels:+The+Effect+of+Cross-Sectional+Shape.">Hydrodynamic Dispersion in Shallow Microchannels: The Effect of Cross-Sectional Shape.</a> Anal. Chem. 78 (2), 387-392 (2006).</li><li>Dutta, D., Ramachandran, A., Leighton, T. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Effect+of+channel+geometry+on+solute+dispersion+in+pressure-driven+microfluidic+systems.">Effect of channel geometry on solute dispersion in pressure-driven microfluidic systems.</a> Microfluid Nanofluid. 2 (4), 275-290 (2006).</li><li>Bontoux, N., Pe&#769;pin, A., Chen, Y., Ajdari, A., Stone, H. 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+characterization+of+hydrodynamic+dispersion+in+shallow+microchannels.">Experimental characterization of hydrodynamic dispersion in shallow microchannels.</a> Lab Chip. 6, 930-935 (2006).</li><li>Vedel, S., Bruus, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Transient+Taylor-Aris+dispersion+for+time-dependent+flows+in+straight+channels.">Transient Taylor-Aris dispersion for time-dependent flows in straight channels.</a> J. Fluid Mech. 691, 95-122 (2012).</li><li>Aminian, M., Bernardi, F., Camassa, R., McLaughlin, R. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Squaring+the+Circle:+Geometric+Skewness+and+Symmetry+Breaking+for+Passive+Scalar+Transport+in+Ducts+and+Pipes.">Squaring the Circle: Geometric Skewness and Symmetry Breaking for Passive Scalar Transport in Ducts and Pipes.</a> Phys. Rev. Lett. 115, 154503 (2015).</li><li>Aminian, M., Bernardi, F., Camassa, R., Harris, D. M., McLaughlin, R. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=How+boundaries+shape+chemical+delivery+in+microfluidics.">How boundaries shape chemical delivery in microfluidics.</a> Science. 354 (6317), 1252-1256 (2016).</li><li>Taylor, G. I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Dispersion+of+soluble+matter+in+solvent+flowing+slowly+through+a+tube.">Dispersion of soluble matter in solvent flowing slowly through a tube.</a> P Roy Soc Lond A Mat. 219 (1137), 186-203 (1953).</li><li>Stone, H. A., Stroock, A. D., Ajdari, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Engineering+Flows+in+Small+Devices:+Microfluidics+Toward+a+Lab-on-a-Chip.">Engineering Flows in Small Devices: Microfluidics Toward a Lab-on-a-Chip.</a> Annu. Rev. Fluid Mech. 36, 381-411 (2004).</li><li>Davis, M. E., Davis, 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=.">.</a> Fundamentals of chemical reaction engineering. , (2003).</li><li>Barton, N. <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+method+of+moments+for+solute+dispersion.">On the method of moments for solute dispersion.</a> J. Fluid Mech. 126, 205 (1983).</li></ol>

The authors have nothing to disclose.

After injecting dye into the pipe, the bolus is transported away from the injection needle using a steady flow. Then, it is necessary to wait long enough for the dye to diffuse across the cross-section of the channel. In this way, a uniform Gaussian-like distribution is obtained and will serve as the initial condition for the experiment. Hence, a laminar background flow is created with the programmable syringe pump. The experimental run lasts for 5 min with photos taken every second.
The most common issues in the setup come from the connection of the parts and the pipes. The various 3D-printed parts need to be sealed properly when connected to avoid leaking. The glass pipes are very delicate and must be handled and installed with care.
An issue we encountered when transitioning from the thin rectangular pipe to the thick square pipe was related to the fact that the pipe volume was reduced by a factor of 10. To maintain the same cross-sectional average flow speed with the mounted 12 mL syringe, the plunger speed in syringe pump A would have needed to be extremely low. At this programmed speed, the plunger velocity was not uniform anymore and a steady flow cannot be guaranteed throughout the experimental run. Therefore, we switched to a much smaller 1 mL syringe when working with the thick square pipe in step 2.5.1.
Also, one should verify that the average intensity along the vertical dimension of the pipe in the initial condition is approximately uniform. If not, a filtering mask needs to be applied across all frames to account for this discrepancy.
The least repeatable part of the experiment is the dye injection (and consequently the width of the initial distribution). As illustrated earlier, it is not a concern for matching with the Monte Carlo simulations, as the experimental initial condition can be recreated using the analysis of the initial photograph. The dye injection and consequent manual withdrawing may not always produce dye plugs of precisely the same width. Particular care needs to be applied when setting up the initial dye bolus. The experiment becomes more repeatable as researchers gain experience in this part of the protocol, but future improvements could certainly be made.
When comparing the setup with microfluidic devices, the only parameter that appear in the governing equation when appropriately nondimensionalized is the P&#233;clet number Pe if the tracer is passive, i.e. the tracer evolution is uncoupled from the flow. Dynamic similarity is implicit in the assumption of low Reynolds (Re &#60;&#60; 1) which ensures stable laminar flows u(y,z). These two parameters are setting the full similarity between microfluidic setups and the scales of our experiment. In practice, the physical length of the pipe only restricts the nondimensional times we can safely reach with our setup. For very late non-dimensional times, the necessary length of the pipe could become prohibitively long for a fixed P&#233;clet number in this large-scale setup.
An obvious limitation of this experimental protocol is that the data collected is a projected 2D representation of 3D geometry as the pictures are taken top-down on the pipe. The current process only allows to obtain the evolution of the cross-sectionally averaged dye distribution. Obtaining a distribution defined at each location in the tube rather than on its cross-sectional average and comparison with theoretical and numerical predictions are the subject of ongoing research.
All the experimental setup parts have technical drawings available for download which makes the setup easily accessible and customizable by any interested researcher. Building on the current results, the same setup will be used to study more complex and unexplored pipe geometries as well as different flow regimes.

In recent years, substantial efforts have been focused on developing microfluidic and lab-on-chip devices that can reduce the costs and increase the productivity of chemical preparation and diagnostics for a range of applications. One of the main features of microfluidic devices is the pressure-driven transport of fluids and dissolved solutes through microchannels. In this context, it has become increasingly important to better understand the controlled delivery of solutes at the microscale. In particular, applications such as chromatographic separation1,2&#160;and microfluidic flow injection analysis3,4&#160;require improved control and understanding of solute delivery. Researchers in microfluidics have studied and documented the influence of the channel&#39;s cross-sectional shape on solute spreading5,6,7,8, and the role of the channel&#39;s aspect ratio9,10.
Analytical and numerical studies of solute spreading along channels have recently lead to the identification of a correlation between the pipe cross-sectional geometry and the shape of the distribution9,10. At early timescales, the distribution strongly depends on the geometry: rectangular pipes break symmetry almost immediately, while elliptical pipes retain their initial symmetry much longer9. On the other hand, progressing into longer timescales the asymmetries in the solute distribution no longer differentiate ellipses from rectangles, and are set solely by the cross-sectional aspect ratio &#955; (ratio of the short to long side). Considering &#34;pipes&#34; of elliptical cross-sections and &#34;ducts&#34; of rectangular cross-sections, predictions from numerical simulations and asymptotic analysis were benchmarked with laboratory experiments. Thin channels (aspect ratio &#60;&#60; 1) produce solutes arriving with sharp fronts and tapering tails, while thick channels (aspect ratio ~ 1) present the opposite behavior10. This robust effect is relatively insensitive to the initial conditions and can be used to help select the solute distribution profile required for any application.
The behavior outlined above of sorting thin versus thick domains happens before the classical &#34;Taylor dispersion&#34; regime is reached. Taylor dispersion refers to the enhanced spreading of passive solutes in laminar flow (stable at low Reynolds number, Re) with a boosted effective diffusivity, inversely proportional to the solute&#39;s molecular diffusivity &#954;11. This enhancement is observed only after long, diffusive timescales, when the solute has diffused across the channel. Such diffusive timescale is defined in terms of the characteristic length scale a of the geometry, as td = a2/&#954;. The P&#233;clet number is a nondimensional parameter which measures the relative importance of fluid advection to diffusion effects. We define this parameter in terms of the shortest length scale as Pe = Ua/&#954;, where U is the characteristic flow speed. (The Reynolds number can be defined in terms of the P&#233;clet number as Re = Pe &#954;/&#957;, where&#160;&#957; is the kinematic viscosity of the fluid.) Typical P&#233;clet number values for microfluidic applications12&#160;vary between 10 and 105, with molecular diffusivities ranging from 10-7 to 10-5 cm2/s. Hence, given the flow speeds and length scales of interest, it is critical to understand the behavior of solutes for intermediate-to-long timescales (relative to the diffusive timescale), well past the initial observations of geometry-driven behavior and into the cross-section-driven regimes universal for a large class of geometries.
Given the interest in microfluidic applications, the choice of a large scale experimental setup may at first seem unnatural. The experiments reported herein are at the millimeter scale, not at the microscale as in true microfluidic devices. However, the same physical behaviors characterize both systems and a quantitative study of the relevant phenomena can still be achieved by properly scaling the governing equations, just as scale models of aircraft are assessed in wind tunnels during the design phase. In particular, matching relevant nondimensional parameters (such as the P&#233;clet number for our experiment) ensures the adaptability of the experimental model. Working at such larger scales, while maintaining a laminar pressure-driven flow, offers several advantages over a traditional microscale setup. In particular, the equipment required to manufacture, perform, and visualize the present experiments is easier to operate and less costly. Furthermore, other common challenges of working with microchannels, such as frequent clogging and the enhanced influence of manufacturing tolerances, are mitigated with the larger setup. Another possible use for this experimental setup is for studies of residence time distribution (RTD) in laminar flows13.
The asymmetries arising in the solute distribution downstream can be analyzed via its statistical moments; in particular, the skewness, which is defined as the centered, normalized third moment, is the lowest order integral statistic measuring the asymmetry of a distribution. The sign of the skewness typically indicates the shape of the distribution, i.e. if it is front-loaded (negative skewness) or back-loaded (positive skewness). Focusing on the aspect-ratios of the channels, there exists a clear correlation of thin geometries with front-loaded distributions, and thick geometries with back-loaded distributions10. Additionally, a critical aspect ratio separating these two opposite behaviors can be calculated for both elliptical pipes and rectangular ducts. Such crossover aspect ratios are remarkably similar for standard geometries, in particular, &#955;* = 0.49031 for pipes, and &#955;* = 0.49038 for ducts, suggestive of the universality of the theory10.
The experimental setup and method described in this paper are used to study the spreading of a pressure-driven passive solute in laminar fluid flows throughout glass capillaries of various cross-sections. The simplicity and reproducibility of the experiment defines a robust method of analysis for understanding the connection between a pipe&#39;s geometrical cross-section and the resulting shape of the injected solute distribution as it is transported downstream. The method discussed in this work has been developed to readily benchmark mathematical and numerical results in a physical laboratory setting.
A simple experimental procedure is described which highlights the definitive role played by a fluidic channel&#39;s cross-sectional aspect-ratio in setting the shape of a solute distribution downstream. The experimental setup requires a programmable syringe pump to produce a laminar steady flow, smooth glass pipes of various cross-sections, a second syringe pump to inject the diffusing solute (e.g. fluorescein dye) into the surrounding laminar flow, and UV-A lights and a camera to record the solute evolution. CAD files are provided for all the custom parts of the setup and such files can be used to 3D-print the experimental parts prior to assembly.

<table border="0" cellpadding="0" cellspacing="0" width="1013">
	<tbody>
		<tr height="36">
			<td height="36" style="height:36px;width:517px;">Flourescein Dye</td>
			<td style="width:175px;">Flinn Scientific&#160;</td>
			<td style="width:127px;">LOT: 118362&#160;&#160;&#160;&#160;&#160;&#160; CAS NO: 518-47-8</td>
			<td style="width:195px;"></td>
		</tr>
		<tr height="41">
			<td height="41" style="height:41px;">PhD ULTRA Hpsi Syringe Pump</td>
			<td>Harvard Apparatus</td>
			<td>703111</td>
			<td style="width:195px;">programmable digital syringe pump</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Compact Infusion Pump Model 975</td>
			<td>Harvard Apparatus</td>
			<td>55-1689</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Form 2 SLA 3D Printer</td>
			<td>Formlabs</td>
			<td>100-240</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Glass pipes</td>
			<td>VitroCom</td>
			<td>4410 and 8100</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">PTFE sealing tape</td>
			<td>Teflon</td>
			<td>4934A12</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">PVC tubing (1/8&#34; ID)</td>
			<td>McMaster</td>
			<td>5231K144</td>
			<td>5 Foot Length</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Reusable Stainless Steel Dispensing Needle 22 Gauge, .016&#34; ID, .028&#34; OD, 1/8&#34; NPT Thrd, 2&#34; Lg&#160;</td>
			<td>McMaster</td>
			<td>7590A45&#160;</td>
			<td>1 Required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">RTV silicone rubber sealant</td>
			<td>McMaster</td>
			<td>74945A69</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Plastic Syringe Manual, w/ Luer Lock Connection, .34 oz Capacity, Packs of 10&#160;</td>
			<td>McMaster</td>
			<td>7510A653&#160;</td>
			<td>1 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Plastic Syringe Manual, w/ Luer Slip Connection, .034 oz Cap, Packs of 10&#160;</td>
			<td>McMaster</td>
			<td>7510A603&#160;</td>
			<td>1 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Plastic Syringe Manual, w/ Luer Lock Connection, 0.1 oz Capacity, Packs of 10&#160;</td>
			<td>McMaster</td>
			<td>7510A651&#160;</td>
			<td>2 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Plastic dispensing tip</td>
			<td>McMaster</td>
			<td>6699A1&#160;</td>
			<td>3 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">6&#34; C-Clamps</td>
			<td>McMaster</td>
			<td>5133A18</td>
			<td>2 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Type 18-8 Stainless Steel Flat Washer Number 6 Screw Size, 0.156&#34; ID, 0.312&#34; OD, Packs of 100&#160;</td>
			<td>McMaster</td>
			<td>92141A008&#160;</td>
			<td>8 required</td>
		</tr>
		<tr height="60">
			<td height="60" style="height:60px;">18-8 SS Pan Head Phillips Machine Screw 6-32 Thread, 2-1/4&#34; Length, Packs of 50&#160;</td>
			<td style="width:175px;">McMaster</td>
			<td>91772A167&#160;</td>
			<td>4 required</td>
		</tr>
		<tr height="65">
			<td height="65" style="height:66px;">Oil-Resistant Buna-N Multipurpose O-Ring 1/16 Fractional Width, Dash Number 016, Packs of 100&#160;</td>
			<td style="width:175px;">McMaster</td>
			<td>9452K6&#160;</td>
			<td>3 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">Type 18-8 Stainless Steel Hex Nut 6-32 Thread Size, 5/16&#34; Wide, 7/64&#34; High, Packs of 100&#160;</td>
			<td style="width:175px;">McMaster</td>
			<td>91841A007&#160;</td>
			<td>4 required</td>
		</tr>
		<tr height="65">
			<td height="65" style="height:66px;">18-8 SS Pan Head Phillips Machine Screw 6-32 Thread, 1/2&#34; Length, Packs of 100&#160;</td>
			<td style="width:175px;">McMaster</td>
			<td>91772A148&#160;</td>
			<td>4 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;">24&#34; Black Light Fixture with bulb</td>
			<td style="width:175px;">American DJ</td>
			<td>B0002F5544</td>
			<td>2 required</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;width:517px;">DSLR camera&#160;</td>
			<td style="width:175px;">Nikon&#160;</td>
			<td>D300</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;width:517px;">24-120 mm lens</td>
			<td style="width:175px;">Nikon</td>
			<td>2193</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;width:517px;">Remote programmable trigger</td>
			<td style="width:175px;">Nikon</td>
			<td>4917</td>
			<td style="width:195px;">remote programmable trigger</td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;width:517px;">Memory Card</td>
			<td style="width:175px;">SanDisk&#160;</td>
			<td>SDCFX-032G-E61</td>
			<td></td>
		</tr>
		<tr height="22">
			<td height="22" style="height:23px;width:517px;">Metric ruler</td>
			<td style="width:175px;">McMaster</td>
			<td>20345A35</td>
			<td></td>
		</tr>
	</tbody>
</table>

the diffusion of passive tracers in laminar shear flow

A simple method to experimentally observe and measure the dispersion of a passive tracer in a laminar fluid flow is described. The method consists of first injecting fluorescent dye directly into a pipe filled with distilled water and allowing it to diffuse across the cross-section of the pipe to obtain a uniformly distributed initial condition. Following this period, the laminar flow is activated with a programmable syringe pump to observe the competition of advection and diffusion of the tracer through the pipe. Asymmetries in the tracer distribution are studied and correlations between the pipe cross-section and the shape of the distribution is shown: thin channels (aspect ratio &#60;&#60; 1) produce tracers arriving with sharp fronts and tapering tails (front-loaded distributions), while thick channels (aspect ratio ~1) present the opposite behavior (back-loaded distributions). The experimental procedure is applied to capillary tubes of various geometries and is particularly relevant to microfluidic applications by dynamical similarity.

The experimental setup after assembly is shown in Figure 1. Images produced in MATLAB show the experimental data above the processed evolution of the concentration curve (Figure 2) for three non-dimensional times. We have verified that there is a linear relation between the tracer&#39;s intensity and concentration. The shape of the distribution changes as time passes and the dye bolus moves downstream. Figure 2 shows such evolution in the case of the thin rectangular duct geometry. The initial dye distribution is narrow and symmetric (Gaussian-like with respect to the longitudinal direction and nearly uniform in the cross-section, Figure 2 left), but the symmetry is broken almost immediately as the background flow starts. The distribution breaks symmetry by presenting a sharp front and long tapering tails (Figure 2, middle and right).
The experimental results are confirmed by Monte Carlo simulations performed matching the initial distribution and flow rate (Figure 3). The fitted value for the dye diffusivity &#954; was determined in an independent experiment (step 2.4 in protocol) and used in this comparison. Monte Carlo methods are often used to simulate the evolution of advection-diffusion problems involving complex geometries as the boundary conditions (homogeneous Neumann in this case) can be simply input as billiard like reflection rules. The approach is to sample realizations of the equivalent stochastic differential equation underlying the advection-diffusion equation in nondimensional form:
<img alt="Equation 1" src="/files/ftp_upload/57205/57205eq1v2.jpg" />
where T(x,y,z,t) is the tracer distribution, &#964; is the nondimensional time normalized by td, x is the longitudinal spatial coordinate, y is the short transverse coordinate, and z is the long transverse coordinate, all normalized by the short side a. The fluid flow u(y,z) is the laminar steady-state solution to the Navier-Stokes equations with no-slip boundary conditions (no flow at the wall), driven by a negative pressure gradient. A Gaussian initial data in the pipe longitudinal direction with a desired variance can be obtained by considering only diffusion (Pe = 0) and evolving the particles for the desired time to match the width of the experimental initial data9,10. These representative results were obtained using the flow rate values specified in the protocol, however we expect the loading phenomena observed to hold in general for the laminar regime10&#160;(Figure 3).
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/57205/57205fig1.jpg" /> 
Figure 1: Experimental setup. (A) Diagram of the experimental setup. This figure has been modified from Aminianet al.10. (B) Presentation of the actual setup. <a href="https://cloudflare.jove.com/files/ftp_upload/57205/57205fig1large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/57205/57205fig2.jpg" /> 
Figure 2: Snapshots of processed data at various times. Top row: photo of the dye concentration diffused along the cross-section of the tube observed normally to the long cross-sectional direction at increasing non-dimensional times. The vertical axis has been scaled 5 times for the sake of clarity. Bottom: intensity of the dye concentration computed summing along the long cross-sectional direction. The peak value is normalized. <a href="https://cloudflare.jove.com/files/ftp_upload/57205/57205fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/57205/57205fig3.jpg" /> 
Figure 3: Comparison of the concentration distribution between Monte Carlo simulations and experiments. The evolution of the cross-sectionally averaged dye concentration along the longitudinal length of the pipe is shown at two instants in time: &#964; = 0.15&#160; and &#964; = 0.30. The dashed lines are the simulation results, while the solid lines represent the experimental data. Top: comparison in the thick (square) channel; bottom: comparison in the thin (rectangular) channel. The area under each curve is normalized to be one and x = 0 corresponds to the center of the initial plug of dye. This figure has been modified from Aminianet al.10. <a href="https://cloudflare.jove.com/files/ftp_upload/57205/57205fig3large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Supplementary File 1. Included CAD drawings of 3D Hexagonal Connector (<a href="https://cloudflare.jove.com/files/ftp_upload/57205/hex_connector_3D.STL" target="_blank">hex_connector_3D.STL</a>)
Supplementary File 2. Included CAD drawings of 3D Injector Post (<a href="https://cloudflare.jove.com/files/ftp_upload/57205/injector_post_3D.STL" target="_blank">injector_post_3D.STL</a>)
Supplementary File 3. Included CAD drawings of 3D Reservoir (<a href="https://cloudflare.jove.com/files/ftp_upload/57205/reservoir_3D.STL" target="_blank">reservoir_3D.STL</a>)
Supplementary File 4. Included CAD drawings of 3D Thick Pipe Plates (<a href="https://cloudflare.jove.com/files/ftp_upload/57205/plate_thick_3D.STL" target="_blank">plate_thick_3D.STL</a>)
Supplementary File 5. Included CAD drawings of 3D Thin Pipe Plates (<a href="https://cloudflare.jove.com/files/ftp_upload/57205/plate_thin_3D.STL" target="_blank">plate_thin_3D.STL</a>)

Watch this Scientific Journal Video about The Diffusion of Passive Tracers in Laminar Shear Flow at JoVE.com

The Diffusion of Passive Tracers in Laminar Shear Flow

A protocol for the study of the diffusion of passive tracers in laminar pressure-driven flow is presented. The procedure is applicable to various capillary pipe geometries.

A simple method to experimentally observe and measure the dispersion of a passive tracer in a laminar fluid flow is described. The method consists of ...

Research

Education

JoVE Journal

JoVE Core

Mechanical Engineering

Anatomy and Physiology

Engineering

Unit 1

Shearing Stresses in Beams and Thin-Walled Members

Cell Membrane Structure and Functions

SCIENTISTS IN ACTION

Experimental Measurement of Settling Velocity of Spherical Particles in Unconfined and Confined Surfactant-based Shear Thinning Viscoelastic Fluids

An experimental study is performed to measure the terminal settling velocities of spherical particles in surfactant based shear thinning viscoelastic ...

Speciation and Bioavailability Measurements of Environmental Plutonium Using Diffusion in Thin Films

The biological uptake of plutonium (Pu) in aquatic ecosystems is of particular concern since it is an alpha-particle emitter with long half-life which can ...

Ultrasonic Welding of Thermoplastic Composite Coupons for Mechanical Characterization of Welded Joints through Single Lap Shear Testing

This paper presents a novel straightforward method for ultrasonic welding of thermoplastic-composite coupons in optimum processing conditions. The ...

In Vitro Culture of Epicardial Cells From Mouse Embryonic Heart

In Vitro Culture of Epicardial Cells From Mouse Embryonic Heart

During embryogenesis, the epicardial contribution to coronary vasculature development has been very well established. Cells derived from the epicardium ...

Experimental Methodology for Estimation of Local Heat Fluxes and Burning Rates in Steady Laminar Boundary Layer Diffusion Flames

Modeling the realistic burning behavior of condensed-phase fuels has remained out of reach, in part because of an inability to resolve the complex ...

Design and Use of a Full Flow Sampling System (FFS) for the Quantification of Methane Emissions

Design and Use of a Full Flow Sampling System FFS for the Quantification of Methane Emissions

The use of natural gas continues to grow with increased discovery and production of unconventional shale resources. At the same time, the natural gas ...

Microfluidic Imaging Flow Cytometry by Asymmetric-detection Time-stretch Optical Microscopy (ATOM)

Microfluidic Imaging Flow Cytometry by Asymmetric-detection Time-stretch Optical Microscopy ATOM

Scaling the number of measurable parameters, which allows for multidimensional data analysis and thus higher-confidence statistical results, has been the ...

Generation of Size-controlled Poly (ethylene Glycol) Diacrylate Droplets via Semi-3-Dimensional Flow Focusing Microfluidic Devices

Generation of Size-controlled Poly ethylene Glycol Diacrylate Droplets via Semi-3-Dimensional Flow Focusing Microfluidic Devices

Uniform and size-controllable poly (ethylene glycol) diacrylate (PEGDA) droplets could be produced via the flow focusing process in a microfluidic device. ...

Visually Based Characterization of the Incipient Particle Motion in Regular Substrates: From Laminar to Turbulent Conditions

Two different experimental methods for determining the threshold of particle motion as a function of geometrical properties of the bed from laminar to ...

Experimental Investigation of the Flow Structure over a Delta Wing Via Flow Visualization Methods

It is well known that the flow field over a delta wing is dominated by a pair of counter rotating leading edge vortices (LEV). However, their mechanism is ...