Confined Fluid Flow:

Microfluidics and Capillarity

Douglas P. Holmes
Mechanical Engineering
Boston University

Sapienza, Università di Roma – Short Course (2015)

 Why Small?
Microfluidic chemostat
…static chemical environment used as bioreactor…
fresh medium is continuously added, culture liquid is continuously removed

Device provided single-cell resolution to study microbial population growth.

Balagaddé, Frederick K., Lingchong You, Carl L. Hansen, Frances H. Arnold, and Stephen R. Quake. "Long-term monitoring of bacteria undergoing
programmed population control in a microchemostat." Science 309, no. 5731 (2005): 137-140.
Whitesides, George M. "The origins and the future of microfluidics." Nature 442, no. 7101 (2006): 368-373.
 Why Small?
Protein cystallization

Droplets containing
proteins are trapped in
microchannel wells.

Dyed droplets in the wells

of the device.

Proteins crystallized within

droplets in a glass capillary.

Device used to efficiently screen for optimal protein crystallization conditions

Zheng, Bo, Joshua D. Tice, L. Spencer Roach, and Rustem F. Ismagilov. "A Droplet‐Based, Composite PDMS/Glass Capillary Microfluidic System for
Evaluating Protein Crystallization Conditions by Microbatch and Vapor‐Diffusion Methods with On‐Chip X‐Ray Diffraction." Angewandte chemie
international edition 43, no. 19 (2004): 2508-2511.
Whitesides, George M. "The origins and the future of microfluidics." Nature 442, no. 7101 (2006): 368-373.
 Why Small?
Bubbles and foams

Bubble generation with tunable, monodisperse sizes

Garstecki, Piotr, Irina Gitlin, Willow DiLuzio, George M. Whitesides, Eugenia Kumacheva, and Howard A. Stone. "Formation of monodisperse bubbles in a
microfluidic flow-focusing device." Applied Physics Letters 85, no. 13 (2004): 2649-2651.
Whitesides, George M. "The origins and the future of microfluidics." Nature 442, no. 7101 (2006): 368-373.
 Why Small?
Capillary origami

Using surface tension to deform thin structures

Py, Charlotte, Paul Reverdy, Lionel Doppler, José Bico, Benoit Roman, and Charles N. Baroud. "Capillary origami: spontaneous wrapping of a droplet with
an elastic sheet." Physical Review Letters 98, no. 15 (2007): 156103.
 Why Small?
Cellular and developmental biology

Study effect of temperature on the development of a fruitfly embryo.

Immobilized embryo has water with cold (left) and warm (right) flow over it.
Lucchetta, Elena M., Ji Hwan Lee, Lydia A. Fu, Nipam H. Patel, and Rustem F. Ismagilov. "Dynamics of Drosophila embryonic patterning network perturbed
in space and time using microfluidics." Nature 434, no. 7037 (2005): 1134-1138.
Whitesides, George M. "The origins and the future of microfluidics." Nature 442, no. 7101 (2006): 368-373.
 Why Small?
Inexpensive diagnostics

Costs as low as
$0.001 per device
(€0.00089)

Low-cost, simple paper-based microfluidics for diagnostics

Carrilho, Emanuel, Andres W. Martinez, and George M. Whitesides. "Understanding wax printing: a simple micropatterning process for paper-based
microfluidics." Analytical chemistry 81, no. 16 (2009): 7091-7095.
 Why Small?
Volume goes as L3
•  Small decrease in size = large
reduction of sample volume.

How much sample volume do we need?

Detection of biomolecules:
Digoxin – heart stimulating drug

Cortisol – stress hormone from adrenal gland

Creatinine – level in blood is measure of kidney

function

Theophylline – drug used to treat respiratory

diseases, e.g. asthma

http://homes.nano.aau.dk/lg/Lab-on-Chip2008_files/Lab-on-chip2008_1.pdf
 Microfluidics

Published items in each year… Citations in each year…

ISI Web of Knowledge: Citation Report “Topic: Microfluidics”, 1990-2009.

 Overview
Confined Fluid Flow:
Microfluidics and Capillarity
Reynolds Number: Inertia vs. Viscous effects
•  Review of characteristic flows…

Péclet Number: Transport phenomena in a continuum

•  Diffusion, separation, and mixing…

Geometric confinement: Controlling and manipulating fluid flow

•  Microfluidic fabrication, valving, pumping…

Capillary Number: Viscosity vs. Surface tension

•  Droplet formation, capillary rise, elasticity…
 Fluid Dynamics
Navier-Stokes Equations (momentum conservation)

Continuity Equation (mass conservation)

 Fluid Behavior
Inertial Forces Viscous Forces

Trailing airplane vortices Coiling honey

Reynolds Number: inertial/viscous

Airplane: http://eis.bris.ac.uk/~glhmm/gfd/Airplane-ChrisWillcox.jpg
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
 Fluid Behavior
Inertial Forces Reynolds Number: inertial/viscous

Viscous Forces

Fluid element accelerating around curve. Fluid element in a channel of

•  During a turn time: contracting length.
•  By mass conservation,
•  Loss of momentum density: velocity increases as:

•  By exerting an inertial centrifugal •  Gain momentum at a rate:

force density:

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
Airplane: http://eis.bris.ac.uk/~glhmm/gfd/Airplane-ChrisWillcox.jpg
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
 Fluid Behavior
Inertial Forces Reynolds Number: inertial/viscous

Viscous Forces

Inertial Forces

Viscous Forces
•  Viscous force densities result from gradients in viscous stress:

Ratio of these two

force densities is the
Reynolds number:
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
Airplane: http://eis.bris.ac.uk/~glhmm/gfd/Airplane-ChrisWillcox.jpg
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
 Fluid Behavior
Inertial Forces Reynolds Number: inertial/viscous
Estimation of Reynolds numbers for common microfluidic devices.
•  Typical fluid – water
•  Viscosity: 1.025 cP @ 25°C
•  Density: 1 g/mL
•  Typical channel dimensions
•  Radius/height (smaller than width): 1 – 100 µm
•  Typical velocities
•  Average velocity: 1 µm/s – 1 cm/s
Viscous Forces
Typical Reynolds number:

Low Reynolds number: viscous forces > inertial forces

•  Flows are linear.
•  Nonlinear terms in Navier-Stokes disappear
•  Linear, predictable Stokes flow

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
Airplane: http://eis.bris.ac.uk/~glhmm/gfd/Airplane-ChrisWillcox.jpg
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
 Fluid Behavior
Inertial Forces Reynolds Number: inertial/viscous

