环形液池内耦合热溶质毛细对流的数值模拟.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：113519 环形液池内耦合热-溶质毛细对流的数值模拟 周永利，唐经文，李友荣，龚振兴 （重庆大学动力工程学院低品位能源利用技术及系统教育部重点实验室，重庆400044 ） （Tel: 023，E-mail: liyourong@cqu.edu.cn ） 摘要：在微重力条件下，通过对环形容器内的耦合热-溶质毛细对流的二维数值模拟，得到了环形液池 内耦合热-溶质毛细对流失稳的临界条件，分析了环形容器内几何尺寸对耦合热-溶质毛细对流的影 响，结果表明：在热毛细力和溶质毛细力之比为1 的环形腔内耦合热-溶质毛细对流相比纯热毛细对流 更易发生。随着深宽比和半径比的增加，发生热-溶质毛细对流的临界毛细雷诺数Rec 会减小，液池内 的热-溶质毛细流动更容易变为周期振荡流动,阐述了其周期振荡的机理。 关键词：热毛细对流，溶质毛细对流，数值模拟 0 前言 耦合热-溶质毛细对流是由于流体中温度和浓度的不均匀引起表面张力不均匀而产 生的一种复杂的对流传热传质过程,当流体中含有两种及两种以上的物质时，溶质浓度在 自由面上的非均匀分布将产生由浓度梯度驱动引起的溶质毛细对流，它与温度梯度在自 由面处引起的热毛细对流相互耦合在一起，这种相互耦合的过程决定着相区内的温度和 浓度再分布，这种复杂物理现象广泛的存在于晶体生长、合金凝固、混合工质的相变传 热等过程中。 目前，一些国内外学者已经开展了热-溶质毛细对流的研究，但大多数文献均是针 对矩形容器内的热-溶质毛细对流：在热毛细力和溶质毛细力之比为-1 时的矩形腔内 Bergman[1]数值模拟发现出现对流的可能，这对后来的研究奠定了一定的基础，Chen[2]、 Zhan[3]分别针对二维和三维模型详细阐述了稳态和非稳态热-溶质毛细对流的机理，分析 了几何尺寸、Pr 数、Le 数对毛细对流的影响；Zhan 还分析了不同毛细力之比对发生毛 细对流传热传质特性的影响。考虑到自然界的实际情形，很多学者分析了更为复杂的对 流过程：Bahloul[4]分析了有质扩散 soret 效应的情形，Chen[5]考虑浮升力和热-溶毛细力 联合驱动的对流情形。以上这些研究均是针对矩形或立方体的模型，对于热-溶质毛细 流动周期振荡的机理也是定性的探讨，鲜有针对环形模型的研究，本文通过有限差分法 求解动量、能量和传质控制方程,数值模拟了环形容器内的热-溶质毛细对流流型结构， 着重探讨和分析了环形容器几何尺寸对毛细对流的影响情况。阐明了环形液池内的热- 溶质毛细流动周期振荡的机理。 1 物理数学模型 物理模型如图1 所示。环形液池内径为r ，外径为r ，液池的宽度为L=r -r ，深度 i o o i 为d ，流体为混合工质，底部为固壁，而顶部为自由界面，底部和顶部均绝热，内、外 壁分别维持恒温T 和T (T Tc) 、恒浓度Cc 和C (C Cc) 。因此，在表面张力梯度的 C H H H H 驱动下，在熔体内会诱发热毛细对流和溶质毛细对流。为简化起见，假定：（1）熔体为 不可压缩的牛顿混合流体，除表面张力外，满足Boussinesp 假设；（2 ）在自由表面考虑 热毛细力和溶质毛细力的作用，而在其它固壁表面满足无滑移条件；（3 ）流速较低，流 动为二维层流；（4 ）表面张力是温度和浓度函数；（5 ）忽略重力的影响，只考虑单纯的 毛细力。 图1 物理模型 2