Effects of viscoelasticity on droplet dynamics and break-up in microfluidic T-Junctions: a lattice Boltzmann study

Anupam Gupta; Mauro Sbragaglia

doi:10.1140/epje/i2016-16006-9

2022 Impact factor 1.8

Soft Matter and Biological Physics

Multi-scale phenomena in complex flows and flowing matter

Eur. Phys. J. E (2016) 39: 6
https://doi.org/10.1140/epje/i2016-16006-9

Regular Article

Effects of viscoelasticity on droplet dynamics and break-up in microfluidic T-Junctions: a lattice Boltzmann study

Anupam Gupta^* and Mauro Sbragaglia

Department of Physics and INFN, University of “Tor Vergata”, Via della Ricerca Scientifica 1, 00133, Rome, Italy

^* e-mail: anupam1509@gmail.com

Received: 31 July 2015
Accepted: 29 December 2015
Published online: 27 January 2016

Abstract

The effects of viscoelasticity on the dynamics and break-up of fluid threads in microfluidic T-junctions are investigated using numerical simulations of dilute polymer solutions at changing the Capillary number (Ca), i.e. at changing the balance between the viscous forces and the surface tension at the interface, up to Ca 3×10^-2. A Navier-Stokes (NS) description of the solvent based on the lattice Boltzmann models (LBM) is here coupled to constitutive equations for finite extensible non-linear elastic dumbbells with the closure proposed by Peterlin (FENE-P model). We present the results of three-dimensional simulations in a range of Ca which is broad enough to characterize all the three characteristic mechanisms of break-up in the confined T-junction, i.e. squeezing, dripping and jetting regimes. The various model parameters of the FENE-P constitutive equations, including the polymer relaxation time and the finite extensibility parameter L², are changed to provide quantitative details on how the dynamics and break-up properties are affected by viscoelasticity. We will analyze cases with Droplet Viscoelasticity (DV), where viscoelastic properties are confined in the dispersed (d) phase, as well as cases with Matrix Viscoelasticity (MV), where viscoelastic properties are confined in the continuous (c) phase. Moderate flow-rate ratios Q O(1) of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, as the flow driving the break-up process upstream of the emerging thread can be sensibly perturbed by the polymer stresses.