Reynolds number is focused on steady inertial

forcing due to the convective derivative:
Viscous Forces •  Almost always unimportant in microfluidic flows.
•  Time dependence arises from U0

Linear unsteady term sets the inertial time

scale to establish steady flows:

Time scale estimated by balancing unsteady

inertial force density with viscous force density

Time required for a vorticity to diffuse a distance L0, with a diffusivity ν=µ/ρ, τi~10ms

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
Airplane: http://eis.bris.ac.uk/~glhmm/gfd/Airplane-ChrisWillcox.jpg
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
 Fluids and Circuits

Pressure-driven flow in a network of parallel Electrical circuit with resistances in parallel.

channels. •  First elementary rule of circuit design
•  Flow rates in three parallel channels. is Ohm’s law
•  Total flow rate: •  Relates electrical potential to
•  Pressure drop and flow rate related by: current:

Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Cylindrical Tube
Continuity Equation (cylindrical coords):

Flow only in z direction:

Navier-Stokes equations:

Assumptions: Reduces to the following:

1.  Neglect gravity:

2.  Steady state: Pressure driven flow – linear change

3.  Axisymmetric:
Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Cylindrical Tube
Navier-Stokes equations - cylindrical pipe:

…integrate once…

Boundary conditions:

1.  Velocity is at a maximum at the

center of the piper:

…integrate twice…
2.  No-slip along the walls.

Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Cylindrical Tube
Parabolic velocity distribution:

Flow rate (volume/time)

The flow rate depends on the fourth power of the radius: important in small systems
•  e.g. For the same pressure gradient, reducing the radius by a factor of two causes
a 16-fold reduction in flow rate.
•  e.g. Consider a blood vessel – a 10% decrease in radius produces more than a
40% decrease in the flow rate of blood (PSA: Eat more kale!)

Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Cylindrical Tube
Average velocity

Recall the Reynolds number:

For a fixed pressure gradient, and so a factor of two change in radius produces
a factor of eight change in Reynolds number.

The average velocity is proportional to the pressure gradient, writing the flow field in
vector form gives Darcy’s law:
- k is the permeability (dimensions length^2)
- linear relation between pressure and velocity
is due to lack of inertia effects
- order-of-magnitude of permeability is
typically the square of smallest dimension.

Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Cylindrical Tube
Hydrodynamic resistance

recall:

Electrical analogy: resistance in fluid

channels depends on the fourth power of
the radius rather than the second power
in the electrical case.
…………………small matters!

Scaling on physical & dimensional arguments…

•  Can arrive at: and from
dimensional analysis.

•  Fluid motion arises as a balance between the pressure drop driving

motion and the frictional (viscous) resistance from the bounding walls.
This balance applies independent of the cross-sectional channel shape.
Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Rectangular
Parabolic velocity distribution:

Average velocity:
Approximate flow rate (channel of width w)

Hydrodynamic
Resistance:

Parabolic velocity distribution (known analytically in terms of a Fourier series):

where: and, from the no-slip boundary conditions:

Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Rectangular
Parabolic velocity distribution:

Average velocity:
Approximate flow rate (channel of width w)

Hydrodynamic
Resistance:

Corresponding flow rate

where:

Flow rate is nearly

linear in aspect ratio
h/w

Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Flow: Rectangular

Flow rate in a rectangular channel

Approximate dimensionless flow rate

Stone, Howard A. "Introduction to fluid dynamics for microfluidic flows." In CMOS Biotechnology, pp. 5-30. Springer US, 2007.
 Workout Problem
While developing a microfluidic chip for your research, you decide a scale model would
help you gain intuition about flow characteristics at the microscale. While most scale
models are smaller that the original, your model will be 100 times larger (so that you can
see it without a microscope).

Viscosity Density
Honey 8750 cP 1.42 g/mL
Olive oil 81 cP 0.92 g/mL
Water 1.025 cP @ 25°C 1 g/mL
Molasses 50,000 cP 1.50 g/mL

For the microfluidic chip, the flow conditions are: fluid velocity 1 mm/s, channel diameter
100 µm, and the liquid is water at T = 25°C. The Reynolds number for both systems is to be
duplicated, and the microchannel is hemicylindrical (use the hydraulic diameter). For the
scale model, which of the condiments listed above from your kitchen would you choose
as the flow liquid, if the velocity in the model is:

a.) 1 mm/s
Recall: hydraulic diameter, D H ,
b.) 60 mm/s relates flow in noncircular channels
c.) 340 mm/s to those in a circular channel:

http://www-bsac.eecs.berkeley.edu/projects/ee245/Homework/
 Workout Problem
Hydraulic diameter of hemicylindrical channel:

Reynold’s number of microfluidic chip:

a. Fluid flowing at 1 mm/s in scaled up model

Use Olive Oil.

b. Fluid flowing at 60 mm/s in scaled up model

Use Honey.

c. Fluid flowing at 340 mm/s in scaled up model

Use Molasses.

http://www-bsac.eecs.berkeley.edu/projects/ee245/Homework/
 Overview
Confined Fluid Flow:
Microfluidics and Capillarity
Reynolds Number: Inertia vs. Viscous effects
•  Review of characteristic flows…

Péclet Number: Transport phenomena in a continuum

•  Diffusion, separation, and mixing…

Geometric confinement: Controlling and manipulating fluid flow

•  Microfluidic fabrication, valving, pumping…

Capillary Number: Viscosity vs. Surface tension

•  Droplet formation, capillary rise, elasticity…
 Fluid Behavior
Convection Diffusion

Convection current as hot water Diffusion of salt ions while curing pork
mixes with cold water

Péclet Number: convection/diffusion

http://quarkyscience.ca/project-of-the-month/convection-demonstration/
http://www.genuineideas.com/ArticlesIndex/diffusion.html
 Fluid Behavior
Convection Péclet Number: convection/diffusion
High Reynolds: eddies chaotically stretch and fold fluid elements.
•  Turbulent mixing & thermal convection move fluids.

Low Reynolds: mixing occurs by diffusion only

•  Diffusion is remarkably slow – long mixing times.

T-Junction in which two

fluids are injected and flow
Diffusion
alongside each other

How far down the channel must the

fluids flow before the channel is
homogenized?

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
http://quarkyscience.ca/project-of-the-month/convection-demonstration/
http://www.genuineideas.com/ArticlesIndex/diffusion.html
 Fluid Behavior
Convection Péclet Number: convection/diffusion

How far down the channel must the

fluids flow before the channel is
homogenized?
Time for particles/molecules to diffuse across channel:

