A viscoelastic flow solver for Rolie-Poly model based on OpenFOAM

The understanding of flow instability of viscoelastic fluid is very important but far from sufficient. By combining the constitutive equation with Navier-Stoke governing equation, we implement a viscoelastic solver based on OpenFOAM for Rolie-Poly model, which is a single-mode constitutive model derived from the tube-based molecular theory. We numerically investigate two-dimensional viscoelastic flow around a circular cylinder with this solver. A second order parallelized finite volume method (FVM) is employed in the simulation in both spatial and the temporal domain. The numerical results illustrate the significant effects of the polymer elasticity on the flow structure. The vortex shedding and recirculation region exist behind the cylinder under low Reynolds numbers due to the elastic effect. Both drag reduction and enhancement have been observed in the simulation. Results are consistent with other literature. And the solver shows good parallel efficiency and scalability.

[1]  Robert C. Armstrong,et al.  Finite element method for viscoelastic flows based on the discrete adaptive viscoelastic stress splitting and the discontinuous Galerkin method : DAVSS-G/DG , 1999 .

[2]  Ma Martien Hulsen,et al.  Simulation of the Doi-Edwards model in complex flow , 1999 .

[3]  Xuejun Yang,et al.  WI-USHER: A grid-based parallel algorithm for particle insertion in hybrid atomistic-continuum method , 2017 .

[4]  Tim Reis,et al.  Rolie-Poly fluid flowing through constrictions: Two distinct instabilities , 2013 .

[5]  Charles-Henri Bruneau,et al.  A numerical study of two dimensional flows past a bluff body for dilute polymer solutions , 2013 .

[6]  Yuichi Masubuchi,et al.  Numerical study of chain conformation on shear banding using diffusive Rolie-Poly model , 2011 .

[7]  Richard S. Graham,et al.  Simple constitutive equation for linear polymer melts derived from molecular theory: Rolie–Poly equation , 2003 .

[8]  Hrvoje Jasak,et al.  Viscoelastic Flow Simulation: Development of a Methodology of Analysis Using the Software OpenFOAM and Differential Constitutive Equations , 2009 .

[9]  F. Nieuwstadt,et al.  Visco-elastic flow past circular cylinders mounted in a channel: experimental measurements of velocity and drag , 2004 .

[10]  Aleksandar Jemcov,et al.  OpenFOAM: A C++ Library for Complex Physics Simulations , 2007 .

[11]  Paulo J. Oliveira,et al.  A numerical study of steady and unsteady viscoelastic flow past bounded cylinders , 2005 .

[12]  N. Phan-Thien,et al.  The flow of an Oldroyd-B fluid past a cylinder in a channel: adaptive viscosity vorticity (DAVSS-ω) formulation , 1999 .

[13]  Wei Yi,et al.  Toward large scale parallel computer simulation of viscoelastic fluid flow: A study of benchmark flow problems , 2015 .

[14]  Mahmood Norouzi,et al.  Numerical investigation of viscoelastic shedding flow behind a circular cylinder , 2013 .

[15]  O. Baskurt,et al.  Blood Rheology and Hemodynamics , 2003, Seminars in thrombosis and hemostasis.

[16]  M. Vannoni,et al.  Relaxation time and viscosity of fused silica glass at room temperature , 2011, The European physical journal. E, Soft matter.