PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Lab Chip. Author manuscript; available in PMC 2010 June 21.
Published in final edited form as:
PMCID: PMC2827298
NIHMSID: NIHMS176015

Evaporation from microreservoirs

Abstract

As a result of very large surface area to volume ratios, evaporation is of significant importance when dealing with lab-on-a-chip devices that possess open air/liquid interfaces. For devices utilizing a reservoir as a fluid delivery method to a microfluidic network, excessive evaporation can quickly lead to reservoir dry out and overall device failure. Predicting the rates of evaporation from these reservoirs is difficult because the position of the air/liquid interface changes with time as the volume of liquid in the reservoir decreases. Here we present a two-step method to accurately predict the rates of evaporation of such an interface over time. First, a simple method is proposed to determine the shape of an air/liquid meniscus in a reservoir given a specific liquid volume. Second, computational fluid dynamics simulations are used to calculate the instantaneous rate of evaporation for that meniscus shape. It is shown that the rate of evaporation is strongly dependent on the overall geometry of the system, enhanced in expanding reservoirs while suppressed in contracting reservoirs, where the geometry can be easily controlled with simple experimental methods. Using no adjustable parameters, the model accurately predicts the position of the inner moving contact line as a function of time following meniscus rupture in poly(dimethylsiloxane) reservoirs, and predicts the overall time for the persistence of liquid in those reservoirs to within 0.5 minutes. The methods in this study can be used to design holding reservoirs for lab-on-a-chip devices that involve no external control of evaporation, such that evaporation rates can be adjusted as necessary by modification of the reservoir geometry.

Introduction

In the past decade, there has been an increasing amount of research dedicated to biological lab-on-a-chip (LOC) applications based on microfluidic unit operations.1 The combined effect of these processes, including volume sampling and pretreatment, fluid pumping, mixing, reaction, separation, and chemical analysis, allow for the completion of biological assays in the span of seconds to minutes.2 The rapid format combined with the small footprint of LOC devices facilitates development of accessible point-of-care diagnostics,3 as well as the simplification of traditional laboratory methods such as polymerase chain reaction (PCR),4 microarray technologies,5 DNA sequencing,6 and chemical analysis.7

Because of the high surface area to volume ratios in micro-fluidic systems, loss of fluid through evaporation and permeation becomes increasingly important.8 Evaporation may rapidly increase the local concentration of solutes in systems that include an open air/liquid interface, which can lead to detrimental effects in an analytical system. If not controlled, evaporation effects in systems containing sub-μL liquid volumes is of critical concern, as times required for the complete evaporation of a liquid sample can be less than a minute. In all microfluidic systems containing an open air/liquid interface, local evaporation rates and their dependence on local conditions and reservoir geometries must be taken into consideration.911

In many open microfluidic systems, attempts are made to either mitigate or enhance evaporation effects. Reducing the evaporation rate is generally required to maintain analyte concentrations when performing chemical analysis on aqueous fluids introduced or extracted in microreservoirs. If no action is taken to replenish the reservoir solvents, sub-μL sized volumes can dry out within a few minutes,12 requiring a sufficiently fast analytical procedure.13 Perhaps the easiest way to control evaporation from these reservoirs is via isolation of the air/liquid interface with either a solid lid14 or an immiscible liquid cap, such as mineral oil15 or another organic liquid.16,17 Other direct approaches involve control of the relative humidity directly above the air/liquid interface, usually performed via sample placement near a water bath to keep the surrounding air saturated with water vapor.9,18 More complicated approaches involve the fabrication of microreservoir arrays from agarose gel.19 These gel-based arrays consist of over 90% water, saturating the local environment above each reservoir and substantially reducing evaporation rates.

Evaporation enhancement in lab-on-a-chip applications is generally associated with passive microfluidic pumping applications. Based on evaporation effects or capillary forces, passive microfluidic pumping eliminates the need for moving parts or external power supplies and is a much sought after topic. Evaporation occurring at the outlet region of a microfluidic network (μFN) serves to drive fluid from an inlet reservoir through the system, usually at substantial volumetric flow rates. A simple and efficient approach to evaporation-based passive pumps can also be realized by exposing the outlet regions of a μFN to the atmosphere. Evaporation in these systems can then be enhanced by either forced convective air motion,20 heating the outlet regions using microfabricated elements,21 or by increasing the air/liquid interface surface area at the μFN outlet with hydrophilic membranes22 or filter paper plugs.23

There are many LOC designs where uncontrolled evaporation can lead to unwanted changes to the initial experimental conditions and, if given enough time, failure of the device. In the case of passive pumping methods based on capillary forces, either between two droplets of water24,25 or liquid in an open faced micro-channel,26 evaporation will continuously reduce the volume of water present in the system. This loss may lead to changes in the capillary forces driving liquid through the system, thus acting to increase (or decrease) the overall observed flow rate. In the case of μNFs used for the immobilization of immunoglobulins to a surface,27 excessive evaporation from the inlet reservoir can lead to reservoir dry out and non-homogeneous deposition of material on to a surface. This behavior may affect the overall reliability of devices that depend on uniform surface coverage of capture ligands, such as those used in various micromosaic immunoassays.28,29

