Mechanism of elastic instability in Couette flow of polymer solutions: Experiment

Abstract

Experiments on flow stability and pattern formation in Couette flow between two cylinders with highly elastic polymer solutions are reported. It is found that the flow instabilities are determined by the elastic Deborah number, De, and the polymer concentration only, while the Reynolds number becomes completely irrelevant. A mechanism of such “purely elastic” instability was suggested a few years ago by Larson, Shaqfeh, and Muller [J. Fluid Mech. 218, 573 (1990)], referred to as LMS. It is based on the Oldroyd-B rheological model and implies a certain functional relation between De at the instability threshold and the polymer contribution to the solution viscosity, ${\eta }_p\mathrm{/\eta ,}$ that depends on the polymer concentration. The elastic force driving the instability arises when perturbative elongational flow in radial direction is coupled to the strong primary azimuthal shear. This force is provided by the “hoop stress” that develops due to stretching of the polymer molecules along the curved streamlines. It is found experimentally that the elastic instability leads to a strongly nonlinear flow transition. Therefore, the linear consideration by LMS is expanded to include finite amplitude velocity perturbations. It is shown that the nature of the elastic force implies major asymmetry between inflow and outflow in finite amplitude secondary flows. This special feature is indeed exhibited by the experimentally observed flow patterns. For one of the flow patterns it is also shown that the suggested elastic force should be quite efficient in driving it, which is important evidence for the validity of the mechanism proposed by LMS. Further, the predicted relation between De and ${\eta }_p\mathrm{/\eta }$ is tested. At fixed ${\eta }_p\mathrm{/\eta }$ the elastic instability is found to occur at constant Deborah number in a broad range of the solution relaxation times in full agreement with the theoretical prediction. The experimentally found dependence of the Deborah number on ${\eta }_p\mathrm{/\eta }$ also agrees with the theoretical prediction rather well if a proper correction for the shear thinning is made. This provides further support to the proposed instability mechanism.