，，

(1. 海軍工程大學(xué) 艦船高溫結(jié)構(gòu)復(fù)合材料研究室，武漢 430033;2. 海軍工程大學(xué) 動(dòng)力工程學(xué)院，武漢 430033)

0 引 言

硬硅鈣石晶體纖維纏繞形成的超輕中空二次粒子是一種新型的多孔絕熱材料，具有超輕的密度、優(yōu)良的絕熱保溫性能、人體友好性和高利用率的特點(diǎn)，近年來(lái)得到高度重視，發(fā)展非常迅速[1]。硬硅鈣石型硅酸鈣絕熱材料按制品密度可分為220 kg/m3左右的普通型和小于135 kg/m3的超輕型兩類，其中孔隙率分別為0.919和大于0.95，屬于高孔隙率的多孔絕熱材料。

多孔絕熱材料的導(dǎo)熱系數(shù)是衡量其絕熱性能的重要參數(shù)。建立起計(jì)算其導(dǎo)熱系數(shù)的模型，并分析各個(gè)因素對(duì)其導(dǎo)熱性能的影響，對(duì)于提高材料的絕熱性能具有十分重要的意義[2]。在研究過(guò)程中，多孔材料內(nèi)部幾何結(jié)構(gòu)的不規(guī)律性導(dǎo)致其物性參數(shù)的不確定性，在理論研究時(shí)，很難用傳統(tǒng)方法對(duì)其導(dǎo)熱性能進(jìn)行準(zhǔn)確的描述。以往的研究大都基于連續(xù)介質(zhì)假設(shè)和容積平均的方法，僅用孔隙率這一個(gè)參數(shù)描述多孔材料的空間結(jié)構(gòu)，但事實(shí)上，孔隙率相同的多孔材料，內(nèi)部固相尺度和分布不同，其導(dǎo)熱系數(shù)也不同。

20世紀(jì)80年代以來(lái)，分形理論開(kāi)始被引入多孔介質(zhì)的傳熱研究中來(lái)[3]。Yongjin Feng、Boming Yu等[4]基于分形理論，利用Sierpinski地毯模擬連續(xù)孔隙的兩相多孔介質(zhì)，推導(dǎo)了其有效熱導(dǎo)率數(shù)學(xué)模型。王唯威等[5]基于分形理論，構(gòu)造了不同分形模型來(lái)模擬分形介質(zhì)，采用有限容積法對(duì)導(dǎo)熱問(wèn)題進(jìn)行了數(shù)值模擬與分析。以上分形模型不含經(jīng)驗(yàn)常數(shù)，參數(shù)較少，計(jì)算簡(jiǎn)便，但不同種多孔介質(zhì)材料內(nèi)部結(jié)構(gòu)的分形特征不同，因此針對(duì)某種特定的材料，推導(dǎo)反映其內(nèi)部微尺度結(jié)構(gòu)特征的傳熱模型具有重要意義。施明恒[6]等利用分形方法對(duì)聚氨酯泡沫塑料微觀泡孔結(jié)構(gòu)和導(dǎo)熱特性進(jìn)行了研究，提出了計(jì)算聚氨酯泡沫塑料導(dǎo)熱的簡(jiǎn)化分形模型。徐平[7]采用分形理論對(duì)泡沫鋁材料的微觀結(jié)構(gòu)進(jìn)行了研究，利用面積-周長(zhǎng)分維算法證明了其微觀結(jié)構(gòu)在一定尺度內(nèi)具有自相似特點(diǎn)。本文的主要工作是針對(duì)所研究的硬硅鈣石超輕中空二次粒子，證明其纖維分布具有分形特征，并建立反映其纖維微觀分布特點(diǎn)的分形導(dǎo)熱模型。

1 硬硅鈣石超輕中空二次粒子纖維壁的微尺度空間結(jié)構(gòu)及分形描述

1.1 微尺度空間結(jié)構(gòu)

硬硅鈣石型硅酸鈣絕熱材料是一種高孔隙率的絕熱材料，是由針狀硬硅鈣石晶體纖維纏繞組成的球狀團(tuán)聚體，外殼密實(shí)，內(nèi)部中空，這種球形聚合體被成為硬硅鈣石二次粒子。

如圖1所示為硬硅鈣石二次粒子表面纖維及孔隙微觀結(jié)構(gòu)的電子顯微鏡掃描圖，可以看出晶體纖維相互交錯(cuò)，形成的幾何結(jié)構(gòu)復(fù)雜多變，且纖維呈隨機(jī)排列，因此，很難用數(shù)學(xué)方法對(duì)其真實(shí)導(dǎo)熱情況加以準(zhǔn)確描述。利用分形的相似性對(duì)其復(fù)雜結(jié)構(gòu)進(jìn)行簡(jiǎn)化，建立簡(jiǎn)化的分形導(dǎo)熱模型，根據(jù)等效熱阻方法，可以推導(dǎo)出其有效導(dǎo)熱系數(shù)的分形計(jì)算式。