In many LOC devices it becomes important to calculate not only the dependence of device performance on the evaporation rate, but also the time at which that device might undergo failure. The rate of evaporation from an air/liquid interface (Qe) is dependent on both the conditions of the environment above the meniscus as well as the geometry of the meniscus itself. Historically, the study of the evaporation rate and shape of a liquid meniscus from μL-sized volumes has focused on the role of evaporation of a sessile drop for deposition of solids onto a flat surface.3033 While the evaporation of sessile drops remains important for many LOC applications, there is little information regarding evaporation of a droplet of water confined to a micro-reservoir constructed from hydrophilic materials. Previously, this topic has been restricted to experimental studies concerning evaporation of droplets in shallow nanoliter wells, with the focus on solute deposition from the droplets,34 evolution of the meniscus with time,35,36 evaporation rate of the droplets,37 or particle tracking studies of flow induced via evaporation.38 In those studies meniscus shape evolution and the transient behavior of the droplet evaporation rate were not estimated or predicted. The absence of these data may very well be due to the complexity of the dynamics: as water in a microreservoir recedes to the point where the bottom of the meniscus reaches the reservoir floor, the meniscus will rupture due to the discontinuity of the surface energies between the initial air/liquid interface and the eventual solid/liquid interface along the reservoir floor (Fig. 1A). In the small geometry of the reservoir, it may be assumed that body forces have negligible impact on the interface shape such that, prior to rupture the meniscus is spherical. However, after rupture the meniscus takes on a more complicated shape.

Fig. 1
(A) As a meniscus recedes towards the bottom of a reservoir, it will rupture into an incomplete state. (B) Geometric parameters describing a meniscus in a microreservoir in the r,z-plane.

In this study we examine the evaporation rate and shape evolution of an air/water meniscus situated in a cylindrical microreservoir modeled after those used in many LOC devices,2729 specifically, with a reservoir height (H), floor diameter (D1), and upper diameter (D2) all on the order of a millimeter. A simple geometric model is proposed to determine the shape of a meniscus representing a given volume of water in a reservoir having a circular, but not necessarily constant, cross-section. With the interface geometry known, the evaporation rate of the water in the well can then be determined utilizing computational fluid dynamics (CFD). The combined geometric/CFD model is shown to reliably predict the evolution of air/water menisci in experimental systems. It is demonstrated that evaporation of water from microreservoirs can be either enhanced or suppressed, depending on the geometry of the reservoir. The methods utilized within this study will prove useful in a wide array of applications such as the secondary measurement of the surface properties of a cylindrical shaped substrate, the concentration of solutes within a microreservoir, and the controlled deposition of solids onto the reservoir floor. Furthermore, knowledge of a meniscus geometry evolution over time can be utilized in the design of a passive microfluidic pumping device, where the capillary forces resulting from the curved meniscus would serve as a controlled driving force between two microreservoirs.

Numerical techniques

Evolution of a meniscus in a microreservoir

The geometries of the menisci and reservoirs used in this study are shown in Fig. 1(B). To predict the rate of evaporation from a microreservoir containing a given volume of liquid, one needs to know the shape of the air/liquid meniscus. Because the shape of the meniscus is dependent on the surface properties of the reservoir, a relationship is needed between those properties and the geometric parameters of the system. For systems in quasi-thermodynamic equilibrium, the relationship between the properties of the surface and the air/liquid interface can be described by Young’s equation,

γsv=γlvcos(θ)+γsl
(1)

where θ is the apparent contact angle between the liquid andsurface and γsv, γlv, and γsl represent the surface and interfacial free energies of the solid–vapor, liquid–vapor, and solid–liquid interfaces, respectively. For reasonably dilute aqueous solutions θ can be assumed constant, and changes in θ will represent the relative hydrophilic/hydrophobic extent of the surface. In this study, θ1 is defined as the contact angle between water and the reservoir floor and θ2 as the contact angle between water and the reservoir sidewalls. The shape of the meniscus is then dependent on the geometry of the reservoir and surface conditions of the reservoir sidewalls and floor (D1, D2, H, θ1, θ2) as well as the given volume of water within the reservoir (V).

Determination of the shape of a free interface in mechanical equilibrium with a solid edge is not a straightforward process. Analytical solutions may only be found for a small number of cases, and numerical methods are usually required.39 Furthermore, experimental measurements concerning the shape of menisci in microreservoirs via confocal37 or interference-contrast microscopy36 suffer from poor resolution near the reservoir walls. The problem can be simplified, however, by noting that the reservoirs will usually have a floor diameter of D1 < 1 mm, such that for water with density ρ the Bond number can be calculated as Bo = ρgD12/4γlv ≤ 0.05; hence, surface tension forces will dominate gravitational forces. In this case it may be assumed that the two principal radii of curvature, R1 and R2, do not vary with respect to position along the interface. In Fig. 1(B), R1 is defined as the radius of curvature of the meniscus in the r,z-plane, having a center of curvature based at the point (ro, zo). R2 is defined (for incomplete menisci) as the radius of curvature in the r, θ-plane, extending from r = 0 to the average radial position of the meniscus. For a given reservoir liquid volume, the shape of the meniscus is then calculated from a curve that intersects the reservoir floor and sidewalls with intersection angles θ1 and θ2, respectively, and whose volume when integrated with respect to r, θ, and z is equivalent to the given volume of liquid in the reservoir. Because there is no explicit expression for ro, zo, and R1 given H, D1, D2, θ1, θ2, and V, an iterative method must be utilized. The details of this process are given in Appendix I in the ESI.