Diffusion

During this time, the stripe moves down the channel:

Number of channel widths required for complete mixing:

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
http://quarkyscience.ca/project-of-the-month/convection-demonstration/
http://www.genuineideas.com/ArticlesIndex/diffusion.html
 Fluid Behavior
Convection Péclet Number: convection/diffusion

Relative importance: convection to diffusion:

Diffusion

Consider a small protein flowing with liquid:

•  Typical size: 5 nm
•  Diffusion constant: 40 µm2/s
•  e.g. Channel width: 100 µm
•  e.g. Flow velocity: 100 µm/s

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
http://quarkyscience.ca/project-of-the-month/convection-demonstration/
http://www.genuineideas.com/ArticlesIndex/diffusion.html
 Sensing & Filtering
Parallel Laminar Flows
Pressure-driven laminar flow of two adjacent
miscible streams.
•  Solution on left contains calcium.
•  Solution on right contains calcium-dependent
fluorophore, Fluo-3.
•  In water, Fluo-3 and calcium form a fluorescent
complex at a diffusion limited rate.

Convective-diffusion equation:

For high Péclet number flows,

•  Diffusive broadening down the channel depends
on both y & x

Away from boundaries, the interfacial region:

Ismagilov, Rustem F., Abraham D. Stroock, Paul JA Kenis, George Whitesides, and Howard A. Stone. "Experimental and theoretical scaling laws for
transverse diffusive broadening in two-phase laminar flows in microchannels." Applied Physics Letters 76, no. 17 (2000): 2376-2378.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Sensing & Filtering
Parallel Laminar Flows
Pressure-driven laminar flow of two adjacent
miscible streams.
•  Solution on left contains calcium.
•  Solution on right contains calcium-dependent
fluorophore, Fluo-3.
•  In water, Fluo-3 and calcium form a fluorescent
complex at a diffusion limited rate.

Away from boundaries, the interfacial region:

Lévêque Problem: Near from boundaries, velocity

varies linearly with distance into the channel.
Diffusion across a linear flow field

Diffusion is the only time scale in both x & y:

Ismagilov, Rustem F., Abraham D. Stroock, Paul JA Kenis, George Whitesides, and Howard A. Stone. "Experimental and theoretical scaling laws for
transverse diffusive broadening in two-phase laminar flows in microchannels." Applied Physics Letters 76, no. 17 (2000): 2376-2378.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Sensing & Filtering
Parallel Laminar Flows
Pressure-driven laminar flow of two adjacent
miscible streams.
•  Solution on left contains calcium.
•  Solution on right contains calcium-dependent
fluorophore, Fluo-3.
•  In water, Fluo-3 and calcium form a fluorescent
complex at a diffusion limited rate.
Edge of channel
Away from boundaries, the interfacial region:

Lévêque Problem: Near from boundaries, velocity

varies linearly with distance into the channel.
Diffusion across a linear flow field

Center of channel

Diffusion is the only time scale in both x & y:

Ismagilov, Rustem F., Abraham D. Stroock, Paul JA Kenis, George Whitesides, and Howard A. Stone. "Experimental and theoretical scaling laws for
transverse diffusive broadening in two-phase laminar flows in microchannels." Applied Physics Letters 76, no. 17 (2000): 2376-2378.
Stone, Howard A., Abraham D. Stroock, and Armand Ajdari. "Engineering flows in small devices: microfluidics toward a lab-on-a-chip." Annu. Rev. Fluid
Mech. 36 (2004): 381-411.
 Sensing & Filtering
H Filter Parallel Laminar Flows
Filter particles by size without a membrane.
•  Pressure-driven laminar flow of two adjacent
miscible streams.
•  One stream is a dilute solution of different sized
particles.
•  Each particle has its own diffusivity and Péclet
number.

Péclet number determines the channel

length required for each component to
diffuse across the channel width.

H filter works best when one Péclet number

is large and the other is small.
•  e.g. small particles diffuse across
channel, large ones do not.

Brody, James P., Paul Yager, Raymond E. Goldstein, and Robert H. Austin. "Biotechnology at low Reynolds numbers." Biophysical journal 71, no. 6 (1996):
3430-3441.
Brody, James P., and Paul Yager. "Diffusion-based extraction in a microfabricated device." Sensors and Actuators A: Physical 58, no. 1 (1997): 13-18.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Sensing & Filtering
H Filter - Variant
Separation by Péclet number
•  Components to be separated spread across
channel at different rates.
•  e.g. Separation of motile vs. nonmotile sperm.

Motile sperm swim rapidly and randomly to

fill the channel, as compared to nonmotile
sperm.
•  Nonmotile sperm spread by diffusion alone.

Motile sperm purity at inlet (blue) vs. motile

sperm purity at outlet (purple)

Cho, Brenda S., Timothy G. Schuster, Xiaoyue Zhu, David Chang, Gary D. Smith, and Shuichi Takayama. "Passively driven integrated microfluidic system for
separation of motile sperm." Analytical Chemistry 75, no. 7 (2003): 1671-1675.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Sensing & Filtering
Fabrication using multiple laminar streams
…beating diffusion…

Fabrication of a three electrode system

•  Multiple streams & large Péclet number.
•  Minimal mixing over large distances.

Nonmixing, high Péclet flows can fabricate/

etch structures within microchannels.

e.g. Parallel streams selectively etch a gold

electrode
•  A silver wire is formed by a precipitation reaction
between the two streams.
•  Result: three-electrode system.

Kenis, Paul JA, Rustem F. Ismagilov, and George M. Whitesides. "Microfabrication inside capillaries using multiphase laminar flow patterning." Science 285,
no. 5424 (1999): 83-85.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Rapid mixing means doing better than diffusion.

Why mix?
•  Study chemical reaction kinetics…
•  Probe protein folding…

Fluid stirring will stretch and fold inhomogenous fluid elements until mixing occurs.
•  Mixing: diffusive migration across streamlines.
•  Stirring motions reduce distance over which mixing must occur.

General design principles for mixing:

Dispersion of tracers occurs first by convective stretching

with the fluid, followed by diffusive homogenization.

Tracers: an object transported by and diffuses into the fluid,

e.g. dyes, analyte molecules, proteins, cells, salt, or heat.

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Taylor Dispersion
Thin stripe of tracers spans circular channel
•  Channel radius: w

Pressure-driven (Poiseuille) flow stretches

strip into a parabola with profile:

Molecular diffusion across the channel

smears the parabola into a plug of width:

Each stripe is convectively stretched then

