石墨的热传导

1、For personal use only in study and research; not for commercial use莂石墨的热传导 (heat conduction of graphite)膈石墨体内存在温度梯度时，热量从高温处向低温处的流动。表征石墨导热能力的参数是热导率。热导率入是单位时间内、单位面积上通过的热量q( 热流密度)与温度梯度 grad T 之间的比例系数。蒄 q=入 grad(1)膅式中负号表示热流方向与温度梯度方向相反。式 (1) 常称为热传导的傅里叶定律。假如垂直于x轴方向的截面积为 S,材料沿x轴方向温度梯度为dT/dx，在时间内，沿x轴正方向流过 S

2、截面的热量为 Q在稳定传热状态下，式(1)具有如下的形式:芈热导率的法定单位是WmSKS对于不稳定传热过程，即物体内各处温 度随时间而变化。与外界无热交换，本身存在温度梯度的物体，随着时间的推 移，温度梯度会趋于零，即热端温度不断降低和冷端温度不断升高，最终达到 一致的平衡温度。在这种不稳定传热过程中，物体内单位面积上温度随时问的3T变化率二为：薂式中T为时间，P为密度，Cp为质量定压热容。入/ p Cp常称为石墨的热扩散率或导温系数，常用单位为cmf/s袀热传导是通过导热载体的运动来实现的。石墨的导热载体有电子、声子（晶格振动波）、光子等。石墨的热导率可表示为各种导热载体的贡献的迭加:莈A=

3、S(l/3)ciZi1芅式中Vi、li、Ci分别为导热载体i的运动速度、平均自由程和单位体积的 比热容。石墨的各种导热载体之间又相互作用、相互制约。例如不同频率的声 子之间互相碰撞、产生散射，声子与晶界、点阵缺陷和杂质之间也产生散射， 影响其平均自由程。因此，石墨的热传导是一个极为复杂的物理过程。理论上 准确预测各种石墨的热导率数值及其随温度的变化，虽然有过长期的艰苦工 作，但仅取得了有限的成绩。粗略地说，在常温和不太高的温度下（小于2000K），声子热导率占压倒优势，电子及光子的热导可以忽略不计。在极低温 度下（小于10K）电子热导才占有一定的分量。光子热导要在很高的温度下 （2000K以上

4、）才开始出现。石墨的热导率随其电导率的增大而升高 （见威德曼?弗朗兹定律）莄石墨单晶 纯净的天然鳞片石墨、高定向热解石墨，这些石墨晶体，缺陷较少而且尺寸较大，一般可认为是较完善的石墨单晶。对这类石墨的热导有过 相当多的研究。在压应力下，经过 3000K 以上处理的热解石墨，其体积密度为 2.25g/cm 3，接近单晶的理论密度 2.266g/cm 3，其(002) 衍射峰半宽角展只有 0.4 °( 镶嵌角 ) ，也十分接近于理论值零度。这种石墨的热导率见表1。这些数值一般认为可代表单晶石墨的相应数值。沿两个主方向的热导率：沿层面的记 为入a,沿垂直于层面的则记为入c0羂在常温下入a比

5、入c大200倍左右。温度升高，这个比值有所下降，但仍然很 大。所以由微晶组成的多晶石墨，其热导为微晶层面热导率入a所控制，入c几乎可不予考虑。天然鳞片石墨的 入a在常温下为280500W/(m?K)之间，比值入丿入c 在35之间，可见其晶体的完善程度远不如咼定向热解石墨。莈晶体结构高度规整的热解石墨，La在2000nm以上，由低温到高温，其导热率随温度的变化呈钟罩形，见图 1、图 20蚆在温度远低于石墨晶体层面热导的特征温度 匚下:螂入 a* exp( 9 ”(5)/bT)蚁式中b约等于2, 9 '有时称做德拜温度，但与表征热容的德拜温度不同(见炭质材料和石墨材料的热容)。在温度远高于

6、9阳寸，则有蒈入 a * T - 1(6)肇按式(5)，在低温下，入a随温度T的增高而上升；按式(6)，在高温下，入则随温度的增高而下降。在低温和高温之间，(5) 、(6) 两式都起作用，在这两种作用互相匹敌时，入a达到最大值。这就是形成钟罩形曲线的原因。蒄在不太低的温度下，石墨晶体的导热载体是声子，式 (3) 可简化为：蒀入=y pcVvl薈式中P为密度，Cv为质量定容热容，V为声子传播速度，I为声子两次散 射或碰撞之间的平均自由程，Y为比例系数。在低温下，I的大小由晶界散射所 制约，I的大小与微晶的尺寸相当。所以入aT曲线峰值的高度和位置为石墨晶 体的尺寸（微晶a向直径La）所控制。热解石

