In gas-liquid mass transfer processes,interfacial turbulence may occur due to the surface tension gradient and the density gradient produced by mass transfer near the interface.The interfacial turbulence can enhance t...In gas-liquid mass transfer processes,interfacial turbulence may occur due to the surface tension gradient and the density gradient produced by mass transfer near the interface.The interfacial turbulence can enhance the mass transfer since it intensifies the movement of interfacial fluid.By means of the shadowgraph optical method,the interfacial turbulence patterns vertical to the interface were observed directly in the volatilization process of binary systems.The images of the amplified interfacial turbulence showed the variation of concentration and the fluid movement under the interface.Two patterns of interfacial turbulence were observed in the experiments:plume and vortex.With the plume,the interfacial fluid moved slowly and penetrated the liquid deeply.With the vortex,the interfacial turbulence occurred in the vicinity of the liquid interface and the fluid moves quite fast.A qualitative analysis was carried out based on the mechanism of Rayleigh-Bénard convection induced by density gradient and Marangoni convection induced by surface tension gradient.展开更多
本文采用湍流热对流的并行直接数值模拟(Parallel Direct Method of Direct Numerical Simulation,PDM-DNS),计算了系列Ra数的二维方腔Rayleigh-Bénard热对流.选取典型Ra=10^(10),Pr=4.3,讨论了流动平均场的近底板区域温度分布.温...本文采用湍流热对流的并行直接数值模拟(Parallel Direct Method of Direct Numerical Simulation,PDM-DNS),计算了系列Ra数的二维方腔Rayleigh-Bénard热对流.选取典型Ra=10^(10),Pr=4.3,讨论了流动平均场的近底板区域温度分布.温度边界层沿横向可分为三个区域:羽流发射区、冲击区和大尺度环流剪切区.在羽流发射区温度剖面存在对数律,冲击区和大尺度环流剪切区没有明显的对数律特征.横向平均得到的温度剖面也存在对数律.研究系列Ra数的横向平均温度剖面,得到温度变化满足A×log(y)+B的关系,温度剖面对数律振幅–A和Ra数呈现标度关系–A^Ra^(–0.145),与实验中发现的关系基本一致.展开更多
This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7&#...This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.展开更多
提出二维湍流热对流DNS模拟的并行直接求解方法(Parallel Direct Method of DNS,PDMDNS),在"天河二号"超级计算机上实现高Ra和极高Ra湍流热对流大规模DNS计算。高分辨率的湍流热对流计算结果表明不同Ra(10~8≤Ra≤10^(13))的...提出二维湍流热对流DNS模拟的并行直接求解方法(Parallel Direct Method of DNS,PDMDNS),在"天河二号"超级计算机上实现高Ra和极高Ra湍流热对流大规模DNS计算。高分辨率的湍流热对流计算结果表明不同Ra(10~8≤Ra≤10^(13))的瞬时温度场的流场特性完全不同。较低Ra流场中有明显的大尺度环流和角涡;较高Ra流场中羽流运动充满随机性;更高Ra流场出现小尺寸漩涡并不断从上下底板产生,这些涡相互影响作用,随大尺度环流一起作绕行运动。二维湍流热对流的Nu与Ra存在标度关系,标度律约为0.3。展开更多
In this paper, we apply a scaling analysis of the maximum of the probability density function（pdf） of velocity increments, i.e., max（） = max（）up p u, for a velocity field of turbulent Rayleigh-Bénard convec...In this paper, we apply a scaling analysis of the maximum of the probability density function（pdf） of velocity increments, i.e., max（） = max（）up p u, for a velocity field of turbulent Rayleigh-Bénard convection obtained at the Taylor-microscale Reynolds number Re60. The scaling exponent is comparable with that of the first-order velocity structure function, （1）, in which the large-scale effect might be constrained, showing the background fluctuations of the velocity field. It is found that the integral time T（x/ D） scales as T（x/ D）（x/ D）, with a scaling exponent =0.25 0.01, suggesting the large-scale inhomogeneity of the flow. Moreover, the pdf scaling exponent （x, z） is strongly inhomogeneous in the x（horizontal） direction. The vertical-direction-averaged pdf scaling exponent （x） obeys a logarithm law with respect to x, the distance from the cell sidewall, with a scaling exponent 0.22 within the velocity boundary layer and 0.28 near the cell sidewall. In the cell＇s central region, （x, z） fluctuates around 0.37, which agrees well with （1） obtained in high-Reynolds-number turbulent flows, implying the same intermittent correction. Moreover, the length of the inertial range represented in decade（）IT x is found to be linearly increasing with the wall distance x with an exponent 0.65 0.05.展开更多
文摘In gas-liquid mass transfer processes,interfacial turbulence may occur due to the surface tension gradient and the density gradient produced by mass transfer near the interface.The interfacial turbulence can enhance the mass transfer since it intensifies the movement of interfacial fluid.By means of the shadowgraph optical method,the interfacial turbulence patterns vertical to the interface were observed directly in the volatilization process of binary systems.The images of the amplified interfacial turbulence showed the variation of concentration and the fluid movement under the interface.Two patterns of interfacial turbulence were observed in the experiments:plume and vortex.With the plume,the interfacial fluid moved slowly and penetrated the liquid deeply.With the vortex,the interfacial turbulence occurred in the vicinity of the liquid interface and the fluid moves quite fast.A qualitative analysis was carried out based on the mechanism of Rayleigh-Bénard convection induced by density gradient and Marangoni convection induced by surface tension gradient.
文摘本文采用湍流热对流的并行直接数值模拟(Parallel Direct Method of Direct Numerical Simulation,PDM-DNS),计算了系列Ra数的二维方腔Rayleigh-Bénard热对流.选取典型Ra=10^(10),Pr=4.3,讨论了流动平均场的近底板区域温度分布.温度边界层沿横向可分为三个区域:羽流发射区、冲击区和大尺度环流剪切区.在羽流发射区温度剖面存在对数律,冲击区和大尺度环流剪切区没有明显的对数律特征.横向平均得到的温度剖面也存在对数律.研究系列Ra数的横向平均温度剖面,得到温度变化满足A×log(y)+B的关系,温度剖面对数律振幅–A和Ra数呈现标度关系–A^Ra^(–0.145),与实验中发现的关系基本一致.
文摘This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.
文摘提出二维湍流热对流DNS模拟的并行直接求解方法(Parallel Direct Method of DNS,PDMDNS),在"天河二号"超级计算机上实现高Ra和极高Ra湍流热对流大规模DNS计算。高分辨率的湍流热对流计算结果表明不同Ra(10~8≤Ra≤10^(13))的瞬时温度场的流场特性完全不同。较低Ra流场中有明显的大尺度环流和角涡;较高Ra流场中羽流运动充满随机性;更高Ra流场出现小尺寸漩涡并不断从上下底板产生,这些涡相互影响作用,随大尺度环流一起作绕行运动。二维湍流热对流的Nu与Ra存在标度关系,标度律约为0.3。
基金supported by the Natural Science Foundation of China(Grant Nos.11102114,11202122 and 11222222)the Innovation Program of Shanghai Municipal Education Commission(Grant No.13YZ008,13YZ124)+1 种基金the Shanghai Shuguang Project(Grant No.13SG40)the Program for New Century Excellent Talents in University(Grant No.NCET-13-0)
文摘In this paper, we apply a scaling analysis of the maximum of the probability density function（pdf） of velocity increments, i.e., max（） = max（）up p u, for a velocity field of turbulent Rayleigh-Bénard convection obtained at the Taylor-microscale Reynolds number Re60. The scaling exponent is comparable with that of the first-order velocity structure function, （1）, in which the large-scale effect might be constrained, showing the background fluctuations of the velocity field. It is found that the integral time T（x/ D） scales as T（x/ D）（x/ D）, with a scaling exponent =0.25 0.01, suggesting the large-scale inhomogeneity of the flow. Moreover, the pdf scaling exponent （x, z） is strongly inhomogeneous in the x（horizontal） direction. The vertical-direction-averaged pdf scaling exponent （x） obeys a logarithm law with respect to x, the distance from the cell sidewall, with a scaling exponent 0.22 within the velocity boundary layer and 0.28 near the cell sidewall. In the cell＇s central region, （x, z） fluctuates around 0.37, which agrees well with （1） obtained in high-Reynolds-number turbulent flows, implying the same intermittent correction. Moreover, the length of the inertial range represented in decade（）IT x is found to be linearly increasing with the wall distance x with an exponent 0.65 0.05.