diffusively smeared
After N time steps of , the initial
stripe grows into a Gaussian with width:

Taylor, Geoffrey. "Dispersion of soluble matter in solvent flowing slowly through a tube." In Proceedings of the Royal Society of London A: Mathematical,
Physical and Engineering Sciences, vol. 219, no. 1137, pp. 186-203. The Royal Society, 1953.
Taylor, Geoffrey. "Conditions under which dispersion of a solute in a stream of solvent can be used to measure molecular diffusion." Proceedings of the
Royal Society of London. Series A. Mathematical and Physical Sciences 225, no. 1163 (1954): 473-477.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Taylor Dispersion
Initial stripe grows to a Gaussian with width:

The tracer distribution grows diffusively with

an effective long-time axial diffusivity:

Axial diffusivity occurs in addition to molecular diffusivity.

Relevant examples for specific geometries:

Taylor, Geoffrey. "Dispersion of soluble matter in solvent flowing slowly through a tube." In Proceedings of the Royal Society of London A: Mathematical,
Physical and Engineering Sciences, vol. 219, no. 1137, pp. 186-203. The Royal Society, 1953.
Taylor, Geoffrey. "Conditions under which dispersion of a solute in a stream of solvent can be used to measure molecular diffusion." Proceedings of the
Royal Society of London. Series A. Mathematical and Physical Sciences 225, no. 1163 (1954): 473-477.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Taylor Dispersion

Convective stretching enhances axial dispersivity

Taylor dispersivity is only valid at

•  Long time scales:

•  Downstream lengths:

Taylor dispersion acts in direction of flow.

•  Not observed in T sensors or H filters

Taylor, Geoffrey. "Dispersion of soluble matter in solvent flowing slowly through a tube." In Proceedings of the Royal Society of London A: Mathematical,
Physical and Engineering Sciences, vol. 219, no. 1137, pp. 186-203. The Royal Society, 1953.
Taylor, Geoffrey. "Conditions under which dispersion of a solute in a stream of solvent can be used to measure molecular diffusion." Proceedings of the
Royal Society of London. Series A. Mathematical and Physical Sciences 225, no. 1163 (1954): 473-477.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Rotary Mixer

Fluid pumped around a

circular channel:
•  radius R
•  height h
•  average velocity U0.

(i.) Diffusion-dominated (ii.) Taylor dispersion-mediated (iii.) Convectively Stirred

•  Mixing occurs when tracers •  Axial spreading increases •  At high flow rates, tracer stripes
diffuse around the circumference diffusivity with Taylor dispersivity: fold into themselves before
of the ring, at time: molecules diffuse across channel.

•  Independent of Péclet number. •  Dominant when molecules diffuse

over h before convection.
H.P. Chou, M.A. Unger, and S. Quake, Biomed. Microdevices, 3, 323, (2001).
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Chaotic advection

Even simple Stokes flows can have

chaotic streamlines that exponentially
Blinking vortex flow stretch and fold.

Steady, incompressible two-dimensional

flows are integrable and cannot exhibit
chaotic trajectories.

Steady three-dimensional flows can

have chaotic streamlines (as can
unsteady two dimensional flows).

Can occur in droplets by superposing

two simple flow fields
•  (e.g. sedimenting drop in a shear flow)

Aref, Hassan. "Stirring by chaotic advection." Journal of fluid mechanics 143 (1984): 1-21.
Aref, Hassan. "The development of chaotic advection." Physics of Fluids, 14,4 (2002): 1315-1325.
http://upload.wikimedia.org/wikipedia/en/b/b8/Blinking_vortex_flow.jpg
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Continuous-flow staggered herringbone mixer
Chaotic advection

Even simple Stokes flows can have

chaotic streamlines that exponentially
stretch and fold.

The staggered herringbone mixer is a

chaotic mixer for continuous flow systems
– independent of inertia.
•  Asymmetric grooves in channel walls induce
axially modulated secondary flow.
•  Counter-rotating fluid rolls.
•  Asymmetry is periodically reversed so that
the distance between stripes halves with
each cycle – exponential stretching/folding.

Aref, Hassan. "Stirring by chaotic advection." Journal of fluid mechanics 143 (1984): 1-21.
Aref, Hassan. "The development of chaotic advection." Physics of Fluids, 14,4 (2002): 1315-1325.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Mixing
Continuous-flow staggered herringbone mixer
Chaotic advection
Asymmetry pattern on channel causes
exponential stretching/folding.

After N cycles of time:

Stripes separated by distance:

Time to diffuse between stripes:

Mixing:
Therefore, the number of cycles:

Time for chaotic mixing:

Aref, Hassan. "Stirring by chaotic advection." Journal of fluid mechanics 143 (1984): 1-21.
Aref, Hassan. "The development of chaotic advection." Physics of Fluids, 14,4 (2002): 1315-1325.
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Overview
Confined Fluid Flow:
Microfluidics and Capillarity
Reynolds Number: Inertia vs. Viscous effects
•  Review of characteristic flows…

Péclet Number: Transport phenomena in a continuum

•  Diffusion, separation, and mixing…

Geometric confinement: Controlling and manipulating fluid flow

•  Microfluidic fabrication, valving, pumping…

Capillary Number: Viscosity vs. Surface tension

•  Droplet formation, capillary rise, elasticity…
Geometric Confinement:
Control & Manipulation of Fluid Flow
 A Giraffe’s Jugular

http://upload.wikimedia.org/wikipedia/commons/e/e6/Giraffes_at_west_midlands_safari_park.jpg
 A Giraffe’s Jugular

Top giraffe: http://cdn.asiancorrespondent.com/wp-content/uploads/2011/11/Giraf2-349x248.jpg

Bottom giraffe: http://www.onekind.org/uploads/a-z/az_giraffe1.jpg
 A Giraffe’s Jugular
Arterial and venous
pressure decrease upon
lowering the head.

Rise in cranial jugular

pressure coincides with
cessation of jugular flow.

Jugular venous flow

ceased for ~30s when
lowered.

Brøndum, E., et al. "Jugular venous pooling during lowering of the head affects blood pressure of the anesthetized giraffe." American Journal of
Physiology-Regulatory, Integrative and Comparative Physiology 297.4 (2009): R1058-R1065.
 A Giraffe’s Jugular
Jugular vein’s cross sectional area changes
upon lifting/lowering of the Giraffe's head.

Brøndum, E., et al. "Jugular venous pooling during lowering of the head affects blood pressure of the anesthetized giraffe." American Journal of
Physiology-Regulatory, Integrative and Comparative Physiology 297.4 (2009): R1058-R1065.
 A Giraffe’s Jugular

