分离结晶过程中熔体热毛细—浮力对流研究.pdf

Research on Thermocapillary-Buoyancy Convection in Melt during Detached Solidification Process A Thesis Submitted to Chongqing University in Partial Fulfillment of the Requirement for the Doctor ’s Degree of Engineering By Li Zhen Supervised by Prof. Peng Lan Specialty：Power Engineering and Engineering Thermophysics College of Power Engineering of Chongqing University, Chongqing, China April 2013 中文摘要 摘 要 新型分离结晶Bridgman 法 （简称：分离结晶法）完美地结合了传统的Cz 提 拉法和垂直Bridgman 法的优点，成为熔体法生长热敏感性晶体材料的最佳方法。 熔体法晶体生长过程中，当自由表面上存在温度梯度时，会在熔体内部形成热毛 细对流。同时，常重力条件下由于熔体密度分布不均匀，还会导致熔体内部产生 浮力对流。热毛细对流与浮力对流的存在，均会引起晶体生长质量的下降，分离 结晶法也不例外。因此，掌握分离结晶过程中熔体内部热毛细－浮力对流的流动 规律，并对其进行主动调控成为分离结晶法生长高质量晶体材料急需解决的科学 问题。本文即以分离结晶过程中熔体流动为研究对象，采用理论分析、数值模拟 与相关实验相结合的方法对其进行了较为深入的研究，重点探究热毛细－浮力对 流的流动规律，得到了熔体内部的速度分布与温度分布，讨论了Marangoni 数、气 缝宽度、熔体高度以及轴向磁场对熔体流动的影响，确定了熔体流动失稳的临界 Marangoni 数，分析了流动失稳的物理机制。 主要研究内容及结果如下： 首先，建立了分离结晶过程中控制晶体-熔体-坩埚气液界面形状的微分方程 组，采用龙格-库塔法对该微分方程组进行了数值求解。结果表明，气液界面两侧 压差和熔体-坩埚接触角是影响气液界面形状和气缝宽度的主要因素。α +θ 180° e c 时，气液界面为一个凸向熔体侧的弯液面，气缝宽度随着气液界面两侧压差的增 大而减小，随着接触角的增大而减小；α +θ 180°时，气液界面为一个凸向气体侧 e c 的弯液面，气缝宽度随着气液界面两侧负压差的减小而增大，随着接触角的增大 而增大。当熔体冷热端气压差调节速率与结晶速率成线性变化时，能够实现稳定 的分离结晶过程。气缝宽度的计算结果，为后续进行熔体热毛细-浮力对流的数值 模拟提供了参考依据。 其次，建立了分离结晶过程中熔体热毛细-浮力对流的物理数学模型，采用有 限差分法对其进行了详细的三维数值模拟。结果表明，在热毛细力和浮力的共同 驱动下，熔体内部存在上下两个流动方向相反的流胞。影响熔体内热毛细-浮力对 流的主要因素有Marangoni 数、气缝宽度和熔体高度等。随着Marangoni 数的增大， 熔体内部流动逐渐增强，当Marangoni 数超过某一临界值后，熔体内部流动会转变 成三维非稳态流动。同时，熔体内部流动随着气缝宽度和熔体高度的减小而增强， 流动失稳的临界Marangoni 数随着气缝宽度的减小而增大，随着熔体高度的减小略 微减小。流动失稳的物理机制可以解释为熔体内温度的变化与速度的变化二者不 同步，导致熔体内流胞流动出现减弱