Droplets in Microchannels: Dynamical Properties of the Lubrication Film

International audienceWe study the motion of droplets in a confined, micrometric geometry, by focusing on the lubrication film between droplet and wall. When capillary forces dominate, the lubrication film thickness evolves non linearly with the capillary number due to viscous dissipation between meniscus and wall. However, this film may become thin enough (tens of nanometers) that intermolecular forces come into play and affect classical scalings. Our experiments yield highly resolved topographies of the shape of the interface and allow us to bring new insights into droplet dynamics in microfluidics. We report the novel characterization of two dynamical regimes as the capillary number increases: (i) at low capillary numbers, the film thickness is constant and set by the disjoinging pressure, while (ii) above a critical capillary number, the interface behavior is well described by a viscous scenario. At a high surfactant concentration, structural effects lead to the formation of patterns on the interface , which can be used to trace the interface velocity that yield direct confirmation of boundary condition in viscous regime. The dynamics of a droplet confined between solid walls and pushed by a surrounding liquid is an old problem, however recent theories are still being developed to describe unexplored regimes and experimental characterizations are still lacking to shed light on these novel developments. A complete understanding of the droplet velocity calls for accurate knowledge of the dissipation mechanisms involved, particularly in the lubrication film. Our understanding of the lubrication properties of menisci travelling in confined geometries has been steadily refined since the pioneering work of Taylor & Saffman [1]. Notably, the influence of the lubrication film left along the wall by the moving meniscus was first taken into account by Bretherton, who investigated the motion of an inviscid bubble in a cylindrical tube [2]. Far from the meniscus, this dynamical film reaches a uniform thickness h ∞ , related to the bubble velocity through the capillary number Ca = µ f U d /γ, where U d is the bubble velocity , µ f the viscosity of the continuous phase, and γ the surface tension. When the capillary pressure dominates over the viscous stress, i.e. in the regime where Ca 1, the thickness of the film follows h Breth = 1.34 r Ca 2/3 , where r is the radius of the capillary tube. Besides bubbles , the case of viscous droplets remains however largely unexplored. A recent theoretical advance in the field by Hodges et al. [3] shows by numerical calculations of the whole flow pattern that significant corrections in the thickness of lubrication films can arise at very low Ca. Furthermore, the regime of the Bretherton theory is only valid for lubrication films thicker than the molecular sizes or than the range of interfacial interactions. The typical velocities and lengthscales involved in droplet-based microfluidics would lead to lubrication films h ∞ ∼

Simulations of Solid-on-Solid Models of Spreading of Viscous Droplets

We have studied the dynamics of spreading of viscous non-volatile fluids on surfaces by MC simulations of SOS models. We have concentrated on the complete wetting regime, with surface diffusion barriers neglected for simplicity. First, we have performed simulations for the standard SOS model. Formation of a single precursor layer, and a density profile with a spherical cap shaped center surrounded by Gaussian tails can be reproduced with this model. Dynamical layering (DL), however, only occurs with a very strongly attractive van der Waals type of substrate potential. To more realistically describe the spreading of viscous liquid droplets, we introduce a modified SOS model. In the new model, tendency for DL and the effect of the surface potential are in part embedded into the dynamics of the model. This allows a relatively simple description of the spreading under different conditions, with a temperature like parameter which strongly influences the droplet morphologies. Both rounded droplet shapes and DL can easily be reproduced with the model. Furthermore, the precursor width increases proportional to the square root of time, in accordance with experimental observations. PACS: 68.10.Gw, 05.70.Ln, 61.20.Ja.Comment: to appear in Physica A (1994), standard LaTex, 20 page

Molecular ordering of precursor films during spreading of tiny liquid droplets

In this work we address a novel feature of spreading dynamics of tiny liquid droplets on solid surfaces, namely the case where the ends of the molecules feel different interactions to the surface. We consider a simple model of dimers and short chain--like molecules which cannot form chemical bonds with the surface. We study the spreading dynamics by Molecular Dynamics techniques. In particular, we examine the microscopic structure of the time--dependent precursor film and find that in some cases it can exhibit a high degree of local order. This order persists even for flexible chains. Our results suggest the possibility of extracting information about molecular interactions from the structure of the precursor film.Comment: 4 pages, revtex, no figures, complete file available from ftp://rock.helsinki.fi/pub/preprints/tft/ or at http://www.physics.helsinki.fi/tft/tft_preprints.html (to appear in Phys. Rev. E Rapid Comm.

Kinetics of Anchoring of Polymer Chains on Substrates with Chemically Active Sites

We consider dynamics of an isolated polymer chain with a chemically active end-bead on a 2D solid substrate containing immobile, randomly placed chemically active sites (traps). For a particular situation when the end-bead can be irreversibly trapped by any of these sites, which results in a complete anchoring of the whole chain, we calculate the time evolution of the probability $P_{ch}(t)$ that the initially non-anchored chain remains mobile until time $t$. We find that for relatively short chains $P_{ch}(t)$ follows at intermediate times a standard-form 2D Smoluchowski-type decay law $ln P_{ch}(t) \sim - t/ln(t)$, which crosses over at very large times to the fluctuation-induced dependence $ln P_{ch}(t) \sim - t^{1/2}$, associated with fluctuations in the spatial distribution of traps. We show next that for long chains the kinetic behavior is quite different; here the intermediate-time decay is of the form $ln P_{ch}(t) \sim - t^{1/2}$, which is the Smoluchowski-type law associated with subdiffusive motion of the end-bead, while the long-time fluctuation-induced decay is described by the dependence $ln P_{ch}(t) \sim - t^{1/4}$, stemming out of the interplay between fluctuations in traps distribution and internal relaxations of the chain.Comment: Latex file, 19 pages, one ps figure, to appear in PR