Calculation of evaporation rates

An accurate relationship between the instantaneous evaporation rate Qe and the liquid volume V(t) is needed to accurately predict the evolution of a meniscus over time. From Fig. 1, it may be seen that for a given reservoir geometry, Qe will depend on the volume of water in the reservoir and cannot be assumed constant as in previous studies.3638 This transient behavior is primarily due to the increase in the average distance needed for water vapor diffusion as well as changes to the interfacial area as the meniscus recedes towards the reservoir floor. To proceed further with the development of the model it is necessary to make several additional, but physically justifiable assumptions.

The microreservoirs used in this study have an average diameter and height on the order of a millimeter. For these geometries it is reasonable to assume the gas (air) and liquid in the reservoir are essentially quiescient, and that water vapor is transported upward via diffusion. The mole fraction of water vapor just above the well (xa2) is assumed constant, determined by the relative humidity (RH), pressure (p), and absolute temperature (T) of the system. The gas phase at the air/water meniscus reaches near-instantaneous equilibrium, and the mole fraction of water vapor at the interface (xa1) may be determined using Raoult’s law. Because the meniscus recedes very slowly as evaporation proceeds the system is at pseudo steady-state, and also it can further be assumed that, due to the high surface-area to volume ratio, the system will be isothermal.

For dilute concentrations of water vapor at constant pressure, the volumetric evaporation rate from the meniscus can then be expressed as

Qe=pMwDρRTAxazdA
(2)

where A is the cross-sectional area of the reservoir in the r,θ-plane, R is the gas constant, Mw is the molecular weight, ρ is the liquid density, and An external file that holds a picture, illustration, etc.
Object name is nihms176015ig1.jpg is the diffusion coefficient of water vapor in air. Because of the curved geometry of the meniscus, exact solutions for xa = xa(r, z) in the gas phase cannot be derived and numerical methods must be applied to determine the evaporation rate from reservoirs with given parameters H, D1, D2, θ1, θ2, and V. Integration of eqn (2) is performed along a plane of constant z, where conservation of mass requires that the plane of integration must be situated at a z-position located higher than the top edge of the meniscus.

Fig. 2 illustrates the domain modeled in this study. Numerical calculations to determine the water vapor mole fraction distribution are performed using a 2-D axisymmetric CFD model. For each reservoir with given surface properties and external conditions a minimum of 25 simulations were performed, each pertaining to a specific reservoir liquid volume, where the shape of the meniscus was predetermined using the geometric model described above. The CFD preprocessing package Gambit® was used to discretize the gas phase within the reservoir into 4-node elements, with typical element dimensions on the order of D1/100. This discretization scheme has been pre-determined in this system to have sufficient resolution to ensure the predicted concentration field is independent of the mesh density. The finite volume CFD package FLUENT® was then used to simulate the 2-D concentration distribution within the r,z-plane, where the mole fractions of water vapor along the meniscus and at the reservoir outlet were set to xa1 and xa2, respectively. The species conservation equation for water vapor was solved using a second-order accurate scheme, and the simulations were considered converged when the normalized residuals fell below 10−7. After the solution for xa has converged the evaporation rate is calculated using eqn (2), where the integration is performed along the reservoir outlet plane at z = H.

Fig. 2
Illustration of the axi-symmetric CFD model utilized in this study. The computational grid is for presentation purposes only.

Experimental techniques

The sidewalls of the microreservoirs were constructed from poly(dimethylsiloxane) (PDMS, Sylgard 184, Dow Corning), whereas the floors consisted of either PDMS or glass. After curing a flat PDMS mold (10 : 1 base : curing agent) on a polished silicon wafer at 95 °C for a minimum of 3 hours, high precision biopsy punches (Technical Innovations) were used to create the reservoir chamber, making sure the reservoir was as close to vertical as possible. The biopsy punch creates a conical cavity, which allows for an expanding (D2 > D1) or contracting (D2 < D1) reservoir, depending on which side of the substrate the reservoir floor is sealed. Prior to sealing, all PDMS substrates were treated via an extraction, which has been shown to stabilize the hydrophilic nature of an oxidized PDMS surface for up to 7 days.40 The PDMS reservoir surfaces were rendered hydrophilic by exposure to air plasma (25 W, 900 mTorr, 20 s), after which the reservoirs were irreversibly sealed to either a glass or PDMS floor.

The reservoirs were typically used within 10 minutes after fabrication. A microliter pipette was used to inject solutions of 0.5 μm diameter rhodamine encapsulated polystyrene beads (Duke Scientific, 106 beads/mL) in nanopure water into the reservoirs. Before filling, the bead solutions were agitated to ensure uniform distribution of the microspheres. Initial liquid volumes (Vo) ranged from 0.2 to 0.75 μL, depending on the reservoir geometry under consideration. The temperature, pressure, and relative humidity of the external environment was recorded at the onset of the experiment (Fluke 971 temperature humidity meter), and typically remained within 1 °C and 1.5% RH of the initial conditions for several hours. The diameters at the top and bottom of the reservoir were measured using a calibrated microscope, and the height of the reservoir was measured using a micrometer.