Saara Tuulia, http://factloops.blogspot.com/2012/01/short-animations-on-a-brief.html

 Flexible Tubes
Almost all vessels carrying fluids with the body are flexible.
Fluid-structure interactions between internal flow and tube deformation often
dictate a vessel’s biological function or dysfunction.

Grotberg, James B., and Oliver E. Jensen. "Biofluid mechanics in flexible tubes." Annual Review of Fluid Mechanics 36.1 (2004): 121.
 Flexible Tubes
Starling Resistor

Grotberg, James B., and Oliver E. Jensen. "Biofluid mechanics in flexible tubes." Annual Review of Fluid Mechanics 36.1 (2004): 121.
 Flexible Tubes
Starling Resistor (1D Model)

Mass conservation:

(Very) Reduced model

•  “Bench-top” model for a deformable airway
•  Consider airways as single, compliant tube
Momentum conservation: •  Prone to instabilities – best with low Reynolds

Compliance of tube gives rise to

Pressure-area relation:
fluid-structure interactions

O.E. Jensen, “Flows through deformable airways” (2002), see: http://www.biomatematica.it/urbino2002/programmi/oejnotes.pdf

Grotberg, James B., and Oliver E. Jensen. "Biofluid mechanics in flexible tubes." Annual Review of Fluid Mechanics 36.1 (2004): 121.
How do we translate these ideas
to microfluidic devices?
 Fabrication

Duffy, David C., et al. "Rapid prototyping of microfluidic systems in poly(dimethylsiloxane)." Analytical chemistry 70.23 (1998): 4974-4984.
 Fabrication

Duffy, David C., et al. "Rapid prototyping of microfluidic systems in poly(dimethylsiloxane)." Analytical chemistry 70.23 (1998): 4974-4984.
M.A. Burns, C.H. Mastrangelo, T.S. Sammarco, F.P. Man, J.R. Webster, et al., “Microfabricated structures for integrated DNA analysis.” Proc. Natl. Acad. Sci.
USA, 93:5556-5561, (1996).
 Fabrication

Duffy, David C., et al. "Rapid prototyping of microfluidic systems in poly(dimethylsiloxane)." Analytical chemistry 70.23 (1998): 4974-4984.
T. Thorsen, S.J. Maerkl, S.R. Quake, “Microfluidic large-scale integration.” Science, 298:580-584, (2002).
 Pneumatic Valves
“Quake” valve
•  Bilayer microfluidic chip.

•  Thin film separating two

flow channels.
•  One channel: Fluid
•  One channel: Air
(controlling)

•  Pressurized air deflects

the thin film and closes
the fluid channel.

Unger, Marc A., Hou-Pu Chou, Todd Thorsen, Axel Scherer, and Stephen R. Quake. "Monolithic microfabricated valves and pumps by multilayer soft
lithography." Science 288, no. 5463 (2000): 113-116.
 Pneumatic Valves
“Quake” valve
•  Bilayer microfluidic chip.

•  Thin film separating two

flow channels.
•  One channel: Fluid
•  One channel: Air
(controlling)

•  Pressurized air deflects

the thin film and closes
the fluid channel.
 Pneumatic Valves
“Quake” valve

Unger, Marc A., Hou-Pu Chou, Todd Thorsen, Axel Scherer, and Stephen R. Quake. "Monolithic microfabricated valves and pumps by multilayer soft
lithography." Science 288, no. 5463 (2000): 113-116.
 Pneumatic Valves
“Quake” valve
Push-down Valves
•  Control lines above flow channels.
•  Pressure flattens membrane valve down to seal.
•  Suitable for low aspect ratio (1:10) & shallow (~10um)
channels.
•  Flow geometry: 100um x 13um
•  Control geometry: 100um x 10-25um
•  Applications: where fluid flow must be in contact with
substrate (spotting DNA, patterned substrate, etc.)

Push-up Valves
•  Control lines belowflow channels.
•  Pressure deflect membrane valve up to seal.
•  Flow geometry: 100um x 13um-50um
•  Control geometry: 100um x 10-25um
•  Applications: suspension of large particles (eukaryotic
cells, large beads, etc.)

Studer, Vincent, Giao Hang, Anna Pandolfi, Michael Ortiz, W. French Anderson, and Stephen R. Quake. "Scaling properties of a low-actuation pressure
microfluidic valve." Journal of Applied Physics 95, no. 1 (2004): 393-398.
https://sharedfacilities.stanford.edu/service_center/show_external/22/microfluidics-foundry
 Pneumatic Valves

Unger, Marc A., Hou-Pu Chou, Todd Thorsen, Axel Scherer, and Stephen R. Quake. "Monolithic microfabricated valves and pumps by multilayer soft
lithography." Science 288, no. 5463 (2000): 113-116.
 Pneumatic Pumps

Multiple control channels –

alternating pressure in/out.

Peristaltic pumping:
•  Frequency dependent
flow rate

Unger, Marc A., Hou-Pu Chou, Todd Thorsen, Axel Scherer, and Stephen R. Quake. "Monolithic microfabricated valves and pumps by multilayer soft
lithography." Science 288, no. 5463 (2000): 113-116.
 Mechanical Valves

Flexible microfluidic device

with single deformable arch.

Stretching /bending the device

reduces the arch height,
partially opens the channel.

Douglas P. Holmes, Behrouz Tavakol, Guillaume Froehlicher, and Howard A. Stone. "Control and manipulation of microfluidic flow via elastic

6
deformations." Soft Matter 9, no. 29 (2013): 7049-7053.
 Mechanical Valves

Douglas P. Holmes, Behrouz Tavakol, Guillaume Froehlicher, and Howard A. Stone. "Control and manipulation of microfluidic flow via elastic
deformations." Soft Matter 9, no. 29 (2013): 7049-7053.
 Mechanical Valves
Buckled arch in microfluidic chamber.

Arch shape obtained from buckling

analysis of elastica. Apex:

Douglas P. Holmes, Behrouz Tavakol, Guillaume Froehlicher, and Howard A. Stone. "Control and manipulation of microfluidic flow via elastic
deformations." Soft Matter 9, no. 29 (2013): 7049-7053.
 Mechanical Valves
Navier-Stokes Equations:
…conservation of momentum…

Continuity Equation:
…conservation of mass…

Stokes Equations:
Neglect inertia & body forces

Incompressible:

