Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFilm thickness distribution in gravity-driven pancake-shaped droplets rising in a Hele-Shaw cell
273
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study here experimentally, numerically and using a lubrication approach; the shape, velocity and lubrication film thickness distribution of a droplet rising in a vertical Hele-Shaw cell. The droplet is surrounded by a stationary immiscible fluid and moves purely due to buoyancy. A low density difference between the two mediums helps to operate in a regime with capillary number $Ca$ lying between $0.03-0.35$, where $Ca=mu_o U_d /gamma$ is built with the surrounding oil viscosity $mu_o$, the droplet velocity $U_d$ and surface tension $gamma$. The experimental data shows that in this regime the droplet velocity is not influenced by the thickness of the thin lubricating film and the dynamic meniscus. For iso-viscous cases, experimental and three-dimensional numerical results of the film thickness distribution agree well with each other. The mean film thickness is well captured by the Aussillous & Quere (2000) model with fitting parameters. The droplet also exhibits the catamaran shape that has been identified experimentally for a pressure-driven counterpart (Huerre $textit{et al}$. 2015). This pattern has been rationalized using a two-dimensional lubrication equation. In particular, we show that this peculiar film thickness distribution is intrinsically related to the anisotropy of the fluxes induced by the droplets motion.
rate research
Read More
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 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.
Droplet migration in a Hele--Shaw cell is a fundamental multiphase flow problem which is crucial for many microfluidics applications. We focus on the regime at low capillary number and three-dimensional direct numerical simulations are performed to investigate the problem. In order to reduce the computational cost, an adaptive mesh is employed and high mesh resolution is only used near the interface. Parametric studies are performed on the droplet horizontal radius and the capillary number. For droplets with an horizontal radius larger than half the channel height the droplet overfills the channel and exhibits a pancake shape. A lubrication film is formed between the droplet and the wall and particular attention is paid to the effect of the lubrication film on the droplet velocity. The computed velocity of the pancake droplet is shown to be lower than the average inflow velocity, which is in agreement with experimental measurements. The numerical results show that both the strong shear induced by the lubrication film and the three-dimensional flow structure contribute to the low mobility of the droplet. In this low-migration-velocity scenario the interfacial flow in the droplet reference frame moves toward the rear on the top and reverses direction moving to the front from the two side edges. The velocity of the pancake droplet and the thickness of the lubrication film are observed to decrease with capillary number. The droplet velocity and its dependence on capillary number cannot be captured by the classic Hele--Shaw equations, since the depth-averaged approximation neglects the effect of the lubrication film.
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}$.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
