1.School of Mechanical Engineering, Ningxia University, Yinchuan 750021, China 2.School of Mathematics and Statistics, Ningxia University, Yinchuan 750021, China 3.Ningxia Key Laboratory of Scientific/Engineering Computing and Data analysis, Yinchuan 750021, China
Abstract:Rayleigh-Bénard (RB) convection in binary fluid mixtures, which shows rich and interesting pattern formation behavior, is a paradigm for understanding instabilities, bifurcations, self-organization with complex spatiotemporal behavior and turbulence, with many applications in atmospheric and environmental physics, astrophysics, and process technology. In this paper, by using a high-order compact finite difference method to solve the full hydrodynamic field equations, we study numerically the RB convection in binary fluid mixtures such as ethanol-water with a very weak Soret effect (separation ratio $\psi=-0.02$) in a rectangular container heated uniformly from below. The direct numerical simulations are conducted in the rectangular container with aspect ratio of $\varGamma=12$ and with four no-slip and impermeable boundaries, isothermal horizontal and perfectly insulated vertical boundaries. The bifurcation and the origin and evolution of pattern in RB convection for the considered physical parameters are studied, and the bifurcation diagram is presented. By performing two-dimensional simulations, we observe three stable states of Blinking state, localized traveling wave and stationary overturning convection (SOC) state, and discuss the transitions between them. The results show that there is a hysteresis in the transition from the Blinking state to the localized traveling wave state for the considered separation ratio, and the evolution of the oscillation frequency, convection amplitude and Nusselt number are discontinuous. Near the lower bound of the Rayleigh number range where the Blinking state exists, a asymmetric initial disturbance is the inducement for the formation of the Blinking state. Inside the range, its inducing effect is weakened, and the oscillatory instability becomes the main reason. It is further confirmed that reflections of lateral walls are responsible for the survival of the stable Blinking state. With the increase of the Rayleigh number, the critical SOC state undergoes multiple bifurcations and forms multiple SOC states with different wave numbers, and then transitions to a chaotic state. There are no stable undulation traveling wave states at both ends of the critical SOC branch. Keywords:Convection/ bifurcation/ binary fluid mixtures/ traveling wave
速度和浓度初始值均为零, 从$ x_0 = {2\varGamma}/{3} $的非对称扰动开始, 流场经长时间的演化和发展形成了稳定的Blinking状态, 其流场的时空结构如图3(a)所示. 可以看到, 随着时间的发展, 流场中左行波和右行波交替成长和衰落, 斑图结构由左行波和右行波交替控制. 当左行波逐渐加强时, 缺陷(左行波和右行波的连接处)的中心向腔体右侧移动, 右行波的成长空间被压缩, 强度逐渐减弱; 当缺陷的中心靠近右端壁时, 在壁面的反射作用下发生转向后向左侧移动; 此后右行波逐渐成长而左行波逐渐减弱, 直到缺陷中心靠近左端壁; 然后在左端壁面的反射作用下再次发生转向, 进入下一个循环, 形成了左右闪动的对流现象. 从观测点一$ (0.13\varGamma, 0.5) $处速度w随时间的变化(图3(b)) 可以看出, 伴随着对流振幅的逐渐增大流场经历了长时间($ \Delta t\approx $3500)的对传波(counterpropagating waves, CPW)状态, 当对流达到一定强度便转变为稳定的Blinking状态. 相应的功率谱密度(图3(c)) 表明Blinking状态具有两个主频, 即闪动频率$ \omega_1 $和行波频率$ \omega_2 $. 图 3$ r=1.015 $时 (a) 流场时空结构, (b) $ w(0.13\varGamma, 0.5) $的时间序列和(c)功率谱密度 Figure3. (a) Spatio-temporal structure, (b) the time series of $ w(0.13\varGamma, 0.5) $ and (c) power spectral density for $ r=1.015 $
图4给出了$ 4000\leqslant t\leqslant 4600 $时段内两个观测点处垂向速度的变化情况. 靠近左侧壁面的观测点一和靠近右侧壁面观测点二$ (0.87\varGamma, 0.5) $的速度w分别反映了左行波和右行波振幅的变化情况. 可以看到, 左行波的振幅增加的时段对应于右行波的振幅衰减的时段, 反之亦然. 在左行波振幅较大的时段, 流场主要由左行波控制, 腔体的左半边对流较强; 相应地, 在右行波振幅较大的时段, 流场主要由右行波控制, 腔体的右半边对流较强. 图 4$ r=1.015 $时, 两个观测点处垂向速度w随时间的发展 Figure4. The time series of the vertical velocity w at two monitoring points (a) $ (0.13\varGamma, 0.5) $ and (b) $ (0.87\varGamma, 0.5) $ for $ r=1.015 $
33.2.2.Blinking状态的特性 -->
3.2.2.Blinking状态的特性
前人研究结果表明, 两端壁面的反射作用在形成和保持Blinking状态中起到了一定的作用, 那么, 端壁反射作用是不是惟一原因, 有没有其他因素在起作用呢?我们知道, 前述$ r = 1.015 $时Blinking状态是从不对称初始扰动出发计算得到的, 此时我们会想到不对称初始扰动或许是形成Blinking状态的原因之一. 为了证实这一点, 在相同的Rayleigh数下, 我们从对称的扰动(即在(6)式中取$ x_0 = {\varGamma}/{2} $)开始计算, 最终获得了与非对称扰动初值时相同的Blinking状态. 但是当$ r = 1.013 $时(为Blinking状态稳定区域下端点), 结果却不同: 从非对称扰动出发可获得稳定的Blinking状态, 而对称初始扰动未能发展起来, 流场回到了传导状态. 由此可见, 在Blinking状态存在的Rayleigh数范围的边界附近, 非对称的初始条件是诱发Blinking的原因之一, 在该范围的内部, 其诱发作用减弱. 在获得稳定的Blinking状态后, 若改变两侧壁边界条件为周期性边界条件后继续计算, 流场便回到了传导状态, 这表明两端壁面反射作用是维持Blinking状态的重要原因[9,11]. 图5给出了 $ r = 1.0171 $时Blinking和LTW两种状态中, 腔体水平中心线上瞬时垂向速度、温度和浓度波的波形, 以及浓度场斑图结构. 可以看出, 两种状态浓度波与速度波和温度波之间均存在相位差, 这是导致波动向前传播的直接原因. 图5(a)和图5(c)为Blinking状态行波向左壁面传播的情形, 此时处于对流成长的中间阶段, 其对流强度依然明显较弱, 与LTW状态的强非线性对流结构(图5(b)和图5(d), 其浓度波为带尖角的台型结构, 且其振幅远大于Blinking状态的振幅)不同, 是一种弱非线性结构. 有趣的是, 由于此Rayleigh数位于Blinking状态稳定区域的上界附近, 行波成长起来后最强处可达到LTW状态的水平, 其浓度波形具有与LTW状态类似的台型结构, 但该状态极不稳定, 对流振幅迅速衰减, 之后又再次成长. 图 5$ r=1.0171 $时, Blinking状态与LTW状态流场典型波形和浓度场的比较 Figure5. Comparison of the lateral profiles and concentration fields between the Blinking and LTW states at $ r=1.0171 $: (a) The lateral profile and (c) concentration field of the Blinking state; (b) The lateral profile and (d) concentration field of the LTW state
在Blinking状态稳定的Rayleigh数范围内, 其闪动频率$ \omega_1 $和行波频率$ \omega_2 $随Rayleigh数的变化如图6所示. 为了比较, 图6(b)中同时给出了LTW状态行波频率随Rayleigh数的变化. 随着Rayleigh数的增大, 闪动频率$ \omega_1 $先减小, 在$ Ra = 1734.9\; (r\approx 1.01575) $处取得极小值, 在此之后逐渐增大并在$ Ra = 1735.8\; (r\approx 1.01628) $达到极大值, 最后再次减小. 行波频率$ \omega_2 $随Rayleigh数的增大开始缓慢递减, 在$ r = 0.015 $之后迅速减小, 行波传播速度逐渐减缓. 当$ Ra = 1737.3\; (r\approx 1.0172) $时, 行波不再左右闪动, 伴随着行波频率的再次调整, 系统从Blinking过渡到LTW. 过渡时, 行波频率的变化是不连续的, 存在明显的跳跃, 突然提升到1.5倍以上的水平 (如图6(b)). 图 6 (a)闪动频率$ \omega_1 $和(b)行波频率$ \omega_2 $随Rayleigh数的变化 Figure6. The variation of (a) blinking frequency $ \omega_1 $ and (b) oscillation frequency $ \omega_2 $ as a function of the Rayleigh number
从Blinking到LTW状态过渡时, Nusselt数也有明显的跳跃, 如图7所示. 相比于强非线性的LTW状态, 弱非线性Blinking状态的对流振幅相对较小, 对流强度较弱, 反映对流传热贡献的Nusselt数较小. 过渡时, 伴随着对流振幅的急剧增长, 对流传热Nusselt数快速增大, 流动达到强非线性对流的水平并维持下来, 便形成了跳跃. 通过大量的模拟计算和分析, 我们进一步发现从Blinking到LTW状态的过渡是不光滑的, 存在迟滞现象. 沿Blinking分支增大Rayleigh数, 从Blinking状态到LTW状态的过渡点为$ r = 1.0172 $, 而沿LTW分支减小Rayleigh数直到$ r = 1.0171 $, 系统才从LTW过渡到Blinking状态. 可以看到, 这种滞后非常微小, 仅为$ \Delta Ra\sim 10^{-1} $, 以致在实验中未能观察到[5]. 可以预测, 在充分小的分离比下, Blinking和LTW状态之间过渡的迟滞现象就会消失, 从而实现光滑过渡. 图 7 Blinking和LTW状态的$ Nu-1 $随r的变化情况. (b)为(a)中虚线标注矩形区域的局部放大 Figure7. The variation of $ Nu-1 $ of the Blinking and LTW states as a function of r. (b) Close-up view of the part delimited by the straight dashed lines depicted in (a)
23.3.LTW状态 -->
3.3.LTW状态
本研究获得的LTW状态的流场结构如图8所示, 它是一种行波前缘与壁面接触(称为Wall-attached, LTW)[24]的局部行波对流状态, 其对流区域宽度约为$ 9.0 $, 不随Rayleigh数改变而变化. 行波传播方向可能是向左也可能向右, 这依赖于流场的初始结构, 但两者Nusselt数$ Nu $、混合参数M、振幅$ |w|_{\max} $及行波频率等特征参数完全相等, 在分岔图图2中不做区分. 事实上, 将左(右)行LTW状态的解沿$ \varGamma/2 $做对称投影$ x\leftrightarrow\varGamma-x $, 即可获得相应的右(左)行LTW状态. 由于浓度场和速度场之间存在相位差, 从而导致了一个覆盖整个对流区域的大尺度时均浓度环流, 在前缘和壁面之间形成了一个与大尺度环流反向的二次环流(图8(b)和图8(c)).大尺度浓度环流引起的浓度再分配使得LTW状态稳定下来; 反向二次环流降低了行波对流的速度, 在环形实验装置中(计算中为周期性边界), 它降低了局部对流区域整体漂移的群速度. 关于LTW状态的稳定机理等前人已做了大量研究, 我们的结论与前人并无明显不同, 本文不再赘述. 图 8 LTW状态的流场结构 Figure8. Structures of the flow field of LTW state: (a) Spatio-temporal structure; (b) a large-scale concentration current; (c) a transient structure of the concentration field
23.4.SOC状态 -->
3.4.SOC状态
33.4.1.SOC$ _{12} $状态的形成 -->
3.4.1.SOC$ _{12} $状态的形成
从零场开始计算, 在$ r = 1.13 $时我们获得了临界SOC状态–SOC$ _{12} $. 图9(a)给出了腔体水平中心线上垂向速度的时空结构, 展示了对流从零场开始到发展为SOC$ _{12} $状态的整个过程, 几个关键时刻的流场瞬时结构如图9(b)所示. 流场瞬时结构为流线和温度扰动(温度值$ \theta $减去热传导状态的解, 表现了温度偏离传导状态的大小)云图的叠加图, 其中实(点)线为正(负)值, 对应涡卷逆(顺)时针旋转, 粗虚线标出了腔体中心$ x = \varGamma/2 $. 从初始零场开始, 随着对流振幅的线性增长, 大约到了$ t = 30 $时, 腔体中心位置的涡卷率先发展起来并逐渐向两侧扩散, 形成了微弱的CPW状态, 该状态的中心驻波位于腔体的中心; 伴随着对流振幅的持续增长, CPW状态持续到$ t\approx 40 $; 之后, 壁面反射波与CPW的左行和右行波在壁面附近产生碰撞和抵消, 使得两侧壁面附近的对流强度降低, 而CPW中心驻波得到了很大的发展($ t = 46 $); 壁面反射波经过碰撞后进一步向腔体中心传播, 与驻波进行叠加并在中心处湮灭, 使得腔体内对流涡卷数量减少, 在$ 42 < t < 50 $的时段内, 中心驻波的相位发生了多次转变, 如$ t = 44.8, 46 $和 $ t = 48 $时刻; 在$ t = 50 $时刻, 反射波与中心驻波竞争结束, 在腔体中心形成了含有两对对流涡卷的结构. 之后, 对流振幅快速增长, 中心处对流涡卷强度加强, 并逐渐带动其附近流体发生对流($ t = 50, 60, 70 $), 直至对流涡卷充满整个腔体($ t = 80 $), 形成了SOC$ _{12} $状态. 图 9$ r=1.13 $时流场时空发展和典型时刻的瞬时结构 Figure9. The spatio-temporal development and transient structures of the flow field at typical times for $ r=1.13 $
对流传热Nusselt数和流体混合参数M随时间的发展如图10(a)所示, 观测点一和观测点三$ (0.5\varGamma, 0.5) $处垂向速度的时间序列由图10(b)给出. 在SOC$ _{12} $状态的形成过程中, Nusselt数曲线整体上呈现先振荡增长后振荡衰减之后再快速增长的变化形式, 期间的振荡由反射波与CPW波之间竞争(有抵消也有叠加)所致; 当二者这种竞争结束后, 腔体中对流涡卷快速发展起来, 系统在短时间内达到饱和状态SOC$ _{12} $, 对应于$ t > 50 $后Nusselt数曲线的指数增长. 相应地, 对流快速增强使得流体的混掺程度升高, M值迅速减小并最终达到一固定值. 图 10$ r=1.13 $时 (a)$ Nu-1 $和M的变化及(b)观测点处垂向速度的时间序列 Figure10. The variation of (a) $ Nu-1 $, M, and (b) the vertical velocity at the monitoring points for $ r=1.13 $
保持$ r = 1.13 $, 到$ t = 358.56 $时刻, 我们得到了稳定的SOC$ _{12} $状态, 其流场结构如图11所示, 通过对比垂向速度、温度和浓度的波形发现, 三者是同相位的, 不存在相位差; 浓度场的波形呈现有尖角的特殊情形, 沿腔体的长度方向, 尖角呈现向上向下的交替分布. 尖角位于相邻对流涡卷的交界面处, 尖角向上对应左侧顺时针(左顺)和右侧逆时针(右逆)的一对涡卷之间 (图11(b)), 左顺右逆的旋转形式将腔体上部的低温高浓度组分向下运输, 使该点的浓度高于、温度低于其临近区域. 从波形图也可以看出, 当温度场内出现波峰时, 浓度场内即出现波谷, 这一特点与温度高的地方浓度反而低的现象是符合的. 从浓度场的分布来看(图11(c)), 流体得到了充分混合, 浓度均匀分布. 图 11$ r=1.13 $时SOC$ _{12} $状态的流场结构 Figure11. The structure of flow field for the SOC$ _{12} $ state at $ r=1.13 $: (a) The lateral profile on the horizontal centerline of the cavity; (b) the streamlines and the structure of the associated temperature field; (c) the structure of the concentration field
模拟发现, 沿着SOC$ _{12} $状态分支逐渐增大Rayleigh数的过程中, 系统会经历数次分岔, 并产生多条具有不同涡卷数的SOC状态. 这些SOC状态在波数-Rayleigh数空间$ (k, r) $上较大的范围内共存, 其稳定性及向混沌状态的过渡等另文讨论. 不同分离比下, SOC$ _{12} $状态Nusselt数关于Rayleigh数变化如图12所示(注意: 在同一坐标尺度下两条曲线靠得很近, $ \psi = -0.02 $曲线略位于$ \psi = -0.1 $的上方. 为便于区分, 将$ \psi = -0.02 $曲线向左平移了0.5个r单位, 其横轴坐标数值标在图上方). 可以看到, Nusselt数随Rayleigh数指数增长, 当$ r\in[1.008, 7.59] $时, SOC$ _{12} $状态是稳定的; 在此范围之上, 对流发生二次不稳定并过渡为混沌状态, 不存在与之对应的稳定UTW状态. 这与分离比$ \psi = -0.1 $下的结论不同: 此分离比下的SOC$ _{12} $状态在$ 1.064\leqslant r\leqslant 5.69 $内稳定, 在该区域的两端均存在稳定的UTW状态. 图 12 SOC$ _{12} $状态Nusselt数随Rayleigh数的变化 Figure12. The variation of $ Nu $ with r for the SOC$ _{12} $ state