After filling (time t = 0), the reservoir was placed on an inverted microscope equipped with an epifluorescence attachment (Nikon TE-200U) and fitted with a CCD camera (Coolsnap fx). A series of images taken immediately after menisci rupture were then used to measure the diameter of the inner meniscus (Dcap) as a function of time. Figs. 3(A) and (B) show the diagram of the experimental setup and images displaying the meniscus before and after rupture. Fig. 3(C) displays Dcap as a function of time for an initial liquid volume of Vo = 0.65 μL in a reservoir with D1 = 1.38 mm, D2 = 0.77 mm, and H = 2.2 mm. The rupture event occurred at tc = 42.9 min with an initial rupture diameter of Dcap(tc) = 0.47 mm, and the reservoir was completely dry at tend = 64.7 min.

Fig. 3
(A) Illustration of the experimental techniques in this study. (B) Experimental images of a ruptured meniscus for a reservoir at T = 28 °C, RH = 31%, with geometry D1 = 1.38 mm, D2 = 0.77 mm, and H = 2.2 mm.

Discussion

Effect of reservoir geometries and surface properties on evaporation rate

Experimentally observed evaporation rates of liquids contained in shallow reservoirs (D1/H [dbl greater-than sign] 1) have been shown to be proportional to the diameter of the reservoir and independent of time.3638 These functional relationships do not hold true for taller reservoirs such as those in this study (D1/H ≈ 1). Fig. 4 displays mole fraction distributions and evaporation rates for water vapor within several reservoir geometries and reservoir liquid volumes. The difference in geometry between the expanding and contracting reservoirs would result from a reversed punch direction regarding a PDMS substrate. Several geometrical relationships affecting Qe can be seen in Fig. 4. In general, a decrease in Qe is seen with decreasing liquid volumes due to the increase in the average distance needed for diffusion as the meniscus recedes towards the reservoir floor. However, the evaporation rate is not only dependent on this diffusion distance, but also on the overall interfacial area as well as the geometry of the reservoir itself. The difference in evaporation rates for different reservoir geometries can be quite dramatic; for example, the evaporation rate for an expanding reservoir (for the conditions shown in Fig. 4) with V = 0.1 μL (Qe = 0.36 pL s−1) is only slightly lower than the case for a contracting reservoir with V = 1 μL (Qe = 0.42 pL s−1), even though there is clearly a difference in the interfacial area and diffusion length between the two cases. For a given volume of liquid, Qe will be higher for an expanding reservoir than a contracting reservoir for several reasons: (1) the position of the meniscus will be higher in an expanding reservoir at a given volume, reducing the distance needed for diffusion and (2) evaporation rates are enhanced in expanding reservoirs due to the increase in the available area for diffusion area as z increases. The latter effect is made clear by examining one-dimensional diffusion in a well, and is examined in more detail in Appendix II, available within the ESI.

Fig. 4
CFD simulations displaying the mole fraction of water vapor in several reservoirs with differing volumes. The simulations pertain to a reservoir with θ1 = 25°, θ2 = 50°, H = 2 mm, T = 28 °C, p = 686 Torr, and RH ...

To determine a functional relationship between Qe and V for a given reservoir with overall volume Vres, a minimum of 25 CFD simulations such as those shown in Fig. 4 were performed, with each simulation pertaining to a different reservoir liquid volume. The simulations are carried out such that a minimum of 10 simulations are evenly spaced along the range of Vc < V < Vo, corresponding to a complete meniscus, and a minimum of 15 simulations are evenly spaced along the range of −3 < log(V) < log(Vc) for an incomplete meniscus, where Vc is the critical volume at which the meniscus ruptures. Interpolation between the resulting data points is accomplished by fitting Qe vs. V to the following relationships,

completemeniscus(V>Vc)Qe=a1V2+a2V+a3
(3)

incompletemeniscus(V>Vc)Qe=b1log(V)+b2log(VVres)b3
(4)

where the two distinct relationships are needed due to the discontinuous nature of the meniscus before and after rupture. The regression analysis within this study yielded R2 coefficient values greater than 0.95 for systems that exhibit stable menisci (to be discussed later). There is little physical meaning associated with the fitting parameters in eqns. (3) and (4) and the parameters describing the system, and as such their values are not reported here.

Fig. 5(a) displays Qe vs. V for a variety of reservoir diameters and heights for a set value of D1 = 1.25 mm and contact angles of θ1 = 25° and θ2 = 50°, representative of a typical reservoir in this study constructed from PDMS. Although the two curves for complete and incomplete menisci are discontinuous at V = Vc, the difference in Qe between the two states at V = Vc is predicted to be small. For a given reservoir, the decrease in Qe with decreasing V is clear in all reservoirs; for example, for the case of D2 = 1.75 mm and H = 1.5 mm, there is a decrease in the evaporation rate by 49% from the liquid volume V = 0.5 μL (Qe = 0.93 pL s−1) to V = 10−3 μL (Qe = 0.47 pL s−1). Furthermore, the effect of the reservoir geometry is very significant, as evaporation rates in expanding reservoirs at a given liquid volume can be much greater than rates in contracting reservoirs: for the liquid volume V = 0.3 μL, the D2 = 1.75 mm, H = 1.5 mm reservoir (Qe = 0.80 pL s−1) has an evaporation rate 260% greater than the D2 = 0.75 mm reservoir of the same height (Qe = 0.30 pL s−1). The overall dependence on reservoir height is also important, as decreases in a reservoir height of H = 2.0 mm to H = 1.5 mm leads to increases in Qe of 38% and 34% for the D2 = 1.75 mm and D2 = 0.75 mm reservoir (V = 0.1 μL), respectively.

