This invention relates to methods and apparatus for microfluidic mixing and refers particularly, though not exclusively, to such methods and apparatus based on instability caused by viscoelastic behavior of fluids.
Throughout this specification a reference to microfluidic mixing is to be taken as including mixing at micro-length scale as well as smaller length scales and larger length scales.
Microfluidic devices and methods have enabled technologies for analytical chemistry and biochemical analysis. In general, analysis processes are carried out at the micro length scale. Because of the larger surface-to-volume ratio, flow in microstructures is laminar and stable. Since mixing is a key process for all chemical processes or most microfluidic applications, effective and fast mixing under laminar conditions is required. A number of micromixer designs have been proposed.
Generally, micromixers can be categorized as being either passive or active. Active micromixers require actuators and involve moving parts. As such, they are not attractive for disposable applications. Passive micromixers have no actuators and no moving parts. Passive mixing concepts rely on molecular diffusion or chaotic advection. Passive mixers based on molecular diffusion utilize concepts such as parallel lamination, serial lamination, and serial segmentation to reduce mixing time and to shorten mixing paths in microchannels.
Passive mixers based on chaotic advection use instabilities caused by geometrical modifications at medium and high Peclet numbers: Pe=
For viscous fluid flow, inertial and viscous forces are relevant and are typified by the Reynolds number: Re=ρ
It is well known that the elasticity of a viscoelastic fluid can introduce elastic stress in addition to viscous stress. The stress experienced by a viscoelastic fluid will not immediately become zero with the cessation of fluid motion and driving forces, but will decay with a characteristic time due to its elasticity. An example of a viscoelastic fluid is a fluid with dilute (i.e. a minute amount of) deformable and high molecular weight polymers. Viscoelastic instability of these non-Newtonian fluids is known.
In general, a viscoelastic fluid with larger and/or higher concentrations of polymer molecules has a longer relaxation time, while a smaller channel has a shorter flow characteristic time. For a viscoelastic fluid flow in a given micro geometry, as the dimensions of a channel decrease, the Re becomes smaller and it is more difficult to have inertia/viscous flow instability; but the Deborah number (De=characteristic relaxation time/flow characteristics time or elastic forces/viscous forces) becomes larger and it is easier to have elastic/viscous instability. Hence, for viscoelastic fluid flows in microchannels, the inertial effects are negligible, and the flow is dominated by viscoelastic forces.
For microfluidics, it is common to have a Reynolds number in the order of unity and a laminar flow is expected for Newtonian fluids and fluids with negligible elasticity.
Microfluidic devices are the key to micro-scale analytical chemistry and biochemical analysis. With large surface to volume ratio and small characteristic length, flow field in microchannels is normally laminar and stable. Without employing viscoelastic fluids, the mixing of two or more streams is normally only able to be achieved by diffusion, and not by the more effective mechanism of flow instability and/or turbulence. However, diffusive mixing will compromise the requirements of short mixing path and time for efficient mixing.
In accordance with a first preferred aspect there is provided apparatus for microfluidic mixing that has at least one first fluid inlet for a first fluid, the at least one first fluid inlet being operatively connected to a first fluid channel. At least one second fluid inlet is provided for at least one second fluid operatively, the at least one second fluid inlet being operatively connected to a second fluid channel. The at least one second fluid channel operatively intersects the first fluid channel for introduction of the second fluid into the first fluid channel. The first fluid channel has at least one outlet end remote from the at least one first fluid inlet. There is at least one contraction in the first fluid channel in at least one location selected from: intermediate the intersection of the first fluid channel with the at least one second fluid channel and the outlet end, and intermediate the at least one first fluid inlet and the intersection of the first fluid channel with the at least one second fluid channel.
According to a second preferred aspect there is provided a method for microfluidic mixing. The method comprises supplying a first fluid to at least one first fluid inlet for flow along a first fluid channel, the first fluid channel having at least one outlet end remote from the at least one first fluid inlet and at least one contraction in the first fluid channel in at least one location selected from: intermediate the intersection with at least one second fluid channel and the outlet end, and intermediate the at least one first fluid inlet and the intersection with the at least one second fluid channel. At least one second fluid is supplied to at least one second fluid inlet for flow along the at least one second fluid channel, the at least one second fluid channel operatively intersecting the first fluid channel for introduction of the at least one second fluid into the first fluid channel for a first stage of mixing of the first fluid and the at least one second fluid. The first fluid and the at least one second fluid are then passed through the contraction for a second stage of mixing of the first fluid and the at least one second fluid.
For both aspects the at least one second fluid channel may comprise two channels, one on each side of the first fluid channel. The two channels may be identical, or non-identical. The two channels may intersect the first fluid channel at the same or different locations along the first fluid channel. The at least one contraction may be an abrupt contraction/expansion and may have a ratio of x:y:z, with x and z greater than y. The ratio may be determined by a flow rate of the mixed first fluid and at the least one second fluid, a viscosity of the mixed first fluid and at the least one second fluid, an elasticity of the mixed first fluid and at the least one second fluid, and to reduce dead volume. The first fluid inlet, first fluid channel, second fluid inlet, second fluid channel, and the outlet end may be formed in an upper portion of a substrate. A lower portion of the substrate may close the first fluid inlet, first fluid channel, second fluid inlet, second fluid channel and the outlet end. The first channel may comprise an upstream portion upstream of the contraction, and a downstream portion downstream of the contraction; the first stage mixing being by viscoelastic instability taking place in the upstream portion, and the second stage mixing being by viscoelastic flow instability and expansive flow taking place in the downstream portion.
In order that the present invention may be fully understood and readily put into practical effect, there shall now be described by way of non-limitative example only preferred embodiments of the present invention, the description being with reference to the accompanying illustrative drawings.
In the drawings:
Each second channel 20 may have its own inlet well 18 so that the fluids can be different. In this way there would be three fluid inlets and channels, for three different fluids. The fluids may be input at different times, and at different flow rates.
The first fluid channel 16 has at least one outlet well 24 at its end. The outlet well 24 may be centered in the substrate 12 or may be located on the periphery of the substrate 12. The outlet well 24 may connect to a second fluid handling device (not shown), and the second fluid handling device may be a duplication of the apparatus 10, or may be a different device. The bottom layer 32 may be transparent to allow optical access. The channels 16, 20 may have a cross-section shape that is rectangular, circular, oval, or trapezoidal, or otherwise as required or desired. The first channel 16 may be larger in cross-sectional area than each of the second channels 20. Alternatively, the first channel 16 may be similar in cross-sectional area as each of the second channels 20. Further alternatively, the first channel 16 may be smaller in cross-sectional area than each of the second channels 20.
The substrate 12 has a top layer 30 and a bottom layer 32 which is generally parallel to the top layer 30. The channels 16, 20 and wells 14, 18, 24 are fabricated onto the top layer 30, and are sealed by the bottom layer 32, which is generally flat. The wells 14, 18, 24 and part of the channels 16, 20 may be located on the bottom layer 32 if desired or required. The substrate 12 may be of any suitable material such as, for example, polymer, silicon, metal, glass, ceramic, or any combination of them.
Intermediate the outlet well 24 and intersection 26 is a contraction 28. The contraction 28 is preferably an abrupt contraction/expansion which may have a ratio of x:y:z, with x and z greater than y. More preferably is at a ratio of at least 4:1:4. Alternatively or additionally, the contraction 28 in the first fluid channel 16 may be intermediate the first fluid inlet 14 and the intersection 26. There may be more than one contraction. If there are more than one, they may be the same or different. They may be located intermediate the intersection 26 and the outlet well 24 and intermediate the first fluid inlet 14 and the intersection 26.
As shown in
The contraction/expansion ratio, a ratio of x:y:z, with x and z greater than y, is determined by:
The contraction ratio (x:y:z) should have x and z greater than y. The ratio used for
The chaotic velocity profile across the section of the channel 16 in both the upstream portion 34 and downstream portion 36 arises from the combination of viscous forces and elastic forces. The viscoelastic forces give rise to the secondary corner vortices and viscoelastic whipping in the upstream portion 34, and result in a first stage mixing. This irregular flow pattern causes flow fluctuation of the main stream through the contraction 28, and fluctuation of flow resistance to the two side streams in second channel 20. The whipping of the main stream facilitates the side streams penetrating deeply into the central flow in the downstream portion 36. These fluctuations result in viscoelastic flow instability downstream of the contraction 38. This viscoelastically induced flow instability together with the sudden expansive flow in the downstream portion 36, promotes effective and efficient mixing.
In the examples of
The mixing was of two dissimilar fluids. The main stream in first channel 16 was 1 wt % polyethyleneoxide (PEO) in 55 wt % glycerol water (1% PGW for brevity). This has a high viscosity and elasticity. The side streams were 0.1 wt % PEO in water (0.1% PW for brevity). They entered the main microchannel 16 through the two side channels 20. The molecular weight (Mw) of PEO employed was approximately 2×106 g/mol. For image acquisition of the flow fields, a fluorescent dye (fluorescein disodium salt C20H10Na2O5) was added to 1% PGW at a weight ratio of 4×10−4:1 to identify the main stream. For the side streams, 3 μm red fluorescent microsphere solution (Duke Scientific Co.) was added to 0.1% PW at a volume ratio of 0.03:1. The addition of fluorescent dye and microspheres has negligible effects on the fluid properties, and the fluid properties were determined with the additives. For each flow rate, flow field images in the same experiment identified by green fluorescent dye (main stream) or red fluorescent microsphere solution (side stream) were captured at different time by changing the filtering lens.
The total volumetric flowrate is {dot over (Q)}, with the mainstream flowrate being 0.5 {dot over (Q)}, and each of the side streams flowrate being 0.25{dot over (Q)}. The {dot over (Q)} investigated were 10, 20 and 40 ml/hr. The sample fluids were primed into the microchannels by driving the syringes, with appropriate size ratio, using the same micro-syringe pump.
Table 1 contains the rheological properties of the fluids. Relaxation time (λ) were measured from the frequency oscillation test and the steady shear viscosities were determined using amplitude sweep test at shear rates 0.01≦{dot over (γ)}≦1000 s−1. The De, Re and Pe are estimated as: De=λ{dot over (γ)}char where the characteristic shear rate is {dot over (γ)}char=
The diffusion coefficient in dye/water solution was determined as D=1.5×10−9 m2/s. Since the diffusion coefficient is inversely proportional to viscosity, D=3.55×10−12 m2/s approximately for 1% PGW. As such, 1% PGW at {dot over (Q)}=40 ml/hr, the dimensionless parameters are Pe=125×106, Re=0.149, De=139 approximately.
Downstream of the contraction, subsequent to the main and the side streams competing and gushing through the contraction at high speed, the viscoelastic flow instability was more significant. At {dot over (Q)}=10 ml/hr, there was some penetration of the side stream fluid into the central portion of the channel, overlapping with the main stream fluid, see
To study the necessity of fluid elasticity for flow instability and mixing, viscous fluids devoid of elastic effects were employed, namely 96% and 26% glycerol/water solutions were employed in the main and side streams respectively, The viscosity ratio is approximately the same as for the viscoelastic flows, see Table 1. The entire flow field was stable for all the flowrates investigated, for example for {dot over (Q)}=40 ml/hr (Re=0.16, De≈0), see
In
Whilst there has been described in the foregoing description preferred embodiments of the present invention, it will be understood by those skilled in the technology concerned that many variations or modifications in details of design or construction may be made without departing from the present invention.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/SG06/00130 | 5/19/2006 | WO | 00 | 6/25/2008 |
Number | Date | Country | |
---|---|---|---|
60701078 | Jul 2005 | US |