Spatio temporal linear stability of viscoelastic free shear flows: dilute regime

Sarthok Sircar,Diksha Bansal

DOI: https://doi.org/10.1063/1.5115455

2019-06-20

Abstract:We report the temporal and spatio-temporal stability analyses of anti-symmetric, free shear, viscoelastic flows obeying the Oldroyd-B constitutive equation in the limit of low to moderate Reynolds number and Weissenberg number. The resulting fourth order Orr-Sommerfeld equation is reduced to a set of six auxiliary equations which are numerically integrated starting from the rescaled far-field conditions, i.~e., via the Compound Matrix Method. Numerical results indicate that with increasing Weissenberg number: (a) the peak of the maximal growth rate (i. e., the maximum value of the imaginary component of the growth rate) is reduced, (b) the entire range of unstable spectrum is shifted towards longer waves (i.~e., the entire region of temporal instability is gradually concentrated near zero wavenumber), (c) the vorticity structure contours are dilated and (d) the residual Reynolds stresses are diminished. All these observations suggest that viscoelasticity reduces the temporal instability but does not completely suppress it. The Briggs idea of analytic continuation is deployed to classify regions of temporal stability, absolute and convective instabilties as well as evanescent (false) modes, in a finite range of Reynolds number, Weissenberg number and the viscosity coefficient. The main result is that, unlike Newtonian fluids, the free shear flow of dilute polymeric liquids are either (absolutely/convectively) unstable for all Reynolds number or the transition to instability occurs at very low Reynolds number, a finding attributed to the fact that viscoelasticity aggravates free surface flow instabilities. Although this transitional pathway connecting the temporally stable state to the elastoinertially unstable state have been identified by some in-vitro experiments, but until now, has not been quantified theoretically via linear stability analysis.

Fluid Dynamics,Plasma Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the spatio - temporal stability of free - shear flows of dilute polymer solutions and their transition mechanisms from laminar to turbulent flows. Specifically, the research objectives include:

1. **Determine the critical conditions of instability**: Through linear stability analysis, find the critical flow/polymer rheological conditions that trigger flow instability.
2. **Explain the effect of elasticity on flow stability**: Explore the influence of viscoelasticity on flow stability and instability under different Reynolds number (Re) and Weissenberg number (We) conditions.
3. **Reveal the transition path**: Provide a theoretical linear spatio - temporal stability analysis to explain the transition path from the temporally stable state to the elastically - inertial unstable state, especially for free - surface flow instabilities at low Reynolds numbers.

### Main problems and methods

The author uses the Oldroyd - B constitutive equation to describe the viscoelastic behavior of dilute polymer solutions and conducts research through the following steps:

- **Mathematical model establishment**: Based on the Navier - Stokes equations and the evolution equations of the extra - stress tensor, a mathematical model describing free - shear flows is constructed.
- **Linear stability analysis**: Linearize the governing equations to obtain the Orr - Sommerfeld equation (OSE), and solve this equation by the Compound Matrix Method.
- **Spatio - temporal stability analysis**: Classify the absolute instability, convective instability and evanescent modes through the Briggs analytic continuation method, and draw the stability phase diagrams in the material parameter space.

### Key findings

1. **Peak decrease of the maximum growth rate**: As the Weissenberg number increases, the peak of the maximum growth rate decreases, and the range of the unstable wave spectrum shifts towards longer wavelengths.
2. **Change in vorticity structure**: The contour of the vorticity structure expands and the residual Reynolds stress weakens, indicating that viscoelasticity reduces temporal instability.
3. **Absolute and convective instabilities**: Unlike Newtonian fluids, the free - shear flows of dilute polymer liquids can be absolutely or convectively unstable at all Reynolds numbers, or undergo instability transitions at very low Reynolds numbers.
4. **Classification of spatio - temporal instabilities**: Through spatio - temporal stability analysis, the regions of absolute and convective instabilities under different parameters are successfully classified, and the existence of spurious modes is revealed.

These results are helpful for in - depth understanding of the behavior of viscoelastic fluids in fields such as microfluidic mixing and bio - fluids, and provide theoretical support for industrial applications.