B. Tavakol, G. Froehlicher, D.P. Holmes, and H.A. Stone. “Extended Lubrication Theory: Estimation of Fluid Flow in Channels with Variable Geometry,” Under
Review: Physics of Fluids, (2015).
 Mechanical Valves
Stokes Equations:
Neglect inertia & body forces

Incompressible:

Dimensionless Parameters:
Lengths:

Velocities:

Shape of the impingement:

Pressure:

B. Tavakol, G. Froehlicher, D.P. Holmes, and H.A. Stone. “Extended Lubrication Theory: Estimation of Fluid Flow in Channels with Variable Geometry,” Under
Review: Physics of Fluids, (2015).
 Mechanical Valves
Dimensionless Equations:

Boundary Conditions:
No slip

Total flow rate is prescribed

Perturbation Expansion:
Shape of the impingement:

B. Tavakol, G. Froehlicher, D.P. Holmes, and H.A. Stone. “Extended Lubrication Theory: Estimation of Fluid Flow in Channels with Variable Geometry,” Under
Review: Physics of Fluids, (2015).
 Extended Lubrication
Fluid flow through channels with variable geometry

Extended Lubrication Theory

B. Tavakol, G. Froehlicher, D.P. Holmes, and H.A. Stone. “Extended Lubrication Theory: Estimation of Fluid Flow in Channels with Variable Geometry,” Under
Review: Physics of Fluids, (2015).
 Extended Lubrication

B. Tavakol, G. Froehlicher, D.P. Holmes, and H.A. Stone. “Extended Lubrication Theory: Estimation of Fluid Flow in Channels with Variable Geometry,” Under
Review: Physics of Fluids, (2015).
 Extended Lubrication

B. Tavakol, G. Froehlicher, D.P. Holmes, and H.A. Stone. “Extended Lubrication Theory: Estimation of Fluid Flow in Channels with Variable Geometry,” Under
Review: Physics of Fluids, (2015).
 Bioinspiration

Large, rapid, reversible diaphragm deformation for fluid pumping.

Davis, N. T. and J. G. Hildebrand (2006). "Neuroanatomy of the sucking pump of the moth, Manduca sexta (Sphingidae, Lepidoptera)." Arthropod

7
Structure & Development.
Eberhard, S. H. and H. W. Krenn (2005). "Anatomy of the oral valve in nymphalid butterflies and a functional model for fluid uptake in Lepidoptera."
Zoologischer Anzeiger - A Journal of Comparative Zoology 243(4): 305-312.
 Electrical Valves
Electrodes

V
Dielectric

B. Tavakol, M. Bozlar, G. Froehlicher, H. A. Stone, I. Aksay, and D. P. Holmes, “Buckling Instabilities of Dielectric Elastomeric Plates for Flexible Microfluidic
Pumps,” Soft Matter, 10(27), 4789–4794, (2014).
 Electrical Valves

Coupled Chambers

B. Tavakol, M. Bozlar, G. Froehlicher, H. A. Stone, I. Aksay, and D. P. Holmes, “Buckling Instabilities of Dielectric Elastomeric Plates for Flexible Microfluidic
Pumps,” Soft Matter, 10(27), 4789–4794, (2014).
 Fluid Electrodes

At what voltage will the plate buckle?

Radial strain in a Plane stress relation Linear stability analysis

dielectric elastomer for a buckling plate

B. Tavakol, A. Chawan, and D. P. Holmes, “Buckling Instability of Thin Films as a Means to Control or Enhance Fluid Flow within Microchannels,”
in preparation, (2015).
 Microfluidic Fabrication

1.  Prepare the top substrate with the

microfluidic channel.
2.  Prepare the bottom substrate with all
controlling channels.
3.  Make the thin, dielectric film.
_
4.  Bond the substrates on both sides of
the dielectric film.
•  Film will buckle at the intersection +
between channels.
5.  Fill the channels with conductive
fluids.
6.  Apply a voltage to induce buckling

B. Tavakol, A. Chawan, and D. P. Holmes, “Buckling Instability of Thin Films as a Means to Control or Enhance Fluid Flow within Microchannels,”
in preparation, (2015).
 Valving

B. Tavakol, A. Chawan, and D. P. Holmes, “Buckling Instability of Thin Films as a Means to Control or Enhance Fluid Flow within Microchannels,”
in preparation, (2015).
 Pumping

B. Tavakol, A. Chawan, and D. P. Holmes, “Buckling Instability of Thin Films as a Means to Control or Enhance Fluid Flow within Microchannels,”
in preparation, (2015).
 Overview
Confined Fluid Flow:
Microfluidics and Capillarity
Reynolds Number: Inertia vs. Viscous effects
•  Review of characteristic flows…

Péclet Number: Transport phenomena in a continuum

•  Diffusion, separation, and mixing…

Geometric confinement: Controlling and manipulating fluid flow

•  Microfluidic fabrication, valving, pumping…

Capillary Number: Viscosity vs. Surface tension

•  Droplet formation, capillary rise, elasticity…
 Fluid Behavior
Viscous Forces Interfacial Forces

Coiling honey Wetting of water on a textured surface

Capillary Number: viscous/interfacial

Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
Droplet: http://www.rycobel.be/en/technical-info/articles/1337/measuring-dynamic-absorption-and-wetting
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Monodisperse droplet generation

•  Droplet emulsions in immiscible fluids
•  Injection of water into stream of oil

Interfacial tension prevents the fluids

from flowing alongside each other.

Interfacial Forces
Surface tension acts to reduce the interfacial area.

Viscous stresses act to extend and drag the interface downstream.

Characteristic droplet size: Capillary number:

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
Thorsen, Todd, et al. "Dynamic pattern formation in a vesicle-generating microfluidic device." Physical review letters 86.18 (2001): 4163.
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
Droplet: http://www.rycobel.be/en/technical-info/articles/1337/measuring-dynamic-absorption-and-wetting
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Drop formation in a flow-

focusing configuration

Water drops are formed in

a silicone oil continuous
Interfacial Forces phase.

Qi is the flow rate of water

Qo is the flow rate of oil

Characteristic droplet size:

Stone, Howard A., Abraham D. Stroock, and Armand Ajdari. "Engineering flows in small devices: microfluidics toward a lab-on-a-chip." Annu. Rev. Fluid
Mech. 36 (2004): 381-411.
Anna et al. 2003
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Interfacial Forces

Sequential drop breakup at T-junctions

•  Generally, it is difficult to form small drops at high-volume fraction.
•  This “passive” route allows size reduction after fabrication.
•  The ratio of two daughter droplets depends on the lengths of the arms off the T-
junction.

