López Carranza S N, Jenny M, Nouar C
LEMTA, UMR 7563 (CNRS), 2, avenue de la forêt de Haye, TSA 60604, 54518 Vandoeuvre Cedex, France.
Phys Rev E Stat Nonlin Soft Matter Phys. 2013 Aug;88(2):023005. doi: 10.1103/PhysRevE.88.023005. Epub 2013 Aug 8.
This study is motivated by recent experimental results dealing with the transition to turbulence in a pipe flow of shear-thinning fluids, where a streaky flow with an azimuthal wave number n=1 is observed in the transitional regime. Here, a linear stability analysis of pipe flow of shear-thinning fluids modulated azimuthally by finite amplitude streaks is performed. The shear-thinning behavior of the fluid is described by the Carreau model. The streaky base flows considered are obtained from two-dimensional direct numerical simulation using finite amplitude longitudinal rolls as the initial condition and by extracting the velocity field at time t(max), where the amplitude of the streaks reaches its maximum, denoted by A(max). It is found that the amplitude A(max) increases with increasing Reynolds number as well as with increasing amplitude E(0) of the initial longitudinal rolls. For sufficiently large streaks amplitude, streamwise velocity profiles develop inflection points, leading to instabilities. Depending on the threshold amplitude A(c), two different modes may trigger the instability of the streaks. If A(c) exceeds approximately 41.5% of the centerline velocity, the instability mode is located near the axis of the pipe, i.e., it is a "center mode." For weaker amplitude A(c), the instability mode is located near the pipe wall, in the region of highest wall normal shear, i.e., it is a "wall mode." The threshold amplitude A(c) decreases with increasing shear-thinning effects. The energy equation analysis indicates that (i) wall modes are driven mainly by the work of the Reynolds stress against the wall normal shear and (ii) for center modes, the contribution of the normal wall shear remains dominant; however, it is noted that the contribution of the Reynolds stress against the azimuthal shear increases with increasing shear-thinning effects.
本研究受近期关于剪切变稀流体管流中向湍流转变的实验结果推动,在过渡区观察到具有方位波数(n = 1)的条纹流。在此,对由有限振幅条纹进行方位调制的剪切变稀流体管流进行了线性稳定性分析。流体的剪切变稀行为由卡罗厄模型描述。所考虑的条纹基流是通过以有限振幅纵向滚动作为初始条件进行二维直接数值模拟,并在时间(t_{max})提取速度场得到的,此时条纹振幅达到最大值,记为(A_{max})。研究发现,(A_{max})随雷诺数的增加以及初始纵向滚动振幅(E_0)的增加而增大。对于足够大的条纹振幅,流向速度剖面会出现拐点,从而导致不稳定性。根据阈值振幅(A_c),两种不同模式可能引发条纹的不稳定性。如果(A_c)超过中心线速度的约(41.5%),不稳定性模式位于管轴附近,即它是一种“中心模式”。对于较弱的振幅(A_c),不稳定性模式位于管壁附近,在壁面法向剪切最高的区域,即它是一种“壁面模式”。阈值振幅(A_c)随剪切变稀效应的增加而减小。能量方程分析表明:(i)壁面模式主要由雷诺应力对壁面法向剪切的作用驱动;(ii)对于中心模式,壁面法向剪切的贡献仍然占主导;然而,需要注意的是,雷诺应力对方位剪切的贡献随剪切变稀效应的增加而增大。