基于蒙特卡洛重构模型的多孔材料热导率计算.docxVIP

基于蒙特卡洛重构模型的多孔材料热导率计算 多孔材料是自然界广泛存在的多相混合体。它在航空、航空航天、冶金、化工、木材、机械等领域占有重要地位。热导率是多孔材料的一个重要热物理参数，它代表了多孔材料的隔热性能。获得正确的热导率在多孔材料的使用和设计中起着重要的领导作用。然而，由于多孔材料内部结构的复杂性以及组成材料的离散性，其内部加热过程中的温度存在许多奇怪现象。因此，传统的传统模式很难准确地描述其内部加热过程。多孔材料的热导率是许多科学家在热容量方面进行研究的，其研究方法主要分为实验测量法和理论分析法。实验测量法测量实际材料样品的热导率所需的实验设备等条件，因此有一定的局限性。因此，许多科学家使用理论方法讨论了多孔材料的热导率。程远贵等人使用分布式理论和等统一热阻法等共同方法，推断出在高温下抗火纤维材料的等效热导率。于子涛等人使用热阻模拟方法，人为木材的横向有效导数系数的公式。一系列多孔材料的内部热阻模型被划分为五种不同的对应热阻模型，并导出了不同模型的热阻率计算公式。这些模型是多孔材料内部结构的简化，部分结构的传统模式可以取代一般的传播模型。根据多孔材料的内部结构特点，这种设计方法可以简化不同类型材料的内部热阻力，并根据不同类型材料的热阻力计算所需的二值化矩阵法。这种方法具有更高的适应性。 1 计算多孔材料的热导率的模型 1.1 固体基质的对流传热 多孔材料内部传热过程主要包括:(1)固体基质间的导热;(2)孔隙中流体之间的导热;(3)孔隙中流体与固体基质之间的对流换热;(4)固体之间及固体与流体之间的辐射换热. 本文研究的多孔材料内部流体主要是气体. 对于本文所关注的多孔介质,其孔隙直径都小于5 mm. 已有研究表明在孔隙当量直径小于5 mm时,对流换热所占比例非常小,可忽略不计. 因此,本文中的热导率计算模型均不考虑对流换热的影响. 1.2 热阻阵列的计算 对于一块多孔材料试样,假设一维热流从试样一个壁面传到相对壁面,保持这两个壁面为恒温,其余四壁绝热. 忽略内部对流换热,可以通过下式计算其内部等效热阻: R=ΔTQR=ΔΤQ. (1) 式中:R为试样的等效热阻,K·W-1;Q为传热方向上的热流量,W;ΔT为热流入口与出口壁面的温差,K. 由此可以得到该试样的整体热导率: k=LARk=LAR. (2) 式中:k为试样的热导率,W·m-1·K-1;L为试样沿传热方向的长度,mm;A为试样沿传热方向的截面积,mm2. 由于忽略试样内部的对流换热,所以试样内部只存在三种换热过程,即固体骨架单元之间的导热换热、孔隙内气隙单元之间的导热换热和固体骨架与孔隙气隙间的辐射换热. 任意取一块多孔介质模型,对其进行网格划分,如图1所示. 进一步进行二值化处理,对模型划分后的网格进行判断,当格子内黑色像素点所占的面积(或体积)超过整个盒子的一半时,则把此盒子看成是孔隙,否则看成是固体基质. 因此模型内部存在两种状态:黑色(代表孔隙)或者白色(代表固体基质). 经过二值化处理和判断后,多孔材料内部孔隙和固体基质的状态可以看成是如图2所示的阵列图. 利用这个阵列图,可以对多孔材料的热导率进行计算,将阵列图中的黑色和白色盒子看成对应气体和基质的等效热阻,利用阵列图中各盒子间的串并联关系,可对其热导率进行计算. 如果所计算的多孔材料内部结构至少在一个方向上有着相同的宏观传热特性,当一维热流垂直于此方向通过多孔材料时,则垂直于热流方向各片材料的热阻阵列为串联关系. 对于各片材料,其气体和基质的等效热阻关系如图3(a)所示,其中白色单元为二值化后判断为代表基质的单元,黑色单元为判断为气体部分的区域单元,单元间连接线代表两个相邻单元之间的传热关系. 假设单片材料在某一方向上(如图中的纵向)有着相同的宏观传热特性,热流方向垂直于纸面,如图3(b)所示,则气体和基质的等效热阻在传热方向为并联关系,在非传热方向为串联关系. 如果单片材料在垂直热流的两个方向上都有着相同的宏观传热特性,则此片材料中气体和基质的等效热阻都为并联关系,如图3(c)所示. 由此,图2中的热阻阵列可简化成图4的形式,能够大大简化计算过程. 由图4中的简化热阻阵列可得 Rc=(Rs/n1)?(Rf/m1)(Rs/n1)+(Rf/m1)+(Rs/n2)?(Rf/m2)(Rs/n2)+(Rf/m2)+?+(Rs/nn)?(Rf/mn)(Rs/nn)+(Rf/mn).(3)Rc=(Rs/n1)?(Rf/m1)(Rs/n1)+(Rf/m1)+(Rs/n2)?(Rf/m2)(Rs/n2)+(Rf/m2)+?+(Rs/nn)?(Rf/mn)(Rs/nn)+(Rf/mn).(3) 式中:Rc为多孔材料整体等效热阻,K·W-1;mn为单片材料气体个数;nn为单片材料基质个数;