环形容器内双扩散对流数值模拟_石科峰.pdf

第22 卷第2 期 中 国科学院研究生院学报 Vol.22 No.2 2005 年3 月 Journal of the Graduate School of the Chinese Academy of Sciences March　2005 文章编号:1002-11 5(2005)02-0152-05 环形容器 双扩散对流的数值模拟＊ 石科峰　卢文强 (中国科学院研究生院物理科学学院, 北京 100049) (2004 年6 月9 日收稿;2004 年 月9 日收修改稿) Shi KF, Lu WQ.Double diffusive convection in an annular vertical cavity.J ournal of the Graduate School of the Chinese Academy of Sciences, 2005, 22(2):152～ 156 摘　要　通过数值模拟, 研究了温度和浓度梯度为水平方向、 壁和外壁上的温度和浓度保持 定值的竖直环形容器 的双扩散对流结构, 其中浮升力之比N =Gr Gr =-1, Gr 和Gr 分 S T S T 别为溶质格拉晓夫数和热格拉晓夫数.首先对环形容器与矩形容器 的双扩散对流状况进行 了比较;然后着重分析了在双扩散对流情况下竖直环形容器 壁面处的热(溶质)边界层以及 平均传热(质)率Nu (Sh)随Le 数的变化情况.结果表明, 溶质边界层的厚度随着Le 数的增大 逐渐减小, 而热边界层随Le 数的变化不大;与Sh 数一直随Le 数的增大呈明显的上升趋势不 同, 当Le10 后Nu 数随Le 数的变化趋于平缓. 关键词　双扩散对流, 晶体生长, 边界层, 有限单元法 中图分类号　O551 1　引言 双扩散对流是由于流体当中热量和溶质具有不同扩散率而产生的一种复杂的对流传热传质过程, [1] 这种物理现象存在于很多领域当中, 如海洋学、地质学以及晶体生长 .在多组元晶体生长过程当中, 液 相区 的热量和溶质的传输是通过扩散和对流来完成的.生长界面的推移所产生的温度梯度和浓度梯 度决定着热量和溶质的扩散传输过程;流体中由于温度梯度和浓度梯度的存在而产生的浮升力则决定 [2] 着热量和溶质的对流传输过程 .这两种传输过程联合作用下流体的流动状况就是一种双扩散对流, 它 决定着相变界面处液相区的温度和浓度分布, 是能否获得高品质的单晶的关键.因此研究双扩散对流, 能够为以后指导制备高品质的单晶奠定基础. 目前, 一些国外学者已经开展了双扩散对流的研究, 但大多数文献均是针对矩形容器 的双扩散对 流, 如文献[3 ～ ]数值模拟了不同边界条件下矩形容器 的双扩散对流结构及其波动情况, 文献[8 ～ 10]则通过实验研究分析了矩形容器 多组元混合物凝固过程中的双扩散对流现象.与此同时, 针对环 形和圆柱形容器 的双扩散对流的研究却比较少, M.Marcoux[ 11] 仅仅对竖直环形多孔介质中的双扩散 对流进行了探讨.本文通过有限单元法求解圆柱坐标系下的涡量、流函数、能量和组分控制方程, 数值模 拟了竖直环形容器 的双扩散对流结构, 着重探讨和分析了竖直环形容器 壁面处热(溶质)边界层和 平均传热(质)率随Le 数的变化情况. ＊国家自然科学基金项目资助 　第2 期 石科峰, 卢文强:环形容器 双扩散对流的数值模拟 153 2 　控制方程及数值方法 2.1　控制方程 我们选取的物理模型和坐标系如 图1 所示.环形容器的 外壁面的半 径分别为r 和r , 径向厚度L =