Fig. 5
(A) Effect of reservoir geometry on evaporation rate for constant contact angles of θ1 = 25° and θ2 = 50° and D1 = 1.25 mm. (B) Effect of reservoir contact angles on evaporation rates from two reservoir geometries with ...

In addition to the effects of the reservoir geometry, the surface properties of the reservoir also have an effect on the evaporation rate, albeit to a lesser degree. Fig. 5(b) displays this effect, as decreases in both θ1 and θ2, that is, increased hydrophilicity, will serve to increase the evaporation rate from both expanding and contracting reservoirs. This effect is more pronounced at higher reservoir liquid volumes associated with complete menisci, where Qe depends only on θ2, and is mainly due to the increase in the interfacial area with decreasing liquid/reservoir contact angles. Fig. 5(b) also displays the dependence of θ2 on Vc, as increases in θ2 lead to decreases in Vc for all reservoir geometries.

For the reservoirs studied here, changes in temperature, pressure, and humidity on the evaporation rate are independent of the reservoir liquid volume and will not be discussed in detail. It should be noted that to ensure accuracy regarding the results in Fig. 5, the maximum reservoir liquid volume was limited so that the average distance from the meniscus to the top of the reservoir was less than H/2, and this will be discussed below.

Prediction of meniscus evolution during evaporation

Given a specific reservoir geometry (D1, D2, H) and specific experimental conditions (T, p, RH, Vo), one only needs to know values of θ1 and θ2 to predict the evolution of an evaporating meniscus over time. Standard experimental methods cannot be used to measure the value of θ2 due to the confined geometry of the reservoir sidewalls. However, because the critical volume of the reservoir depends only on θ2, knowledge about the initial rupture diameter can be used to gain insight regarding predicted values of θ1 and θ2. Fig. 6 shows a contour plot of the predicted dependence of θ1 and θ2 on Dcap(tc) for a reservoir with D1 = 1.2 mm, D2 = 1.6 mm, and H = 1.8 mm. A specific value of θ2 determines Vc for a given reservoir geometry, which in turn establishes a relationship between θ1 and Dcap(tc). Using experimental values for Dcap(tc), such as that shown in Fig. 3(c), a functional relationship for θ1 in terms of θ2 may then be determined. Thus, one needs only experimental values for Dcap(tc) and an initial guess for θ2 to predict the parameters V, ro, zo, R1 and R2 as a function of time.

Fig. 6
Contour plot displaying the predicted initial rupture diameter Dcap(tc) of a meniscus as a function of θ1 and θ2 for a reservoir with D1 = 1.2 mm, D2 = 1.6 mm, and H = 1.8 mm.

To predict the evolution of a meniscus during evaporation, the time rate of change for the volume within a reservoir is discretized using a first order forward difference

Vn+1=VnQeΔt
(5)

where Vn represents the volume of the liquid at time t = nΔt. The simulation is initialized with a liquid volume Vo, and values of Δt are adjusted such that QeΔt < 10−4 μL. An initial guess for θ2 and the experimental value for Dcap(tc) are used to determine a matching value for θ1. These values are then used along with the experimental values of D1, D2, H, T, p, and RH as inputs to a series of CFD simulations to determine the parameters (a1, a2, a3, b1, b2, b3) associated with eqns (3) and (4). These parameters determine the relationship between Qe and V, which in turn is used in eqn (5) to calculate V as a function of time. The predicted values of Dcap as a function of time (for t > tc) can then be calculated using the geometric model (Appendix I in the ESI). Typically, several iterations involving a guess for θ2 are required for a good fit to the experimental data, where each iteration involves calculation of the parameters associated with eqns (3) and (4).

Fig. 7 compares experimental measurements to predicted values of Dcap as a function of time for several reservoir geometries and constructs (PDMS reservoir, glass or PDMS floor). The parameters associated with the results shown in Fig. 7 are displayed in Table 1. Aside from small scale differences in Dcap persisting for roughly 3 minutes after the meniscus rupture, the model does very well with respect to predicting experimental measurements. The geometric/CFD model predicts the time of meniscus rupture and overall time of the experiment (tend) within 1.5 and 0.5 minutes, respectively, for all cases (16 total experiments). The effect of reservoir geometry on the overall evaporation rate may be seen from the results in Fig. 7. For example, the difference between an expanding and contracting reservoir can clearly be seen between the case (C) and (F), where times for reservoir dry out for a 0.5 μL droplet in an expanding reservoir was 16.2 min compared with 64.7 minutes for dry out of a 0.65 μL droplet in a smaller contracting reservoir. The reservoir in case (F) had 30% more liquid, yet took nearly 300% more time to dry out under roughly the same conditions. The initial slopes of Dcap shown in Fig. 7 are a complicated relationship between the reservoir geometries, surface conditions, and external conditions. This behavior is unlike the evaporation of water from shallow wells, where time derivatives of Dcap are independent of the reservoir geometry and volume.35 The discrepancy between the two cases is most likely due to the pinning of the meniscus around the upper rim of a shallow reservoir while the lower meniscus edge recedes freely. Neither the upper or lower edge of menisci in the reservoirs in this study has been observed to be static. If either of the two contact lines of the meniscus were to be pinned in a microreservoir, the combined CFD/geometric model would likely not show the accuracy displayed in Fig. 7, since the geometric model assumes a contact line that can move freely along both the reservoir floor and sidewalls.