Binary separation in very thin nematic films: thickness and phase coexistence

The behavior as a function of temperature of very thin films (10 to 200 nm) of pentylcyanobiphenyl (5CB) on silicon substrates is reported. In the vicinity of the nematic/isotropic transition we observe a coexistence of two regions of different thicknesses: thick regions are in the nematic state while thin ones are in the isotropic state. Moreover, the transition temperature is shifted downward following a 1/h^2 law (h is the film thickness). Microscope observations and small angle X-ray scattering allowed us to draw a phase diagram which is explained in terms of a binary first order phase transition where thickness plays the role of an order parameter.Comment: 5 pages, 3 figures, submitted to PRL on the 26th of Apri

Spreading of a Macroscopic Lattice Gas

We present a simple mechanical model for dynamic wetting phenomena. Metallic balls spread along a periodically corrugated surface simulating molecules of liquid advancing along a solid substrate. A vertical stack of balls mimics a liquid droplet. Stochastic motion of the balls, driven by mechanical vibration of the corrugated surface, induces diffusional motion. Simple theoretical estimates are introduced and agree with the results of the analog experiments, with numerical simulation, and with experimental data for microscopic spreading dynamics.Comment: 19 pages, LaTeX, 9 Postscript figures, to be published in Phy. Rev. E (September,1966

Dynamics of Spreading of Chainlike Molecules with Asymmetric Surface Interactions

In this work we study the spreading dynamics of tiny liquid droplets on solid surfaces in the case where the ends of the molecules feel different interactions with respect to the surface. We consider a simple model of dimers and short chainlike molecules that cannot form chemical bonds with the surface. We use constant temperature Molecular Dynamics techniques to examine in detail the microscopic structure of the time dependent precursor film. We find that in some cases it can exhibit a high degree of local order that can persist even for flexible chains. Our model also reproduces the experimentally observed early and late-time spreading regimes where the radius of the film grows proportional to the square root of time. The ratios of the associated transport coefficients are in good overall agreement with experiments. Our density profiles are also in good agreement with measurements on the spreading of molecules on hydrophobic surfaces.Comment: 12 pages, LaTeX with APS macros, 21 figures available by contacting [email protected], to appear in Phys. Rev.

Molecular Weight Dependence of Spreading Rates of Ultrathin Polymeric Films

We study experimentally the molecular weight $M$ dependence of spreading rates of molecularly thin precursor films, growing at the bottom of droplets of polymer liquids. In accord with previous observations, we find that the radial extension R(t) of the film grows with time as R(t) = (D_{exp} t)^{1/2}. Our data substantiate the M-dependence of D_{exp}; we show that it follows D_{exp} \sim M^{-\gamma}, where the exponent \gamma is dependent on the chemical composition of the solid surface, determining its frictional properties with respect to the molecular transport. In the specific case of hydrophilic substrates, the frictional properties can be modified by the change of the relative humidity (RH). We find that \gamma \approx 1 at low RH and tends to zero when RH gets progressively increased. We propose simple theoretical arguments which explain the observed behavior in the limits of low and high RH.Comment: 4 pages, 2 figures, to appear in PR

Controlling the temperature sensitivity of DNA-mediated colloidal interactions through competing linkages

We propose a new strategy to improve the self-assembly properties of DNA-functionalised colloids. The problem that we address is that DNA-functionalised colloids typically crystallize in a narrow temperature window, if at all. The underlying reason is the extreme sensitivity of DNA-mediated interactions to temperature or other physical control parameters. We propose to widen the window for colloidal crystallization by exploiting the competition between DNA linkages with different nucleotide sequences, which results in a temperature-dependent switching of the dominant bond type. Following such a strategy, we can decrease the temperature dependence of DNA-mediated self assembly to make systems that can crystallize in a wider temperature window than is possible with existing systems of DNA functionalised colloids. We report Monte Carlo simulations that show that the proposed strategy can indeed work in practice for real systems and specific, designable DNA sequences. Depending on the length ratio of the different DNA constructs, we find that the bond switching is either energetically driven (equal length or `symmetric' DNA) or controlled by a combinatorial entropy gain (`asymmetric' DNA), which results from the large number of possible binding partners for each DNA strand. We provide specific suggestions for the DNA sequences with which these effects can be achieved experimentally

DNA-templated assembly of droplet-derived PEG microtissues

Patterning multiple cell types is a critical step for engineering functional tissues, but few methods provide three-dimensional positioning at the cellular length scale. Here, we present a “bottom-up” approach for fabricating multicellular tissue constructs that utilizes DNA-templated assembly of 3D cell-laden hydrogel microtissues. A flow focusing-generated emulsion of photopolymerizable prepolymer is used to produce 100 μm monodisperse microtissues at a rate of 100 Hz (10[superscript 5] h[superscript −1]). Multiple cell types, including suspension and adherently cultured cells, can be encapsulated into the microtissues with high viability ([similar]97%). We then use a DNA coding scheme to self-assemble microtissues “bottom-up” from a template that is defined using “top-down” techniques. The microtissues are derivatized with single-stranded DNA using a biotin–streptavidin linkage to the polymer network, and are assembled by sequence-specific hybridization onto spotted DNA microarrays. Using orthogonal DNA codes, we achieve multiplexed patterning of multiple microtissue types with high binding efficiency and >90% patterning specificity. Finally, we demonstrate the ability to organize multicomponent constructs composed of epithelial and mesenchymal microtissues while preserving each cell type in a 3D microenvironment. The combination of high throughput microtissue generation with scalable surface-templated assembly offers the potential to dissect mechanisms of cell–cell interaction in three dimensions in healthy and diseased states, as well as provides a framework for templated assembly of larger structures for implantation