7、墨的退火温度越高，晶体越完善， La随之增大，因而热导率入a增高，峰值增大，峰位向低温侧移动（图3）。膄两种石墨晶体，晶粒a向直径分别为La.1和La.2，热导率峰位分别为Tm.1和 Tm.2，这些参数之间有如下关系：(8)艿提供了一种由热导率数据估算La的方法。由这种方法得到的La数值与由X光衍射法的大体相当蚇热导椭球晶体两个主方向的热导率为入a和入c,沿任一方向的热导率为入,为这一方向与晶轴c的交角，有薅入© =入asin22+入cCOS(9)蚄式（9）pT形象地用以长径为旋转轴的一个旋转椭球来表示 （图4）。椭球的 半长径为入c-1/2，半短径为入a-1/2。这一椭球称为石墨的

8、热导椭球。在任一方向的 热导率入®,可由椭球在该方向上的半径丫。来表示：节入 ©=1/ 丫 r（10）螇在该方向上的半径越短，热导率越大。羆层面热导率理论 石墨晶体热导率的理论，十分繁杂，依靠计算机的帮助 取得了不少进展，但还有不少问题有待进一步的探讨。兹仅以无缺陷理想石墨 晶体的层面热导率入a为例，把晶格振动波加以量子化，形象地把振动波称为声 子，振动波是向量，可称为波矢。波矢的能量和状态是晶体倒易点阵的函数。整个晶体的倒易点阵可用一个小区域来代表；这一区域叫做布里渊区。只要把 声子在这一区域内的能量和状态搞清楚，声子在整个晶体内的情况也就了如指 掌了。膂石墨晶体的布里渊

9、区是一个六角棱柱体(图5)。如果只讨论石墨晶体层面 的热导率，作为一种简化模型，只讨论声子在图5的正六角形面上的运动情况就够了。这种二维情况使问题大为简化，处理较为方便。用n代表波数，在nx，ny平面上，六角形截面的面积，可用一个半径为 nm的圆面来代表，由图5 得出：肁=2kZ/3 =0.246X10®(11)袇式(11)中a是石墨一个晶格参数，a=0.246 x 10-8cm nm就是声子振动的最大波数，即声子在单位长度上的振动次数。声子运动速度 v与波数n的乘积 是声子的频率，声子的能量与频率成正比。声子的最大角频率 wm=2nvnm，而2n nm称为最大角波数，常记为 qm。

10、qm=1.55X108cmt莇把声子的运动情况加以分类，每一类称为一个声子分支，每一分支给予一 个代号。在布里渊区的正六角形层面上有好几个声子分支，主要的有3个：1.LA，纵向分支，最大频率为 37THz速度为vl=2.36 X 106cm/s； 2.TA，横向分 支，最大频率为25THz速度为VT=1.59X106cm/s; 3低TA分支，又称为弯曲振 动分支，最大频率为14THz速度为Vb=0.53 X 106cm/s。此外还有折叠LA分 支、横向光学分支TO等，这些非主要分支的频率都低于 4THz而且与其他分 支发生强烈的相互作用，因此小于 4THz即角频率小于wt=2.5 X1013S

11、“的这些 分支，在热量传输中不起什么作用，可以忽略不计。wc称为声子角频率下限。低TA分支的速度与LA、TA相比低很多，也可不予考虑。在这种大为简化的情 况下，只考虑LA、TA这两个分支，并且只考虑热导，不涉及热容。这就是所谓 二维声子气模型。由此可定义一个德拜速度Vd：螀由以上列举的数据得到：德拜速度 vd=1.86 xi06cm/s,声子最大角频率wm=VDDcm=2.88x1014s-1。袇在热导载体为声子所垄断，即在常温和不太高的温度下，理想石墨晶体的层面热导率为入a.id，贝U(13)芁式中p为理想石墨晶体的密度2.266g/cm3, y为格林爱森系数（见石墨的热 容），可取Y =2

12、,由此得到袃入 a.id=5.73/T X105羇此式简捷明了，又显然为式（6）的T1关系提供了理论依据。由此式算得的热导率与高度完善的高定向热解石墨实测数值的对比见表2。羄实测值与理论值大体相适应，由十分简化的理论模型得到的结果竟然与实 际符合得如此之好。两者之比平均为 0.94，这表明即使如此高度完美的石墨晶 体，其完善程度与理想晶体相比仍有不足之处。肃多晶石墨 多晶石墨的热导率为众多因素所左右：骨料与黏结剂的种类和 配比、成型条件、热处理温度等制造工艺有显著的影响；微晶的尺寸与分布、 孔隙的数量和形状等结构因素，其影响尤为突出。不同石墨品种之间，热导率 千差万别，即使同一种石墨，不同批次

