，，

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

0 引 言

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

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

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

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

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

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

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

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

1.2 分形描述

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

N(r) ∝rdf。

(1)

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

N(r)=Crdf。

(2)

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

logN(r)=logC+dflogr。

(3)

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

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

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

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

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é)合文獻[9]、[10]給出的計算分形體孔隙率的方法，孔隙率φ為：

(5)

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

在纖維長度方向、熱流方向、熱流垂直方向上同時取長度為L的尺度，構(gòu)成立方體導(dǎo)熱單元體，于是硬硅鈣石二次粒子纖維壁的熱傳遞模型可簡化為如圖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)纖維截面的等效正方形的邊長；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)。因此在對材料的熱物性進行優(yōu)化設(shè)計時，不僅要考慮材料自身的熱物性參數(shù)，還要考慮材料的孔隙率與微觀結(jié)構(gòu)。

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

3.1 參數(shù)的確定

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

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，計算公式為：

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

1.608 2×10-8T2。

(16)

將各參數(shù)的表達式帶入式(12)即可計算出硬硅鈣石二次粒子纖維壁有效導(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)系

計算得到材料的分形維數(shù)為df=2.51 68，當溫度為T=300 K時，可得到硬硅鈣石二次粒子纖維壁有效導(dǎo)熱系數(shù)隨孔隙率變化曲線如圖6所示。由圖6所示，材料的有效導(dǎo)熱系數(shù)隨孔隙率的增大而減小，這是因為氣相的導(dǎo)熱系數(shù)遠小于固相的導(dǎo)熱系數(shù)。當材料的孔隙率較小時，其有效導(dǎo)熱系數(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時，材料有效導(dǎo)熱系數(shù)隨分形維數(shù)的變化曲線。由圖7可以得出即使是孔隙率相同的同種材料，當內(nèi)部固相的尺度或者分布(分形維數(shù)描述)不同時，其有效導(dǎo)熱系數(shù)也不同。

4 結(jié)束語

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

(2)基于硬硅鈣石二次粒子殼體表面纖維及孔隙微尺度電鏡掃描照片，采用差分盒維數(shù)算法，得到纖維的分布符合分形特征，并計算了材料的分形維數(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ì)輸運性質(zhì)的分形分析研究進展[J]. 力學(xué)進展, 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é)報, 2007, 28(5): 835-837.

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

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

[8] 闞安康, 韓厚德, 湯 偉.基于分形理論的開孔聚氨酯泡沫等效導(dǎo)熱系數(shù)研究[J]. 材料導(dǎ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] 魏高升. 納米孔超級絕熱材料的熱物性研究[D]. 北京:北京科技大學(xué)機械工程學(xué)院, 2006.