Fig. 7
Experimental and predicted data showing the inner ruptured meniscus (Dcap) as a function of time. The parameters for the predictions are shown in Table 1.
Table 1
Experimental and predicted parameters for the results shown in Fig. 7

A total of 16 experiments were fit to predictions such as those shown in Fig. 7 corresponding to a range of parameters from 0.77≤D1,D2 ≤ 2.2 mm, 1.0 ≤ H ≤ 3.5 mm, 9.8 ≤θ1 ≤ 30°, and 41 ≤ θ2 ≤ 58°. The predicted contact angle of water and a glass floor was θ1,glass = 14 ± 5.6°, whereas the predicted value pertaining to a PDMS floor was θ1,PDMS = 26 ± 3.4° (based on 8 experiments each). These values compare very well to experimental contact angle measurements of clean glass and PDMS oxidized with the conditions in this experiment (25 W, 900 mTorr air, 20 s).41 The predicted contact angle of water and the PDMS reservoir was θ2,PDMS = 49 ± 4.7°, where there was no significant difference between θ2,PDMS values pertaining to PDMS/PDMS and PDMS/glass reservoirs. The difference between θ1,PDMS and θ2,PDMS may be attributed to several factors, including the increase in surface roughness of the reservoir sidewalls as a result of the punching process42 and the lower efficiency of plasma penetration to the smaller diameter cavity. The small variance in the fitted parameters θ1,PDMS, θ1, glass, and θ2,PDMS with respect to the large range of experimental parameters in this study highlight the effectiveness of the combined geometric/CFD model.

Regardless of the discrepancy between θ1,PDMS and θ2,PDMS, these averaged parameters may be used to predict the evolution of a meniscus with high accuracy. Table 2 displays the errors associated with tc and tend for the evolution of a meniscus in the reservoirs shown in Fig. 7 using the averaged parameters θ1,PDMS, θ1,glass, and θ2,PDMS given above. As expected, the values of Δtc and Δtend are strongly dependent on Δθ2, the difference between the averaged value θ2,PDMS and the value θ2 shown in Table 1. Because of the small dependence of Qe on θ1 and θ2, as shown in Fig. 5(b), the error in Δtend is much smaller than Δtc, accurately predicting reservoir dry out within a minute for all cases.

Table 2
Errors of prediction for the cases shown in Fig. 7 by utilizing the averaged parameters θ1,glass = 14°, θ1,PDMS = 26°, and θ2,PDMS = 49°

Meniscus stability and model accuracy

It is hypothesized that there exists a stability criterion for a ruptured meniscus (as shown in Fig. 1) similar to that of a 2D meniscus lying on a solid edge with angle 2α. Following the example presented by Langbein,43 it is expected that for values of 2α + θ1 + θ2 < π the ruptured meniscus will be stable in the concave form shown in Fig. 1. Beyond this stability criterion, for values of 2α + θ1 + θ2 > π, the ruptured meniscus is expected to exist in either a convex form or a shape not readily described mathematically. For experimental systems (PDMS/PDMS) made hydrophilic with a lower quality plasma (200 mTorr, 15 W, 20 s), the meniscus indeed ruptures to a state not defined by the geometric model proposed here; however, the boundary between the two states of stability has not been resolved. Animations of an unstable rupture can be viewed in the ESI.

Because of the dilute nature of the fluorescent bead solution used in this study, it is doubtful that θ1 or θ2 will change as a function of reservoir liquid volume. However, because the liquid within the reservoir reaches a point of infinite concentration as t → tend, it is unlikely that solutions containing realistic solute concentrations would display constant apparent contact angles over the duration of the experiment. To model these systems, one would only need a relationship between θ1 or θ2 and the solute concentration, which in principle are obtainable through experiment. As long as solute precipitation does not interfere with the position of the mensiscus, or the stability criterion shown above is not violated, these systems could most likely be predicted with the same accuracy as the cases shown in Fig. 7.

In this study, care was taken to ensure that the average meniscus position was at a z-position lower than H/2 at the onset of the experiment. As the meniscus height approaches the top of the reservoir, or VVres, the evaporation rate calculated by eqn (2) increases in inverse proportion to the distance from the meniscus to the reservoir outlet. Under those conditions eqn (2) will over predict the evaporation rate from the meniscus, resulting in an under prediction of tc and tend. Because eqn (2) achieves significant accuracy as the liquid volume decreases, the model can be used for all systems with VVres to estimate the minimum time for reservoir dry out. To accurately describe the evolution of a meniscus in a microreservoir for conditions described above, one must (1) develop a model that accounts for the steady-state diffusion of water vapor into the space above the reservoir and (2) develop a function that gives a smooth transition between the two models as a function of the liquid volume in the reservoir.