13、之间也有相当大的差异。影响因素虽 多，但控制热导率的基本规律不变。在以声子热导为主的温度区界内，仍为式 （7） 所表明的规律所控制。薁多晶石墨由众多的微晶组成。多晶石墨的热导通过微晶的层面传递（a 向热导），因为微晶的入a比入c约大两个数量级，C向热导可忽略而不计，如图6所 示。在中等温度下，微晶的入a主要为两种散射过程所控制：1.晶界散射所控制的热导入B，微晶尺寸La越大，入B越大。2.声子间互相碰撞引起的散射所控制的 热导入u,温度越高，这种散射越强烈，入u随温度的增高而减小。入a、入B、入u之 间有如下关系：肇 1/ 入 a=1/ 入 B+1/ 入 u(15)莅在任一方向(X方向)的热导

14、率入X取决于多晶石墨中微晶的取向和分布。由于热量传递的路径蜿蜒曲折，微晶之间还可能存在非晶态及不完善的晶态炭素 物质，过渡性炭素物质，入X与入a之间的关系中应列入一个校正系数a X，即：(1莀由理论分析，入u随温度的变化数据列在表3中。再把不同温度下热导率的实测数据与理论式(16)比较，即可得到入B和ax。对一种挤压成型的核石墨PGA和模压成型的ZTA石墨，其热导率的实测值与计算值的对比表示在图7腿表3入u随温度的变化螆温度/ K膃100腿150芆200腿250袅300膂350莆400芄500莃600羁700蒆800蚅900肅1000螀入u/W-螀(cm K)-1肆39120453.626.7

15、20.114.912.19.298.006.876.205.615.15热导率随温度而变化的情况，对几种模压石墨，分别表示在图8图9上，入-T曲线都呈钟罩形。高热导石墨 挤压成型的宇航石墨ATJ- S,密度为1.84g/cm3,以及各向同性的细颗粒高密度石墨，密度达 2.0g/cm3HDG HDFG用短纤维增强的HDG都是高热导多晶石墨。这些石墨的热导率随温度而变化的情况见图10。热导率与密度早在19世纪中叶，著名物理学家、电磁波理论的创始人J.C.麦克斯韦（Maxwell）。在其名著电磁波理论（1873）中就指出：对含有孔 隙的材料，设孔隙是以等径小球的形状均匀分散在材料中，材料的传导率 （

16、电导或热导），从理论上可由下式计算:l p1 + 0.5P(17)式中P为孔隙率，入o为无孔(P=0)时的热导率。此式具有历史意义。对于石墨，孔隙并非呈球状，更非等径，此式当然不适用。但它表明孔隙率越大(即密度越小)，热导率越小。这一定性结论却正确无误。一种挤压成型的、经过不 同浸渍处理的核石墨，在常温下，其热导率入/随孔隙率的变化符合如下关系：入 / =入 oexp(-bP)(18)式中入o=128OW/(n?K)，为无孔隙时的极限热导率，常数 b=7.00同一类型的石墨，热导率随其密度的增大而上升，图11表示HDFGC性石墨的入与密度的关系。热处理温度 多晶石墨大多是由焙烧毛坯经高温热处理

17、制成，热处理温度越高，微晶的发育越完善，La增大，热导率也随之增大。用煅后石油针状焦及中温煤沥青，经挤压成型做成的焙烧小棒，经不同热处理(HTT)后，其La的数值见表4。其轴向热导率入/随温度变化的情况见图12。热导率的倒数1/入称为 热阻。在不同热处理温度下，这种石墨的轴向热阻 1/入/与其l/L a的关系见图13。也是用石油焦和中温煤沥青做成的另一种挤压石墨，图14显示出其入/依赖于La的情况。对于一种模压石墨，其 入丄与HTT之间的关系见图15。热扩散系数a又称为导温系数，a =入/ P Cp 0 （见式（3）。它表征材料在加热或冷却过程中，各部分温度趋向于一致的能力；是在不稳定传热过程

18、中，说 明温度变化速度的一个特性参数。材料的导温系数越高，材料内部温度的传播 速度越大，材料内的温差就越小。一种高密度，P = 1 .8 1 g/Cm 3、各向同性细颗粒石墨EK- 98,其a随温度的变化情况见图16 上0热散逸系数£表征石墨材料热性能的一个综合参数，与热导率密切相关， 其定义为：£ =（入 Cpp ）1/2（19）在法定单位制中，£的单位是WS2 ?m 2?K：它表征材料表面散热或吸热 能力的大小。EK- 98石墨的热散逸系数随温度变化情况示于图17。热导异向度 石墨材料的各向异性在热导上表现为沿平行对称轴方向的热导 率入/与沿垂直方向的热导率 入丄的差异上。一般，对挤压石墨 入/ >入丄，把 入/ /入丄这一比值称为热导异向度；对模压石墨，入丄>入/,则把比值入丄/入/称为热导异向度；即异向度最小为 1（同向性）。设沿石墨对称轴oz的取向 参数为R