过冷水快速冻结过程中的相界面温度确定.PDF

航 空 学 报 Acta Aeronautica et Astronautica Sinica hkxb@ DOI: 10.7527/S1000-6893.2016.0315 过冷水快速冻结过程中的相界面温度确定 1 2,* 王天宝 ，陈民 清华大学 航天航空学院，北京 100084 摘 要：过冷液体的冻结过程中温度分布与相界面的位置可以通过求解Stefan 问题来预测。经典Stefan 问题中认为相界面 的温度为凝固点，但在快速冻结发生时由于非平衡条件的影响，相界面的温度会低于凝固点，此时如何确定相界面的温 度是一个难点。通过分子动力学模拟方法研究了过冷水快速冻结过程，得到了相界面在不同移动速度下相界面温度的变 -1 -1 化规律。相界面过冷度为0~5K范围内，界面动力学系数可以近似为0.55m·s ·K ，即相界面移动速度每增加0.55m/s ，相 界面的温度降低1K。对于更高的相界面过冷度，相界面移动速度和相界面过冷度的关系可用一个二次函数来描述。该相 界面温度确定方法可应用于过冷水快速冻结的数值模拟。 关键词：相界面温度；相界面移动速度；界面动力学系数；过冷水冻结；Stefan 问题 中图分类号：V211.3 文献标识码：A 文章编号： 飞机在穿过含有过冷水滴的云层时，过冷水 中认为相界面的温度为凝固点。而实际上，在快 滴会撞击在飞机表面导致结冰，这是造成飞机安 速冻结过程中，由于非平衡性的影响，相界面的 [1] 全性事故的重要因素之一 。过冷水滴撞击在飞 温度与相界面的移动速度有关。相界面移动速度 机表面上的结冰是一个碰撞、变形及冻结的复杂 与相界面过冷度的比值被定义为界面动力学系数 过程。从冻结的角度来看，传统的数值模拟方法 [4] 。近年来的结冰数值模拟研究中已经引入了相 [5-8] 通常假设相界面的温度恒等于环境条件下的凝固 界面移动速度来修正相界面的温度 ，但其中 点。此假设即便对于静态过冷水滴的冻结过程， 界面动力学系数具体如何确定仍然有待研究。 仍有待于进一步考察。 目前分子动力学 （Molecular Dynamics ，简 静态过冷水滴的冻结过程宏观上表现为液固 称MD ）方法被广泛地应用于物理、化学、材料 相界面的移动，可以用Stefan 问题来描述。Stefan 科学、生命科学等各个领域，其中也包括固液相 在1891年研究了冰的熔化过程并建立了对应的数 变过程的研究。Nada等人用分子动力学方法研究 [2] 学模型 ，该数学模型也可用于处理过冷水冻结 了冰晶的各向异性生长，并且发现从六角冰的底 问题。此后类似的带相变的传热问题或带自由边 面生长出来的冰-水界面是光滑的，而从六角冰 [3] 界的数学问题都被称为Stefan 问题 。 棱面生长得到的冰-水界面则是粗糙的，这种粗