No Arabic abstract
We present an analytical study, validated by numerical simulations, of electroosmotic flow in a Hele-Shaw cell with non-uniform surface charge patterning. Applying the lubrication approximation and assuming thin electric double layer, we obtain a pair of uncoupled Poisson equations which relate the pressure and the stream function, respectively, to gradients in the zeta potential distribution parallel and perpendicular to the applied electric field. We solve the governing equations for the fundamental case of a disk with uniform zeta potential and show that the flow-field in the outer region takes the form of a pure dipole. We illustrate the ability to generate complex flow-fields around smooth convex regions by superposition of such disks with uniform zeta potential and a uniform pressure driven flow. This method may be useful for future on-chip devices, allowing flow control without the need for mechanical components.
We adopt a boundary integral method to study the dynamics of a translating droplet confined in a Hele-Shaw cell in the Stokes regime. The droplet is driven by the motion of the ambient fluid with the same viscosity. We characterize the three-dimensional (3D) nature of the droplet interface and of the flow field. The interface develops an arc-shaped ridge near the rear-half rim with a protrusion in the rear and a laterally symmetric pair of higher peaks; this pair of protrusions has been identified by recent experiments (Huerre et al., Phys. Rev. Lett., vol. 115 (6), 2015, 064501) and predicted asymptotically (Burgess and Foster, Phys. Fluids A, vol. 2 (7), 1990, pp. 1105-1117). The mean film thickness is well predicted by the extended Bretherton model (Klaseboer et al., Phys. Fluids, vol. 26 (3), 2014, 032107) with fitting parameters. The flow in the streamwise wall-normal middle plane is featured with recirculating zones, which are partitioned by stagnation points closely resembling those of a two-dimensional droplet in a channel. Recirculation is absent in the wall-parallel, unconfined planes, in sharp contrast to the interior flow inside a moving droplet in free space. The preferred orientation of the recirculation results from the anisotropic confinement of the Hele-Shaw cell. On these planes, we identify a dipolar disturbance flow field induced by the travelling droplet and its $1/r^2$ spatial decay is confirmed numerically. We pinpoint counter-rotating streamwise vortex structures near the lateral interface of the droplet, further highlighting the complex 3D flow pattern.
We study the spreading and leveling of a gravity current in a Hele-Shaw cell with flow-wise width variations as an analog for flow {in fractures and horizontally heterogeneous aquifers}. Using phase-plane analysis, we obtain second-kind self-similar solutions to describe the evolution of the gravity currents shape during both the spreading (pre-closure) and leveling (post-closure) regimes. The self-similar theory is compared to numerical simulations of the partial differential equation governing the evolution of the currents shape (under the lubrication approximation) and to table-top experiments. Specifically, simulations of the governing partial differential equation from lubrication theory allow us to compute a pre-factor, which is textit{a priori} arbitrary in the second-kind self-similar transformation, by estimating the time required for the current to enter the self-similar regime. With this pre-factor calculated, we show that theory, simulations and experiments agree well near the propagating front. In the leveling regime, the currents memory resets, and another self-similar behavior emerges after an adjustment time, which we estimate from simulations. Once again, with the pre-factor calculated, both simulations and experiments are shown to obey the predicted self-similar scalings. For both the pre- and post-closure regimes, we provide detailed asymptotic (analytical) characterization of the universal current profiles that arise as self-similarity of the second kind.
We study microfluidic self digitization in Hele-Shaw cells using pancake droplets anchored to surface tension traps. We show that above a critical flow rate, large anchored droplets break up to form two daughter droplets, one of which remains in the anchor. Below the critical flow velocity for breakup the shape of the anchored drop is given by an elastica equation that depends on the capillary number of the outer fluid. As the velocity crosses the critical value, the equation stops admitting a solution that satisfies the boundary conditions; the drop breaks up in spite of the neck still having finite width. A similar breaking event also takes place between the holes of an array of anchors, which we use to produce a 2D array of stationary drops in situ.
Micro and nanodroplets have many important applications such as in drug delivery, liquid-liquid extraction, nanomaterial synthesis and cosmetics. A commonly used method to generate a large number of micro or nanodroplets in one simple step is solvent exchange (also called nanoprecipitation), in which a good solvent of the droplet phase is displaced by a poor one, generating an oversaturation pulse that leads to droplet nucleation. Despite its crucial importance, the droplet growth resulting from the oversaturation pulse in this ternary system is still poorly understood. We experimentally and theoretically study this growth in Hele-Shaw like channels by measuring the total volume of the oil droplets that nucleates out of it. In order to prevent the oversaturated oil from exiting the channel, we decorated some of the channels with a porous region in the middle. Solvent exchange is performed with various solution compositions, flow rates and channel geometries, and the measured droplets volume is found to increase with the Peclet number $Pe$ with an approximate effective power law $Vpropto Pe^{0.50}$. A theoretical model is developed to account for this finding. With this model we can indeed explain the $Vpropto Pe^{1/2}$ scaling, including the prefactor, which can collapse all data of the porous channels onto one universal curve, irrespective of channel geometry and composition of the mixtures. Our work provides a macroscopic approach to this bottom-up method of droplet generation and may guide further studies on oversaturation and nucleation in ternary systems.
The flow in a Hele-Shaw cell with a time-increasing gap poses a unique shrinking interface problem. When the upper plate of the cell is lifted perpendicularly at a prescribed speed, the exterior less viscous fluid penetrates the interior more viscous fluid, which generates complex, time-dependent interfacial patterns through the Saffman-Taylor instability. The pattern formation process sensitively depends on the lifting speed and is still not fully understood. For some lifting speeds, such as linear or exponential speed, the instability is transient and the interface eventually shrinks as a circle. However, linear stability analysis suggests there exist shape invariant shrinking patterns if the gap $b(t)$ is increased more rapidly: $displaystyle b(t)=left(1-frac{7}{2}tau mathcal{C} tright)^{-{2}/{7}}$, where $tau$ is the surface tension and $mathcal{C}$ is a function of the interface perturbation mode $k$. Here, we use a spectrally accurate boundary integral method together with an efficient time adaptive rescaling scheme, which for the first time makes it possible to explore the nonlinear limiting dynamical behavior of a vanishing interface. When the gap is increased at a constant rate, our numerical results quantitatively agree with experimental observations (Nase et al., Phys. Fluids, vol. 23, 2011, pp. 123101). When we use the shape invariant gap $b(t)$, our nonlinear results reveal the existence of $k$-fold dominant, one-dimensional, web-like networks, where the fractal dimension is reduced to almost one at late times. We conclude by constructing a morphology diagram for pattern selection that relates the dominant mode $k$ of the vanishing interface and the control parameter $mathcal{C}$.