环形浅液池内热流体波的本质特征-计算力学学报.PDFVIP

第26卷第1期 　 计 算 力 学 学 报 　 Vol . 2 6 ,No . 1 2009 年 2 月 Chinese Journal of Computational Mechanics February 2009 文章编号 (2009) 0 1005907 环形浅液池内热流体波的本质特征 1 1 1 1 2 石万元 , 　李友荣 , 　彭 　岚 , 　曾丹苓 , 　今石宣之 ( 1. 重庆大学 动力工程学院 重庆 400030 ;2 . 日本九州大学 先导物质化学研究所 , 日本福冈 8168580) 摘 　要 :为了了解环形浅液池内热毛细对流的特征 ,利用有限容积法对外壁受热 、内壁冷却 、厚度为 1 mm 的环形 浅液池内硅油的热毛细对流进行了三维数值模拟 。结果表明 ,当 Marangoni (Ma) 数小于临界值时 ,随着 Ma 数的 增加 , 内、外壁附近的温度梯度上升 ,稳定的二维轴对称流动的径向速率增加 ;超过临界 Ma 数后 ,流动转化为三维 振荡对流 ,形成热流体波 。沿径向运动的同时 ,伴随着热流体波的传播流体质点成对地绕顺时针和逆时针方向旋 转 。热流体波的周向传播速度较快 ,而流体质点的周向速度很小 。分析发现 ,热流体波为对数螺线形波纹 ,其传 播角为常数 ;随着 Ma 数的增大 ,传播角增大 。 关键词 :热毛细对流 ;热流体波 ;硅油 ;数值模拟 中图分类号 : T K124 　　　文献标识码 :A 065c St 硅油 。底部为固壁, 顶部为不变形的自由边 1 　引 言 界, 底部和顶部均绝热, 内、外壁分别维持恒温 Tc 当气液界面存在切向的温度梯度时 ,热毛细力 ( ) ( ) 和 T h Th Tc 。为简化起见, 假定 : 1 硅油为不 驱使流体运动 ,形成热毛细对流 ,这一现象存在于 可压缩流体, 除表面张力外, 其他物性都为常数; 薄膜喷涂[ 1 ] 、液滴蒸发[2 ] 、熔体单晶生长[3 ] 等过程 ( ) ( ) 2 流体流速较低, 流动为层流; 3 在顶部 自由 中。当界面的温差超过阈值时 ,流动呈不稳定态 , 表面考虑热毛细力的作用, 而在所有固壁满足无滑 形成振荡的热毛细对流 ,温度的波动表现为热流体 移条件 。 波 。在 Czochral ski 法硅单晶生长过程中 ,它使熔 体内物质成分分布不均匀而影响单晶体质量[49 ] 。 不稳定的热毛细对流是流动与传热相互作用的结 果 ,其物理机理较为复杂 , 已有的报道对热流体波 现象以及发生热流体波的临界条件进行了较多的 研究 ,而对不稳定流体流动本身的行为 ,特别是流 动与热流体波传播的关系等方面的研究报道则较 图 1 　物理模型 少 。本文通过数值模拟深度为 1 mm 的环形硅油 Fig . 1 　Phy sical model 液池内热毛细对流 ,计算了流体质点的运动轨迹 , μ