Please use this identifier to cite or link to this item:

TitleA closure model for the drag coefficient of a sphere translating in a viscoelastic fluid
Author(s)Faroughi, S. A.
Fernandes, C.
Nóbrega, J. M.
McKinley, G. H.
KeywordsConfined and unconfined flow
Drag coefficient
Faxén model
Inertial effects
Oldroyd-B model
Square channel
Viscoelastic fluids
Wake effects
Wall correction factor
Wall effects
Issue dateMar-2020
JournalJournal of Non-Newtonian Fluid Mechanics
Abstract(s)In many large-scale industrial applications dealing with particle-laden viscoelastic fluids, the ensemble-averaged behavior of the mixture is of most interest. The first step to parametrize this behavior is to develop an accurate expression to rapidly evaluate the drag coefficient over a broad range of kinematic parameters. The drag coefficient of a spherical particle translating in a viscoelastic matrix is strongly affected by the viscoelasticity of the fluid. In this study, we aim to parametrize the effects of fluid elasticity, especially the relaxation and retardation times, as well as inertia on the drag coefficient of a sphere translating in a viscoelastic fluid described by the Oldroyd-B model. To this end, we employed three-dimensional direct numerical simulations of viscoelastic flow past a stationary sphere. The accuracy of the numerical formulation is thoroughly tested against a number of benchmark problems consisting of steady flow past a sphere in a bounded circular or square domain filled with either a Newtonian or viscoelastic fluid. Initially, the numerical computations for the drag coefficient over a wide range of geometric and flow parameters are validated by comparison with existing data and drag correction models from the literature. The drag coefficient correction is then evaluated for unconfined flow past a sphere at different Reynolds number, Re, over a wide range of Deborah number, De < 9, and polymer viscosity ratio, 0 < ζ < 1. For small Deborah number (De < 1), the drag coefficient decreases with respect to the Stokes drag coefficient, whereas, at large Deborah number (De > 1), the drag is enhanced due to the large elastic stresses that develop on both the surface and wake of the sphere. These canonical behaviors, observed in the inertia-less flow regime (Re ≤ 1) are amplified as the polymer viscosity ratio approaches unity. At higher Reynolds numbers (Re > 1), the drag coefficient correction arising from viscoelasticity is found to be always bigger than unity, but sm
Publisher version
AccessRestricted access (UMinho)
Appears in Collections:IPC - Artigos em revistas científicas internacionais com arbitragem

Files in This Item:
File Description SizeFormat 
  Restricted access
3,68 MBAdobe PDFView/Open

Partilhe no FacebookPartilhe no TwitterPartilhe no DeliciousPartilhe no LinkedInPartilhe no DiggAdicionar ao Google BookmarksPartilhe no MySpacePartilhe no Orkut
Exporte no formato BibTex mendeley Exporte no formato Endnote Adicione ao seu ORCID