Dynamics of particle sedimentation in viscoelastic fluids: A numerical study on particle chain in two-dimensional narrow channel

TL;DR

This study introduces a numerical method to simulate particle sedimentation in viscoelastic fluids, revealing how polymer extension limits influence chain formation and particle interactions in narrow channels.

Contribution

A novel numerical approach for simulating particle sedimentation in FENE-CR viscoelastic fluids, analyzing the effects of polymer extension limits on chain formation and particle dynamics.

Findings

Polymer extension limit L does not affect vertical chain formation for 2-3 particles.

For 6 or more particles, smaller L prevents formation of long chains, resulting in separated clusters.

Shorter chains and smaller clusters form with decreasing L.

Abstract

In this article we present a numerical method for simulating the sedimentation of circular particles in two-dimensional channel filled with a viscoelastic fluid of FENE-CR type, which is generalized from a domain/distributed Lagrange multiplier method with a factorization approach for Oldroyd-B fluids developed in [J. Non-Newtonian Fluid Mech. 156 (2009) 95]. Numerical results suggest that the polymer extension limit L for the FENE-CR fluid has no effect on the final formation of vertical chain for the cases of two disks and three disks in two-dimensional narrow channel, at least for the values of L considered in this article; but the intermediate dynamics of particle interaction before having a vertical chain can be different for the smaller values of L when increasing the relaxation time. For the cases of six particles sedimenting in FENE-CR type viscoelastic fluid, the formation of chain of 4 to 6 disks does depend on the polymer extension limit L. For the smaller values of L, FENE-CR type viscoelastic fluid can not bring them together like the case of these particles settling in a vertical chain formation in Oldroyd-B fluid; but two separated chains of three disks are formed. Similar results for the case of ten disks are also obtained. The numerical results of several more particle cases suggest that for smaller values of L, the length of the vertical chain is shorter and the size of cluster is smaller.