Effects of Shear Intensity on the Linear Instability of Viscoelastic Rayleigh-Bénard-Poiseuille Flow

Zhen-Ze Yao,Cai-Lei Lu,Chu -Tong Zhou,Kang Luo,Hong-Liang Yi,He -Ping Tan
DOI: https://doi.org/10.1016/j.ijheatfluidflow.2024.109336
IF: 2.643
2024-01-01
International Journal of Heat and Fluid Flow
Abstract:We study the two-dimensional Rayleigh-Bénard-Poiseuille flow in a viscoelastic fluid using linear stability analysis and energy analysis. Different from the case of Newtonian fluids, three instability regimes are identified when gradually increasing the shear intensity (characterized by the Reynolds number Re). In regime I with a very small Re, the viscoelasticity has destabilizing effects on the flow, leading to a smaller critical Rayleigh number Rac than that in the Newtonian case. This is due to the generation of oscillatory instability for high elasticity. Energy analysis indicates that PD (the rate of work done by the polymeric stresses) makes a negligible contribution to the total energy balance and viscous dissipation (VD) becomes the main driving force to reduce the perturbation kinetic energy. In regime II with moderate Re, in contrast to regime I, the viscoelasticity has stabilizing effects on the flow, and this results in a higher Rac in comparison to the Newtonian case. This can be evidenced by the increase of the negative work rate of the polymeric stresses. In higher Re regime III, the viscoelasticity is again enhancing flow instability, and this is owing to the increase of PD. Among all polymer stress components, the streamwise normal stress plays the most important role. In addition, the effects of Weissenberg number W, viscosity ratio β and Prandtl number Pr on linear instability threshold have been investigated. By fixing Ra and streamwise wavenumber α, energy analysis indicates that suppressing the heat energy transferred from perturbed temperature field is the main way to inhibit flow instability at lower Re, while the disturbance kinetic energy absorbed by shear flow is the major stabilizing factor at larger Re. In very high Re or W cases, Tollmien–Schlichting instability induced by strong shear and center mode instability caused by strong elasticity are dominant respectively, which are nearly unaffected by the temperature gradient.
What problem does this paper attempt to address?