Elasto-visco-plastic flows in benchmark geometries: II. Flow around a Confined Cylinder

TL;DR

This study computationally investigates elastoviscoplastic fluid flow around a cylinder, revealing complex unyielded regions, elongated yielded zones, and effects of material properties on flow behavior using the Saramito-Herschel-Bulkley model.

Contribution

It provides new insights into flow patterns and unyielded regions in EVP fluids around a cylinder, highlighting the influence of elasticity and yield stress on flow dynamics.

Findings

Unyielded regions form around the symmetry plane and near the cylinder.

Elongated yielded zones appear downstream under certain conditions.

Material elasticity affects the size of unyielded regions and flow oscillations.

Abstract

We examine computationally the two-dimensional flow of elastoviscoplastic (EVP) fluids around a cylinder symmetrically placed between two plates parallel to its axis. The Saramito-Herschel-Bulkley fluid model is solved via the finite-volume method using the OpenFOAM software. As in viscoplastic materials, unyielded regions arise around the plane of symmetry well ahead or behind the cylinder, as two small islands located above and below the cylinder and as polar caps at the two stagnation points on the cylinder. Most interestingly, under certain conditions, an elongated yielded area around the midplane is predicted downstream of the cylinder, sandwiched between two unyielded areas. This surprising result appears, for example, with Carbopol 0.1% when considering a blockage ratio of 0.5 (the ratio of the cylinder's diameter to the channel's width) and above a critical elastic modulus (G>30 Pa). An approximate semi-analytical solution in the region mentioned above, reveals that it is caused by the intense variation of the stress magnitude there, which may approach the yield stress asymptotically either from above or below, depending on material elasticity. The drag coefficient on the cylinder increases with yield stress and blockage ratio but decreases with material elasticity. The unyielded regions expand as the yield stress increases. They also expand when material elasticity increases because this allows the material to elastically deform more before yielding. Furthermore, by decreasing the elastic modulus or increasing the yield stress beyond a critical value, the yield surface may exhibit damped oscillations, or irregular shapes even without a plane of symmetry, all under creeping flow conditions.