Applications to lab-on-a-chip systems

When the situation arises for which long incubation times (> 10 min) are required for a continuous flow μFN based assay, it is important that analyte concentrations within the inlet reservoir remain relatively constant. For continuous flow assays, the volume of liquid in an inlet reservoir can be accurately determined using eqn (5) modified with a term to account for the loss of fluid via volumetric flow through a μFN,

Vn+1=Vn(Qe+Qf)Δt
(6)

where Qf is the volumetric flow rate drawn from the inlet reservoir.

Rather than utilizing complex fabrication strategies or inclusion of unnecessary experimental protocols, the techniques in this study can be used to design inlet reservoirs such that the inlet concentration is within acceptable limits over a predetermined time. The simplicity of this approach is in significant contrast with some previous technologies; for example, the assay system developed by Zimmermann et al. requires an inlet reservoir to the μFN to be cooled to 2 °C above the dew point in order for a 0.6 μL water droplet to persist for over an hour.21 A reservoir such as that shown in Fig. 7, case (F), achieves the same result without the need to control the external environment. We have also fabricated reservoirs with D1 = 1.1 mm, D2 = 1.6 mm, and H = 3.5 mm that allow a 0.6 μL droplet to persist for over 90 minutes. For that case, it is predicted that (provided Vres/Qf [dbl greater-than sign]tend) at t = 20 min, V/Vo = 0.9, which is sufficient for most analytical measurements. This same approach may also be used to enhance evaporation in systems that require analyte pre-concentration for detection. An inlet reservoir designed with dimensions that produce higher evaporation rates used in microfluidic based assays would lead to systems with lower detection limits, as the reservoir solution could be concentrated by an order of magnitude within several minutes.

The predictive techniques used in this study can also be used for passive pumping applications based on capillary forces. The capillary pressure difference between a static air/liquid interface is described by the Young–Laplace equation

Δp=pp1=γlv(1R1+1R2)
(7)

where pl represents the pressure in the liquid phase at the interface. Two reservoirs connected by a common μFN may be used to drive fluid between one another, provided that the two principal radii of curvature of liquid in the two reservoirs are different. This system would provide several advantages over the type of passive pump developed by Walker and Beebe,25 namely (1) the air/liquid interface is located entirely within the reservoir, reducing the chance of unwanted contact with foreign materials and (2) the values of R1 for reservoirs with low liquid volumes can be much smaller than the characteristic size of the reservoir, allowing for very high differential pressures. From the geometric model proposed here R1 and R2 can be calculated for a given reservoir liquid volume, thus the volumetric flow rate between two reservoirs can be accurately determined using eqn (7), while the volume of liquid in each reservoir can be determined using eqn (6). Initial studies show fluid velocities in a μFN between two reservoirs greater than 1 mm s−1, corresponding to pressure drops Δp > 8600 Pa. The ease of fabrication and straightforward prediction of fluidic response may prove to be extremely valuable for lab-on-a-chip applications, and is the focus of future work.

Conclusion

There are many lab-on-a-chip designs that utilize a micro-reservoir either for direct chemical analysis or as a method for fluid introduction into a microfluidic network. Owing to both the small sample size and the likelihood that the liquid volume is open to the external environment, evaporation can quickly lead to detrimental results within the device. Knowledge of the rates of liquid evaporation within a microreservoir, and how that rate depends on the local conditions, is essential when considering the performance reliability of such a device.

To gain insight into the rate of evaporation in a given system, a simple geometric model is proposed that accounts for the shape of an air/liquid meniscus present in a cylindrical microreservoir with variable cross section. This model is utilized as a template for studying the rates of evaporation from the liquid volume via computational fluid dynamics. The combined geometric/CFD approach has been shown to quantitatively predict the evolution of a meniscus over time in experimental systems, and with high accuracy. This method can be used with any liquids as long as the geometry of the reservoir and solution properties obey the relation 2α + θ1 + θ2 < π. For liquid present in a microreservoir, simple methods to enhance or moderate evaporation have been discussed. Evaporation rates can be increased in microreservoirs by reducing the height of reservoir or increasing the overall diameter of the reservoir such that D2 > D1. Similarly, evaporation is readily decreased by increasing the height of the reservoir or decreasing the overall diameter such that D1 > D2.

The techniques used in this study may be used to optimize a variety of lab-on-a-chip applications that involve a holding reservoir open to the environment. Furthermore, the geometric model may be used to design passive pumping between multiple reservoirs, whose characteristics can be predicted over time as the meniscus in each reservoir evolves. The passive pumping mechanisms related to this study are the focus of future work.

Supplementary Material

Appendices

Acknowledgments

This work was supported by NIH grant number EB00726.

Footnotes

Electronic supplementary information (ESI) available: Two videos showing stable and unstable rupture. Appendices I and II. See DOI: 10.1039/b900556k

References

