Characterizing elastic turbulence in the three-dimensional von Karman swirling flow using the Oldroyd-B model

TL;DR

This study uses three-dimensional simulations with the Oldroyd-B model to analyze the transition to elastic turbulence in von Karman swirling flow, revealing a subcritical transition, hysteresis, and power-law scaling consistent with experiments.

Contribution

It provides the first detailed 3D numerical characterization of elastic turbulence onset and development in von Karman flow using the Oldroyd-B model, including flow transition details and spectral analysis.

Findings

Transition to elastic turbulence occurs at Wi_c=12.

Hysteresis observed in flow resistance and order parameter.

Power-law exponents match experimental values.

Abstract

We present the full three-dimensional numerical investigation of the von Karman swirling flow between two parallel plates using the Oldroyd-B model and characterize the onset and development of elastic turbulence. We quantify the flow state with the secondary-flow strength, a measure of the average strength of the velocity fluctuations, and then define an order parameter as the time average of the secondary-flow strength. The order parameter displays a subcritical transition from the laminar to a bistable flow that switches between weakly chaotic flow and elastic turbulence. The transition to the bistable flow occurs at the critical Weissenberg number $\mathrm{Wi_c}=12$. Above $\mathrm{Wi_c}$, in the elastic turbulent state, we observe a strong increase in velocity fluctuations and flow resistance, which we define as the total work performed on the fluid. Upon starting simulations in the turbulent state and subsequently lowering $\mathrm{Wi}$ below its critical value, we observe hysteretic behavior in the order parameter and the flow resistance, which is a common feature of a subcritical transition. Hysteresis has also been found in experiments. Additionally, we find power-law scaling in the spatial and temporal power spectra of the velocity fluctuations, characteristic for elastic turbulence. The maximum values of the power-law exponents in our simulations are $\alpha_t= 3.41$ for the temporal exponent and $\alpha_s= 3.17$ for the spatial exponent, which are remarkably close to the values obtained in experiments.