Thermocapillary motion of a Newtonian drop in a dilute viscoelastic fluid

• Published in: Journal of non-Newtonian fluid mechanics 2019-08, Vol.270, p.8-22
• Main Authors: Capobianchi, Paolo, Pinho, Fernando T., Lappa, Marcello, Oliveira, Mónica S.N.
• Format: Article
• Language: English
• Subjects: Deborah number; Deformation; Droplet dynamics; Droplets; Formability; Gravitational effects; Interfacial stresses; LS-VOF; Mathematical models; Newtonian fluids; Non Newtonian fluids; Parameters; Stagnation point; Temperature distribution; Temperature gradients; Thermocapillary migration; Viscoelastic effects; Viscoelastic fluids; Viscoelasticity; Viscosity; Viscosity ratio
• ISSN: 0377-0257; 1873-2631
• DOI: 10.1016/j.jnnfm.2019.06.006

•The thermocapillary motion of a Newtonian drop in a viscoelastic fluid is investigated numerically.•Viscoelastic stresses concentrate at the rear stagnation point of the drop.•The drop is found to deform into a prolate ellipsoid. For sufficiently large Deborah numbers, a pointed end is also observed.•The drop migration speed is found to be affected by the presence of viscoelastic stresses. In this work we investigate the role played by viscoelasticity on the thermocapillary motion of a deformable Newtonian droplet embedded in an immiscible, otherwise quiescent non-Newtonian fluid. We consider a regime in which inertia and convective transport of energy are both negligible (represented by the limit condition of vanishingly small Reynolds and Marangoni numbers) and free from gravitational effects. A constant temperature gradient is maintained by keeping two opposite sides of the computational domain at different temperatures. Consequently the droplet experiences a motion driven by the mismatch of interfacial stresses induced by the non-uniform temperature distribution on its boundary. The departures from the Newtonian behaviour are quantified via the “thermal” Deborah number, DeT and are accounted for by adopting either the Oldroyd-B model, for relatively small DeT, or the FENE-CR constitutive law for a larger range of DeT. In addition, the effects of model parameters, such as the concentration parameter c=1−β (where β is the viscoelastic viscosity ratio), or the extensibility parameter, L2, have been studied numerically using a hybrid volume of fluid-level set method. The numerical results show that the steady-state droplet velocity behaves as a monotonically decreasing function of DeT, whilst its shape deforms prolately. For increasing values of DeT, the viscoelastic stresses show the tendency to be concentrated near the rear stagnation point, contributing to an increase in its local interface curvature.