霜晶体结构对霜层热导率的影响.docxVIP

霜晶体结构对霜层热导率的影响 热导率的模型 湿空气中的蒸馏水在寒冷表面上凝固，是一个与传热传递质量和特殊多孔介质运动轨迹有关的问题。由于问题的复杂性, 许多研究工作都需要进一步深入。结霜问题的重点是霜层物理性质的研究, 其中最重要的参数之一是霜层的热导率。 霜层的热导率由结霜的环境条件所决定, 并且与结霜过程的历史密切相关, 随霜层的生长而连续变化。许多研究表明, 热导率主要取决于密度, 同时也受霜层微观结构的影响。因此, 热导率的变化相当复杂, 要得到一个基于严密理论的关于热导率的关系式是非常困难的。 许多学者都提出了霜层的结构模型, 从理论上探讨热导率的变化规律, 并进行了实验研究。然而, 大多数模型都认为霜层是一个均匀结构, 建立的模型大多将冰晶体作为柱状。但事实上, 霜层结构非常复杂, 霜层在不同温度下生长时, 会呈现树枝状、柱状以及平板状形态。因而目前的各种模型都在不同程度上存在缺陷。为此, 本文研究了霜层中树枝状冰晶体的生长过程以及微观结构, 以期能够弥补目前文献的不足, 进一步加深对霜层生长过程的认识, 并为预测霜层的生长奠定基础。 1 冰晶体的形成过程 图1 (a) 、 (b) 是冰晶体在不同的环境条件下生长时所拍摄的放大图片。从图中可以看出, 冰晶体的结构呈树枝状, 并且随着时间的增长, 不断有新的树枝状晶体出现。图1 (c) 是霜层在充分生长后拍摄的图片, 此时霜层密度比较大, 已经很难将一个霜晶体与另一个区分开来。 目前, 研究结霜现象的工作者大都认同Hayashi等将结霜过程分为3个阶段的提法, 即霜柱的发生期、霜层生长期以及霜层成熟期。首先是冰晶体在冷表面上形成霜柱, 这些霜柱是霜层构造的骨架。此后, 霜柱表面不断有新的晶体生长, 形成树枝状的霜层, 图1所显示的就是这个时期的图片。此后霜层的厚度不断加大, 表面温度持续升高, 由于水蒸气凝华放热, 此时霜层会融断或者融化, 融断的霜晶体会在重力作用下倒伏, 而融化形成的液态水会渗入霜层, 在霜层内部重新凝结为冰。同时, 水蒸气仍然不断扩散进入霜层, 这就是霜层的成熟期。对于霜层的成熟期, 由于霜层密度较大, 采用显微拍摄的方法很难获得清晰的图片。因此, 融化形成的液态水在何处发生凝结难以预计。 从以上叙述可以看出, 冰晶体的形成过程极为复杂, 并且具有阶段性。因此, 在模拟霜层的生长时, 有必要分阶段描述。 2 witen-sada模型 Witten和Sander提出了一种模型来解释烟尘微粒的聚集。这个模型的规则很简单:先在一个平面的中心放置一个粒子, 称其为种子。然后在远离种子的地方释放一个粒子, 该粒子可以随机行走 (扩散) , 当运动到种子旁边时, 停止行走, 留在该处不动。接着释放第二个粒子, 规则与第一个相同。如此反复进行, 就可以获得一个大的凝聚体, 该凝聚体的结构具有有限扩散条件下生成物的特征, 即一个树枝状的标度不变的复杂结构。这个树枝状的机构具有统计意义下的相似性。这个树枝状结构的最尖端可以最有效地捕获扩散过来的粒子。这样, 使得小的涨落被强化。这个不稳定性与模型固有的随机性结合起来, 形成了一个具有综合特性的凝聚体。图2是30000个粒子采用Witten-Sander模型模拟得到的结果。 从图2可以看出, Witten-Sander得到的结果具有多个分支, 其中每个分支都与实验图片的冰晶体非常相似。对Witten-Sander方法适当改进, 便有可能得到与实验结果吻合的霜柱发生期以及霜层生长期的冰晶体形状。 为了模拟结霜过程中的上述两个时期, 对Witten-Sander模型提出以下改变。首先, 选取单位长度的一个正方形。然后将其n等分, 这样这个正方形具有n×n个网格。使用0来代表空气网格点, 1代表霜晶体网格点。选取最底层的一行网格将其标记为1, 这行网格就是种子, 其余网格标记为0。然后在正方形最上面一行随机选取一个作为起始点, 释放一个粒子。这个粒子随机行走, 遇到种子的时候停止, 此时粒子所在的网格标记为1。然后继续释放下一个粒子, 直到这个粒子遇到标记为1的网格。这样依次进行, 最后形成树枝状的结构。 图3 (a) 是采用上述方法模拟得到的霜层结构。从图中可以看出, 释放360个粒子以及540个粒子所得到的图片与实验所得的图片吻合非常好。而释放720个粒子所得的图片在网格顶部发生了霜晶体的聚集, 与实验存在较大的差别。这说明此时霜层已经进入成熟期, 需要采用另外的方法。 在进入霜层成熟期后, 霜层表面会发生融化, 此过程形成的液态水会渗入霜层内部而凝固。至于在什么地方凝固, 很难准确预计。而早期生成的霜柱也可能会融断倒伏, 对于这个现象同样难以准确预测。因此假定在网格没有霜晶体的部分, 都有可能有霜晶体的加入, 并且