圖1 硬硅鈣石二次粒子及纖維壁纏繞的纖維形貌

1.2 分形描述

分形幾何學(xué)是一門以非規(guī)則幾何形狀為研究對(duì)象的幾何學(xué)，如土壤、海岸線、多孔材料等，其內(nèi)部發(fā)生的過(guò)程不能用簡(jiǎn)單的線性近似方法來(lái)認(rèn)識(shí)和描述。分形理論為描述自然界中非線性現(xiàn)象和幾何形狀、物體內(nèi)部的復(fù)雜結(jié)構(gòu)和空間分布提供了一種行之有效的手段，從而為精確研究復(fù)雜結(jié)構(gòu)內(nèi)部客觀規(guī)律開(kāi)辟了一條新路[8]。根據(jù)分形理論，存在于df維歐氏空間中的一個(gè)多孔分形體的面積、體積或長(zhǎng)度等物理量N(r)與物體度量尺度r有以下自相似標(biāo)度關(guān)系：

N(r) ∝rdf。

(1)

式中：N為分形物體的空間占有積(線、面或體)；r為度量尺度；df為分形維數(shù)，可以是整數(shù)也可以是非整數(shù)，兩物體只要滿足分形維數(shù)相等，那么這兩個(gè)物體是自相似的。根據(jù)分形理論，硬硅鈣石超輕中空二次粒子纖維壁的纖維分布若滿足分形特征，須滿足以下關(guān)系式：

N(r)=Crdf。

(2)

式中：C為比例常數(shù)。等號(hào)兩邊取對(duì)數(shù)得：

logN(r)=logC+dflogr。

(3)

查找硬硅鈣石纖維的場(chǎng)發(fā)射掃描電鏡照片可以發(fā)現(xiàn)，二次粒子纖維直徑最大特征尺寸B=100 nm，最小特征尺寸A=35 nm，對(duì)于在B→A的不同長(zhǎng)度度量尺度r，對(duì)硬硅鈣石二次粒子殼體表面纖維電鏡掃描照片進(jìn)行計(jì)盒維數(shù)計(jì)算。采用差分計(jì)盒維數(shù)算法，此算法在計(jì)算灰度圖像時(shí)，簡(jiǎn)單、快速、精確。經(jīng)計(jì)算可以看出硬硅鈣石二次粒子纖維壁的纖維分布滿足分形特征，對(duì)于被測(cè)樣品(表觀密度ρ=125 kg·m-3)其分形維數(shù)為df=2.516 8。

實(shí)際硬硅鈣石二次粒子纖維及孔隙的幾何結(jié)構(gòu)是復(fù)雜多變的，纖維排列既可能是有序的，又可能是無(wú)序排列。利用分形的相似性對(duì)其復(fù)雜的幾何結(jié)構(gòu)進(jìn)行簡(jiǎn)化，假定硬硅鈣石二次粒子纖維壁的纖維排列呈多邊形結(jié)構(gòu)，建立分形等效模型如圖3所示，圖3為硬硅鈣石二次粒子纖維壁截面簡(jiǎn)化模型，進(jìn)而可把每一個(gè)多邊形簡(jiǎn)化為具有相同截面積的正方形網(wǎng)格，即分形單元體。

圖2 硬硅鈣石二次粒子纖維壁電鏡掃描照片分形維數(shù)計(jì)算

圖3 二次粒子纖維壁截面簡(jiǎn)化模型與分形單元體模型

2 分形模型的有效導(dǎo)熱系數(shù)

硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)Ke除與組成材料的各相自身導(dǎo)熱系數(shù)有關(guān)外，還與其空間結(jié)構(gòu)(孔隙率、分形維數(shù))有關(guān)，這里用孔隙率和分形維數(shù)作為描述材料的空間結(jié)構(gòu)的主要參數(shù)。根據(jù)分形理論，有效導(dǎo)熱系數(shù)可以表示為：

Ke=f(∑Ki,Φ,df,T)。

(4)

式中：Ki為硬硅鈣石二次粒子纖維壁的氣相、固相的導(dǎo)熱系數(shù)；φ為硬硅鈣石二次粒子纖維壁纖維的體孔隙率；df為硬硅鈣石二次粒子纖維壁孔隙分形維數(shù)；T為材料的溫度。

根據(jù)孔隙率定義，孔隙率為孔隙體積和總體積的比值，結(jié)合文獻(xiàn)[9]、[10]給出的計(jì)算分形體孔隙率的方法，孔隙率φ為：

