A numerical study of droplet splitting in branched T-shaped microchannel using the two-phase level-set method

2021 ◽  
Vol 13 (11) ◽  
pp. 168781402110454
Author(s):  
Mohammad Raad ◽  
Sajad Rezazadeh ◽  
Habib Jalili ◽  
Davod Abbasinezhad Fallah

Droplet splitting as a significant feature of droplet-based microfluidic systems has been widely employed in biotechnology, biomedical engineering, tissue engineering, and it has been preferred over continuous flow systems. In the present paper, two-dimensional numerical simulations have been done to examine the asymmetrical droplet splitting process. The two-phase level set method (LSM) has been predicted to analyze the mechanism of droplet formation and droplet splitting in immiscible liquid/liquid two-phase flow in the branched T-junction microchannel. Governing equations on flow field have been discretized and solved using finite element-based COMSOL Multiphysics software (version 5.3a). Obtained numerical results were validated by experimental data reported in the literature which show acceptable agreement. The model was developed to simulate the mechanism of droplet splitting at the branched T-junction microchannel. This study provides a passive technique to asymmetrically split up microdroplets at the downstream T-junctions. The results show that outlet branches’ pressure gradient affects the droplet splitting. Specifically, it has been shown that the splitting ratio increases by increasing the length ratio, and equal droplet splitting can be achieved where the ratio is LL/ Lu = 1. We have used two outlet branches having the same width but different lengths to create the required pressure gradient. As the length ratio of the outlet branches increases, the diameter ratio increases as well.

2008 ◽  
Vol 11 (4-6) ◽  
pp. 221-235 ◽  
Author(s):  
S. P. van der Pijl ◽  
A. Segal ◽  
C. Vuik ◽  
P. Wesseling

2014 ◽  
Vol 100 ◽  
pp. 138-154 ◽  
Author(s):  
Lanhao Zhao ◽  
Jia Mao ◽  
Xin Bai ◽  
Xiaoqing Liu ◽  
Tongchun Li ◽  
...  

2016 ◽  
Vol 18 (1) ◽  
pp. 31-52 ◽  
Author(s):  
Peter Frolkovič ◽  
Dmitriy Logashenko ◽  
Christian Wehner

Geofluids ◽  
2017 ◽  
Vol 2017 ◽  
pp. 1-11 ◽  
Author(s):  
Yunfeng Dai ◽  
Zhifang Zhou ◽  
Jin Lin ◽  
Jiangbo Han

To describe accurately the flow characteristic of fracture scale displacements of immiscible fluids, an incompressible two-phase (crude oil and water) flow model incorporating interfacial forces and nonzero contact angles is developed. The roughness of the two-dimensional synthetic rough-walled fractures is controlled with different fractal dimension parameters. Described by the Navier–Stokes equations, the moving interface between crude oil and water is tracked using level set method. The method accounts for differences in densities and viscosities of crude oil and water and includes the effect of interfacial force. The wettability of the rough fracture wall is taken into account by defining the contact angle and slip length. The curve of the invasion pressure-water volume fraction is generated by modeling two-phase flow during a sudden drainage. The volume fraction of water restricted in the rough-walled fracture is calculated by integrating the water volume and dividing by the total cavity volume of the fracture while the two-phase flow is quasistatic. The effect of invasion pressure of crude oil, roughness of fracture wall, and wettability of the wall on two-phase flow in rough-walled fracture is evaluated.


2016 ◽  
Vol 19 (2) ◽  
pp. 329-353 ◽  
Author(s):  
Florian Kummer ◽  
Tim Warburton

AbstractIn two-phase flow simulations, a difficult issue is usually the treatment of surface tension effects. These cause a pressure jump that is proportional to the curvature of the interface separating the two fluids. Since the evaluation of the curvature incorporates second derivatives, it is prone to numerical instabilities. Within this work, the interface is described by a level-set method based on a discontinuous Galerkin discretization. In order to stabilize the evaluation of the curvature, a patch-recovery operation is employed. There are numerous ways in which this filtering operation can be applied in the whole process of curvature computation. Therefore, an extensive numerical study is performed to identify optimal settings for the patch-recovery operations with respect to computational cost and accuracy.


Sign in / Sign up

Export Citation Format

Share Document