Signature of elastic turbulence of viscoelastic fluid flow in a single pore throat

When a viscoelastic fluid, such as an aqueous polymer solution, flows through a porous medium, the fluid undergoes a repetitive expansion and contraction as it passes from one pore to the next. Above a critical flow rate, the interaction between the viscoelastic nature of the polymer and the pore configuration results in spatial and temporal flow instabilities reminiscent of turbulentlike behavior, even though the Reynolds number Re << 1. To investigate whether this is caused by many repeated pore body-pore throat sequences, or simply a consequence of the converging (diverging) nature present in a single pore throat, we performed experiments using anionic hydrolyzed polyacrylamide (HPAM) in a microfluidic flow geometry representing a single pore throat. This allows the viscoelastic fluid to be characterized at increasing flow rates using microparticle image velocimetry in combination with pressure drop measurements. The key finding is that the effect, popularly known as "elastic turbulence," occurs already in a single pore throat geometry. The critical Deborah number at which the transition in rheological flow behavior from pseudoplastic (shear thinning) to dilatant (shear thickening) strongly depends on the ionic strength, the type of cation in the anionic HPAM solution, and the nature of pore configuration. The transition towards the elastic turbulence regime was found to directly correlate with an increase in normal stresses. The topology parameter, Q(f), computed from the velocity distribution, suggests that the "shear thickening" regime, where much of the elastic turbulence occurs in a single pore throat, is a consequence of viscoelastic normal stresses that cause a complex flow field. This flow field consists of extensional, shear, and rotational features around the constriction, as well as upstream and downstream of the constriction. Furthermore, this elastic turbulence regime, has high-pressure fluctuations, with a power-law decay exponent of up to vertical bar-2.1 vertical bar which is higher than the Kolmogorov value for turbulence of vertical bar-5/3 vertical bar