Stone, Howard A., Abraham D. Stroock, and Armand Ajdari. "Engineering flows in small devices: microfluidics toward a lab-on-a-chip." Annu. Rev. Fluid
Mech. 36 (2004): 381-411.
Link et al. 2003
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Interfacial Forces

Jiandi Wan and Howard A. Stone. "Microfluidic generation of a high volume fraction of bubbles in droplets." Soft Matter 6, no. 19 (2010): 4677-4680.
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
Droplet: http://www.rycobel.be/en/technical-info/articles/1337/measuring-dynamic-absorption-and-wetting
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Large surface-to-volume ratios in microfluidic devices

•  Makes surface effects increasingly important.
•  Important when free fluid surfaces are present.

Surface tensions can exert significant stress

•  Result in free surface deformations.
Interfacial Forces •  Can drive fluid motion.

Capillary forces tend to draw fluid into wetting

microchannels
•  Occurs when solid-liquid interfacial energy is lower than
the solid-gas interfacial energy.

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
Honey: http://www.honeyassociation.com/webimages/honey-dipper.jpg
Droplet: http://www.rycobel.be/en/technical-info/articles/1337/measuring-dynamic-absorption-and-wetting
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Capillary forces tend to draw fluid into wetting

microchannels

Dynamics of surface tension-driven intrusion into a pipe of

radius w occurs from balance of capillary and viscous forces.

Interfacial Forces Curved meniscus at fluid-gas interface – Laplace pressure

Difference in surface energies (solid-liquid to solid-gas)

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Capillary forces tend to draw fluid into wetting

microchannels

Dynamics of surface tension-driven intrusion into a pipe of

radius w occurs from balance of capillary and viscous forces.

Interfacial Forces Stoke’s flow: Scaling of the flow rate:

From the pressure – surface tension relationship:

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Capillary forces tend to draw fluid into wetting

microchannels
From the pressure – surface tension relationship:

Interfacial Forces Capillary number determines the dynamics:

The column length changes as the fluid front moves: ( )

•  Fluid invades the channel at an ever slowing rate:

Washburn
equation
Washburn, 1921
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Surface patterning can create “wall-less”

microchannels which confine fluids.

Chemically treat a surface to make some

areas hydrophilic and others hydrophobic.
Spreading along hydrophilic stripe is similar to the Washburn
Interfacial Forces analysis, however the height of the fluid stripe:

•  Rr is the radius of curvature of the inlet reservoir drop.

•  Spreading requires:
Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
Gau et al., 1999
Darhuber et al., 2001
 Fluid Behavior
Viscous Forces
Capillary Number: viscous/interfacial

Chemically treat walls of microfluidic channel

to make pressure-sensitive valves.
Pressure sensitive pumping
•  From pressure drop across fluid channel (different from
Quake valve)
•  Fluid is confined to hydrophilic region if:

Interfacial Forces

Channel is patterned with hydrophilic stripe (center), weakly hydrophobic

surface (bottom channels), and strongly hydrophobic (top channels).

Zhao et al., 2001

Squires, Todd M., and Stephen R. Quake. "Microfluidics: Fluid physics at the nanoliter scale." Reviews of modern physics 77.3 (2005): 977.
 Capillary Rise
Classical Problem:
•  Noted as early as 15th by Leonardo
da Vinci.

•  Attributed to circulation in plants in

17th century by Montanari.

Balance: Surface Tension & Gravity

J.M. Bell and F.K. Cameron, "The flow of liquids through capillary spaces," J. Phys. Chem. 10, 658-674, (1906).
 Elastocapillarity
Fluid-structure
interaction:

•  Droplet bends and folds

the sheet.

•  Droplet is minimizing the

amount of its surface in
contact with air.

•  Liquid-air surface area is

minimized at the expense
of bending the sheet.

Py, Charlotte, Paul Reverdy, Lionel Doppler, José Bico, Benoit Roman, and Charles N. Baroud. "Capillary origami: spontaneous wrapping of a droplet with
an elastic sheet." Physical Review Letters 98, no. 15 (2007): 156103.
 Elastocapillarity
Fluid-structure
interaction:

Elastic energy of a plate – bending:

Relation between in-plane strain to

out-of-plane bending:

Bending energy:

Py, Charlotte, Paul Reverdy, Lionel Doppler, José Bico, Benoit Roman, and Charles N. Baroud. "Capillary origami: spontaneous wrapping of a droplet with
an elastic sheet." Physical Review Letters 98, no. 15 (2007): 156103.
 Elastocapillarity
Fluid-structure
interaction:
Bending energy:

Surface energy:

Elastocapillary length:

Elastocapillary bending of sheet:

Py, Charlotte, Paul Reverdy, Lionel Doppler, José Bico, Benoit Roman, and Charles N. Baroud. "Capillary origami: spontaneous wrapping of a droplet with
an elastic sheet." Physical Review Letters 98, no. 15 (2007): 156103.
 Elastocapillarity
Capillary rise between flexible fibers.

Curvature of meniscus

Assume gap is much smaller

than width:

Initial gap:

Approximation of
meniscus curvature:
J. Bico et al (2004)

H.-Y. Kim and L. Mahadevan, "Capillary rise between elastic sheets, " J. Fluid Mech. 548, 141-150, (2006).
J.M. Aristoff, C. Duprat, and H.A. Stone, "Elastocapillary Imbibition," Int. J Nonlinear Mech. 48, 648-656, (2011).
C. Duprat, J.M. Aristoff, and H.A. Stone, "Dynamics of elastocapillary rise," J. Fluid Mech. 679, 641-654, (2011).
 Elastocapillarity
Capillary rise between flexible fibers.

Approximation:

Scaling:

J. Bico et al (2004)
Balance: Bending & Surface Tension
H.-Y. Kim and L. Mahadevan, "Capillary rise between elastic sheets, " J. Fluid Mech. 548, 141-150, (2006).
J.M. Aristoff, C. Duprat, and H.A. Stone, "Elastocapillary Imbibition," Int. J Nonlinear Mech. 48, 648-656, (2011).
C. Duprat, J.M. Aristoff, and H.A. Stone, "Dynamics of elastocapillary rise," J. Fluid Mech. 679, 641-654, (2011).
 Elastocapillarity
Capillary rise between flexible fibers.

Potential energies:
Elastic energy

Gravitational potential energy

Surface energy