1. Haeberle S, Zengerle R. Lab on a Chip. 2007;7:1094–1110. [PubMed]
2. Yi CQ, Li CW, Ji SL, Yang MS. Analytica Chimica Acta. 2006;560:1–23.
3. Toner M, Irimia D. Annual Review of Biomedical Engineering. 2005;7:77–103. [PMC free article] [PubMed]
4. Zhang CS, Xu JL, Ma WL, Zheng WL. Biotechnology Advances. 2006;24:243–284. [PubMed]
5. Situma C, Hashimoto M, Soper SA. Biomolecular Engineering. 2006;23:213–231. [PubMed]
6. Aborn JH, El-Difrawy SA, Novotny M, Gismondi EA, Lam R, Matsudaira P, Mckenna BK, O’Neil T, Streechon P, Ehrlich DJ. Lab on a Chip. 2005;5:669–674. [PubMed]
7. Liu Y, Vickers JA, Henry CS. Analytical Chemistry. 2004;76:1513–1517. [PubMed]
8. Randall GC, Doyle PS. Proceedings of the National Academy of Sciences of the United States of America. 2005;102:10813–10818. [PubMed]
9. Berthier E, Warrick J, Yu H, Beebe DJ. Lab on a Chip. 2008;8:852–859. [PMC free article] [PubMed]
10. Berthier E, Warrick J, Yu H, Beebe DJ. Lab on a Chip. 2008;8:860–864. [PMC free article] [PubMed]
11. Eijkel JCT, van den Berg A. Lab on a Chip. 2005;5:1202–1209. [PubMed]
12. Litborn E, Emmer A, Roeraade J. Electrophoresis. 2000;21:91–99. [PubMed]
13. Neugebauer S, Evans SR, Aguilar ZP, Mosbach M, Fritsch I, Schuhmann W. Analytical Chemistry. 2004;76:458–463. [PubMed]
14. Moerman R, Knoll J, Apetrei C, van den Doel LR, van Dedem GWK. Analytical Chemistry. 2005;77:225–231. [PubMed]
15. Bratten CDT, Cobbold PH, Cooper JM. Analytical Chemistry. 1997;69:253–258.
16. Litborn E, Roeraade J. Journal of Chromatography B. 2000;745:137–147. [PubMed]
17. Gratzl M, Yi C. Analytical Chemistry. 1993;65:2085–2088.
18. Litborn E, Emmer A, Roeraade J. Analytica Chimica Acta. 1999;401:11–19.
19. Minarik M, Kleparnik K, Gilar M, Foret F, Miller AW, Sosic Z, Karger BL. Electrophoresis. 2002;23:35–42. [PubMed]
20. Goedecke N, Eijkel J, Manz A. Lab on a Chip. 2002;2:219–223. [PubMed]
21. Zimmermann M, Bentley S, Schmid H, Hunziker P, Delamarche E. Lab on a Chip. 2005;5:1355–1359. [PubMed]
22. Effenhauser CS, Harttig H, Kramer P. Biomedical Microdevices. 2002;4:27–32.
23. Xu Z, Zhong C, Guan Y, Chen X, Wang J, Fang Z. Lab on a Chip. 2008;8:1658–1663. [PubMed]
24. Berthier E, Beebe DJ. Lab on a Chip. 2007;7:1475–1478. [PubMed]
25. Walker GM, Beebe DJ. Lab on a Chip. 2002;2:131–134. [PubMed]
26. Chen YF, Tseng FG, ChangChien SY, Chen MH, Yu RJ, Chieng CC. Microfluidics and Nanofluidics. 2008;5:193–203.
27. Delamarche E, Bernard A, Schmid H, Michel B, Biebuyck H. Science. 1997;276:779–781. [PubMed]
28. Bernard A, Michel B, Delamarche E. Analytical Chemistry. 2001;73:8–12. [PubMed]
29. Murphy BM, He XY, Dandy D, Henry CS. Analytical Chemistry. 2008;80:444–450. [PMC free article] [PubMed]
30. Deegan RD, Bakajin O, Dupont TF, Huber G, Nagel SR, Witten TA. Physical Review E. 2000;62:756–765. [PubMed]
31. Deegan RD, Bakajin O, Dupont TF, Huber G, Nagel SR, Witten TA. Nature. 1997;389:827–829.
32. Hu H, Larson RG. Journal of Physical Chemistry B. 2002;106:1334–1344.
33. Deng Y, Zhu XY, Kienlen T, Guo A. Journal of the American Chemical Society. 2006;128:2768–2769. [PubMed]
34. Chen CT, Chieng CC, Tseng FG. Journal of Microelectromechanical Systems. 2007;16:1209–1218.
35. Chen CT, Tseng FG, Chieng CC. Sensors and Actuators a-Physical. 2006;130:12–19.
36. van den Doel LR, van Vliet LJ. Applied Optics. 2001;40:4487–4500. [PubMed]
37. Hjelt KT, van den Doel R, Lubking W, Vellekoop MJ. Sensors and Actuators a-Physical. 2000;85:384–389.
38. Rieger B, van den Doel LR, van Vliet LJ. Physical Review E. 2003;68 [PubMed]
39. Ryskin G, Leal LG. Journal of Fluid Mechanics. 1984;148:1–17.
40. Vickers JA, Caulum MM, Henry CS. Analytical Chemistry. 2006;78:7446–7452. [PubMed]
41. Bhattacharya S, Datta A, Berg JM, Gangopadhyay S. Journal of Microelectromechanical Systems. 2005;14:590–597.
42. Shastry A, Case MJ, Bohringer KF. Langmuir. 2006;22:6161–6167. [PubMed]
43. Langbein D. Journal of Fluid Mechanics. 1990;213:251–265.