(5)

式中：φ為材料的體孔隙率；L為度量尺度。

在纖維長(zhǎng)度方向、熱流方向、熱流垂直方向上同時(shí)取長(zhǎng)度為L(zhǎng)的尺度，構(gòu)成立方體導(dǎo)熱單元體，于是硬硅鈣石二次粒子纖維壁的熱傳遞模型可簡(jiǎn)化為如圖4所示。假設(shè)為一維穩(wěn)態(tài)傳熱，熱流方向自上而下，單元體上下邊界溫度分別為T1和T2，其他邊界取絕熱條件。

圖4 單元體導(dǎo)熱模型和等效熱阻網(wǎng)絡(luò)

從圖4可以看出，左側(cè)通道的氣體熱阻R1為：

R1=L/(KgA1)=2/[Kg(L-D)]。

(6)

中間通道的總熱阻R2為：

R2=R21+R22+R23= (L-D)/(KgA2)+

D/(KsA2)=(L-D)/(KgDL)+1/(KsL)。

(7)

右側(cè)通道的熱阻R3為：

R3=R1=2/[Kg(L-D)]。

(8)

單元體的總熱阻為：

R總=(R1R2R3)/(R2R3+R1R3+

R1R2)= [L(L-D)Ks+LDKg]/{LDKsKg+

Kg(L-D)[Ks(L-D)+KgD]}。

(9)

式中：D為單元體內(nèi)纖維截面的等效正方形的邊長(zhǎng)；A1，A2分別左側(cè)通道和中間通道在熱流方向上的傳熱面積；Kg為氣相導(dǎo)熱系數(shù)；Ks為固相導(dǎo)熱系數(shù)。

由傅立葉定律，總熱流Q為：

Q=(T1-T2)/總=KcA(T1-T2)/L。

(10)

式中：A為單元體總傳熱面積。

分形單元體模型的有效導(dǎo)熱系數(shù)為：

Ke=L/(AR總)=1/(LR總) =(1-D/L)Kg+

1/[(L/D-1)( 1/Kg) +1/Ks)]。

(11)

將(5)式帶入(11)，可得：

Ke=[1-(1-φ)1/df]Kg+1/ {(1/Kg)

[ (1-φ)-1/df-1] +1/Ks)}。

(12)

從(12)式可以看出，硬硅鈣石二次粒子纖維壁的有效導(dǎo)熱系數(shù)與各相自身的導(dǎo)熱系數(shù)有關(guān)，也與材料的微觀結(jié)構(gòu)(分形維數(shù)、孔隙率)有關(guān)。因此在對(duì)材料的熱物性進(jìn)行優(yōu)化設(shè)計(jì)時(shí)，不僅要考慮材料自身的熱物性參數(shù)，還要考慮材料的孔隙率與微觀結(jié)構(gòu)。

3 影響硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)的因素分析

3.1 參數(shù)的確定

固相導(dǎo)熱系數(shù)Ks為硬硅鈣石針狀晶體的導(dǎo)熱系數(shù)，采用來(lái)估計(jì)。這里f1、f2、f3分別是二氧化硅、氧化鈣和水的摩爾百分?jǐn)?shù)。魏高升[11]對(duì)文獻(xiàn)[11]中的數(shù)據(jù)進(jìn)行多項(xiàng)式擬合，得到了二氧化硅、氧化鈣和水的導(dǎo)熱系數(shù)關(guān)于溫度的表達(dá)式：

KSiO2=0.752 33+0.003 13T-4 052 647×

10-6T2+3.526 37×10-9T3。

(13)

KCaO=1.290 23-9.241 18×10-5T-

1.271 45×10-7T2。

(14)

KH2O=-1.088 92+0.011 13T-2.266 28×

10-5T2+1.477 69×10-8T3。

(15)

式中:Kg為環(huán)境壓力下空氣隨溫度變化的導(dǎo)熱系數(shù)，單位為W·m-1·K-1，計(jì)算公式為：

Kg=3.954×10-3+7.720 8×10-5T-

1.608 2×10-8T2。

(16)

將各參數(shù)的表達(dá)式帶入式(12)即可計(jì)算出硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)Ke。

3.2 溫度與硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)的關(guān)系

圖5 硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)隨溫度變化的關(guān)系

取硬硅鈣石二次粒子纖維壁的孔隙率為0.45[12]，取溫度變化范圍為200～1 000 K,得到材料有效導(dǎo)熱系數(shù)隨溫度變化曲線如圖5所示，可以看出隨溫度的升高，有效導(dǎo)熱系數(shù)增大。

3.3 孔隙率與硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)的關(guān)系

孔隙率與硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)的關(guān)系，如圖6所示。

圖6 硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)隨孔隙率變化的關(guān)系

