自由表面形状对熔体内部流动和温度分布的影响.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123444 自由表面形状对熔体内部 流动和温度分布的影响一 彭岚，周洋，李友荣，王飞 （重庆大学动力工程学院，低品位能源利用技术及系统教育部重点实验室，重庆400044 ） Tel: 023 E-mail penglan@cqu.edu.cn （ : ） 摘 要：本文借助有限体积法，在常重力和给定熔体下壁面均匀热流条件下，研究矩形液池自由表面 形状对晶体生长过程中熔体内部流动与温度分布的影响。模拟结果如下：1）最低温区域介于中心线和 2 右侧壁之间，熔体内部具有两个流向相反的流胞； ）凸形液面使熔体内部流动增强，最低温度值增大， 凸度增大使最低温区域向右侧壁方向移动；3 ）凹形液面使熔体内部流动减弱，最低温度值减小，凹度 增大使最低温区域向左侧壁方向移动。 关键词：液面形状；熔体流动；垂直热流；表面辐射；有限体积法 0前言 LiNbO3 [1]在微波技术、大规模集成光学系统等领域中具有重要用途，是常用 的晶体生长物质。本文针对晶体生长查阅文献如下所述，2003 年，由You-Rong Li 等学者[2]针对环形液池中硅熔体内部流动的不稳定性进行了一系列三维数值模拟，发现 根据工业实际情况确定的3w/cm2 的下壁面均匀热流对流动影响不大。同年，You-Rong Li, Peng L 等学者[3]对与前者主题有相关性的问题进行了研究；2004 年，Q.Kang 等学者[4]讨论了矩形腔中流体的自由表面变形；2007 年，由Peng L, You-Rong Li 等 [5] 学者 针对Cz 法硅熔体浅液层中瞬态的热对流进行了三维数值模拟，以弄清流型的瞬态 机理。2011 年，K.Li 等学者[6]在微重力条件下模拟了Marangoni 对流和热毛细对流的耦 合作用，并研究了不同液面形状下的耦合情况。本文借助fluent 软件模拟研究晶体生长 的内部流动问题，国内外学者对晶体生长方法的研究具有两个共同点：1）通常着重于横 向温度梯度对熔体流动的影响；2 ）通常研究的是假定为平面的自由表面。而实际上，熔 体内部存在竖直温度梯度，并且固液气界面的接触角和内部流动的影响使液面并非平面。 因此本文在给定下壁面均匀热流的条件下，针对不同液面形状的熔体内部流动以有限体 积法[7]进行二维数值模拟，研究工质采用典型物质LiNbO . 3 1物理数学模型 本文采用的数值模拟物理模型如图 1 所示，模型分为内部熔体，左、右等高侧壁以 及下壁面，上部为熔体自由表面，下壁面长度L=55.5mm ，侧壁高度为d=8mm 。左侧壁 温度 T =1606K ，右侧壁温度 T =1586K ，液面与环境之间辐射换热，环境温度 h c 一 基金项目：国家自然科学基金（No）资助项目 T =0.9T =1374.3K ，T 为工质熔点温度，在下壁面具有竖直方向上的均匀热流q 。考虑 a m m 到利用Gambit 软件建立模型的方便性，以及本文着重于对研究液面实际变形铺垫的前 提，因此假定液面形状为椭圆的一半，通过椭圆中心沿长轴对开两半，即为凸液面和凹 液面，改变椭圆的短半轴R 即可实现凹凸程度的变化，以此改变液面形状。 根据模拟的需要，作如下假定：(1)液体为不可压缩的牛顿流体，除了表面张力以外， 相关物性参数均为常数，满足Boussinesq 近似；(2)表面张力对温度的梯度值为常数；(3) 自由表面满足给定形状，固液界面满足无滑移边界条件；(4)流动为层流； 1