H.-Y. Kim and L. Mahadevan, "Capillary rise between elastic sheets, " J. Fluid Mech. 548, 141-150, (2006).
J.M. Aristoff, C. Duprat, and H.A. Stone, "Elastocapillary Imbibition," Int. J Nonlinear Mech. 48, 648-656, (2011).
C. Duprat, J.M. Aristoff, and H.A. Stone, "Dynamics of elastocapillary rise," J. Fluid Mech. 679, 641-654, (2011).
 Elastocapillarity
Capillary rise between flexible fibers.

Characteristic length scales:

Elastocapillary length

Capillary gravity length

Elastogravity length

H.-Y. Kim and L. Mahadevan, "Capillary rise between elastic sheets, " J. Fluid Mech. 548, 141-150, (2006).
J.M. Aristoff, C. Duprat, and H.A. Stone, "Elastocapillary Imbibition," Int. J Nonlinear Mech. 48, 648-656, (2011).
C. Duprat, J.M. Aristoff, and H.A. Stone, "Dynamics of elastocapillary rise," J. Fluid Mech. 679, 641-654, (2011).
 Elastocapillarity
Capillary rise between flexible fibers.

Dimensionless parameters
Bond number:

Elastocapillary number

Elastogravity number

H.-Y. Kim and L. Mahadevan, "Capillary rise between elastic sheets, " J. Fluid Mech. 548, 141-150, (2006).
J.M. Aristoff, C. Duprat, and H.A. Stone, "Elastocapillary Imbibition," Int. J Nonlinear Mech. 48, 648-656, (2011).
C. Duprat, J.M. Aristoff, and H.A. Stone, "Dynamics of elastocapillary rise," J. Fluid Mech. 679, 641-654, (2011).
 Elastocapillarity

H.-Y. Kim and L. Mahadevan, "Capillary rise between elastic sheets, " J. Fluid Mech. 548, 141-150, (2006).
J.M. Aristoff, C. Duprat, and H.A. Stone, "Elastocapillary Imbibition," Int. J Nonlinear Mech. 48, 648-656, (2011).
C. Duprat, J.M. Aristoff, and H.A. Stone, "Dynamics of elastocapillary rise," J. Fluid Mech. 679, 641-654, (2011).
 Elastocapillarity

H.-Y. Kim and L. Mahadevan, "Capillary rise between elastic sheets, " J. Fluid Mech. 548, 141-150, (2006).
J.M. Aristoff, C. Duprat, and H.A. Stone, "Elastocapillary Imbibition," Int. J Nonlinear Mech. 48, 648-656, (2011).
C. Duprat, J.M. Aristoff, and H.A. Stone, "Dynamics of elastocapillary rise," J. Fluid Mech. 679, 641-654, (2011).
Elastocapillarity
Elastocapillarity
 Capillarity & Swelling

D.P. Holmes, A. Pandey, P.-T. Brun, and S. Protière, In Preparation, (2015).

 Capillarity & Swelling

D.P. Holmes, A. Pandey, P.-T. Brun, and S. Protière, In Preparation, (2015).

 Capillarity & Swelling

1. Elastocapillary rise between flexible fibers.

D.P. Holmes, A. Pandey, P.-T. Brun, and S. Protière, In Preparation, (2015).

 Capillarity & Swelling

1. Elastocapillary rise between flexible fibers.

Stationary meniscus height rises linearly with

elastocapillary length.
D.P. Holmes, A. Pandey, P.-T. Brun, and S. Protière, In Preparation, (2015).
 Polymers & Swelling

Swelling
Bad Solvent Good Solvent
 Capillarity & Swelling

2. Swelling-induced bending.

Dynamics of bending ~ Diffusion of temperature

D.P. Holmes, A. Pandey, and S. Protière, In Preparation, (2015).

D.P. Holmes, M. Roché, T. Sinha, and H.A. Stone. “Bending and Twisting of Soft Materials by Non-homogenous Swelling” Soft Matter, 7, 5188, 2011.
 Capillarity & Swelling

2. Swelling-induced bending.

•  Thermal diffusion through the beam thickness.

•  Shape obtained by minimizing the bending
moment in the beam.
•  Beam curvature as temperature diffuses.

•  Poroelastic time scale μ = Solvent viscosity

h = Thickness
k = Permeability (k ≈ 10-18 m2/s)
E = Elastic modulus (E = 106 Pa)
D.P. Holmes, A. Pandey, and S. Protière, In Preparation, (2015).

D.P. Holmes, M. Roché, T. Sinha, and H.A. Stone. “Bending and Twisting of Soft Materials by Non-homogenous Swelling” Soft Matter, 7, 5188, 2011.
 Capillarity & Swelling

2. Swelling-induced bending.

D.P. Holmes, A. Pandey, and S. Protière, In Preparation, (2015).

D.P. Holmes, M. Roché, T. Sinha, and H.A. Stone. “Bending and Twisting of Soft Materials by Non-homogenous Swelling” Soft Matter, 7, 5188, 2011.
 Capillarity & Swelling

3. Bending dominates surface tension.

Bending energy

Surface energy
Peeling occurs if the curvature
exceeds the bending capillary length.
D.P. Holmes, A. Pandey, and S. Protière, In Preparation, (2015).

Roman, Benoit, and José Bico. "Elasto-capillarity: deforming an elastic structure with a liquid droplet." Journal of Physics: Condensed Matter 22.49
(2010): 493101.
 Capillarity & Swelling

1. Elastocapillary rise between flexible fibers.

At short times, elastocapillary rise dominates the

deformation.

2. Swelling-induced bending.

Bending is constrained by surface tension, as the beam

bends with a lower curvature than a free swelling beam.

3. Bending dominates surface tension.

Separation occurs as the “natural” curvature of the

beam exceeds the fluids ability to confine it.

D.P. Holmes, A. Pandey, and S. Protière, In Preparation, (2015).

 Capillarity & Swelling

Stretching energy Bending energy Surface energy

D.P. Holmes, A. Pandey, and S. Protière, In Preparation, (2014).

 Capillarity & Swelling

D.P. Holmes, A. Pandey, and S. Protière, In Preparation, (2014).

 Baobab Flowering

BBC - Planet Earth: Seasonal Forests

 Summary
Confined Fluid Flow:
Microfluidics and Capillarity
Reynolds Number: Inertia vs. Viscous effects
•  Review of characteristic flows…

Péclet Number: Transport phenomena in a continuum

•  Diffusion, separation, and mixing…

Geometric confinement: Controlling and manipulating fluid flow

•  Microfluidic fabrication, valving, pumping…

Capillary Number: Viscosity vs. Surface tension

•  Droplet formation, capillary rise, elasticity…