計(jì)算得到材料的分形維數(shù)為df=2.51 68，當(dāng)溫度為T=300 K時(shí)，可得到硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)隨孔隙率變化曲線如圖6所示。由圖6所示，材料的有效導(dǎo)熱系數(shù)隨孔隙率的增大而減小，這是因?yàn)闅庀嗟膶?dǎo)熱系數(shù)遠(yuǎn)小于固相的導(dǎo)熱系數(shù)。當(dāng)材料的孔隙率較小時(shí)，其有效導(dǎo)熱系數(shù)隨之增加而減小的幅度大；當(dāng)材料的孔隙率較大時(shí)，其有效導(dǎo)熱系數(shù)隨之增加而減小的幅度小。

3.4 分形維數(shù)與硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)的關(guān)系

分形維數(shù)與硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)的關(guān)系，如圖7所示。

圖7 硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)隨分形維數(shù)變化的關(guān)系

圖7給出了T=300 K，φ=0.45時(shí)，材料有效導(dǎo)熱系數(shù)隨分形維數(shù)的變化曲線。由圖7可以得出即使是孔隙率相同的同種材料，當(dāng)內(nèi)部固相的尺度或者分布(分形維數(shù)描述)不同時(shí)，其有效導(dǎo)熱系數(shù)也不同。

4 結(jié)束語(yǔ)

(1)采用分形理論分析硬硅鈣石二次粒子纖維壁的有效導(dǎo)熱系數(shù)，建立有效導(dǎo)熱系數(shù)關(guān)于材料孔隙率、分形維數(shù)以及氣固導(dǎo)熱系數(shù)的函數(shù)關(guān)系式。

(2)基于硬硅鈣石二次粒子殼體表面纖維及孔隙微尺度電鏡掃描照片，采用差分盒維數(shù)算法，得到纖維的分布符合分形特征，并計(jì)算了材料的分形維數(shù)，結(jié)合材料的孔隙率，較為精確地描述了材料的微觀空間結(jié)構(gòu)。

(3)硬硅鈣石二次粒子纖維壁的有效導(dǎo)熱系數(shù)不僅與材料固相的自身性能有關(guān)，還與其內(nèi)部的微尺度空間結(jié)構(gòu)(纖維的分布、孔隙率)、材料的使用溫度等因素相關(guān)。

[1] 陳小佳. 超輕硬硅鈣石的制備和性能研究[D]. 武漢: 武漢理工大學(xué)材料學(xué)院, 2007.

[2] 蘇高輝, 楊自春, 王洪宇. 微孔硅酸鈣的有效導(dǎo)熱系數(shù)和最佳使用密度[J]. 應(yīng)用能源技術(shù), 2011,4: 1-5.

[3] 郁伯銘. 多孔介質(zhì)輸運(yùn)性質(zhì)的分形分析研究進(jìn)展[J]. 力學(xué)進(jìn)展, 2003, 33(3): 333-346.

[4] Yongjin Feng, Boming Yu, Mingqing Zou, Duanming Zhang. A generalized model for the effective thermalconductivity of porous media based on self-similarity[J]. Appl. Phys, 2007, 37: 3030-3040.

[5] 王唯威, 淮秀蘭. 分形多孔介質(zhì)導(dǎo)熱數(shù)值模擬分析[J]. 工程熱物理學(xué)報(bào), 2007, 28(5): 835-837.

[6] 施明恒, 李小川,陳永平. 利用分形方法確定聚氨酯泡沫塑料的有效導(dǎo)熱系數(shù)[J]. 中國(guó)科學(xué)E, 2006, 36(5): 560-568.

[7] 徐 平, 楊 昆, 于英華. 泡沫鋁微結(jié)構(gòu)的分形特征及與孔隙率的關(guān)系研究[J]. 材料導(dǎo)報(bào)B, 2011, 25(11): 140-143.

[8] 闞安康, 韓厚德, 湯 偉.基于分形理論的開(kāi)孔聚氨酯泡沫等效導(dǎo)熱系數(shù)研究[J]. 材料導(dǎo)報(bào)B, 2012, 26(2): 143-146.

[9] Dathe A, Thullner M. The relationship between fractal properties of solid matrix and pore space in porous media[J]. Geoderma, 2005, 129(3-4): 279.

[10] Hunt A G. Basic transport properties in natural porous media: Continuum percolation theory and fractal model[J]. Complexity, 2005, 10(3): 22.

[11] Y S Touloukian, P E Livey, S C Saxena. Thermal conductivity: nonmetallic solids. New York: IFI/Plenum press, 1970.

[12] 魏高升. 納米孔超級(jí)絕熱材料的熱物性研究[D]. 北京:北京科技大學(xué)機(jī)械工程學(xué)院, 2006.