基于有限元分析的真空玻璃傳熱性能數(shù)值模擬研究

李 宏，李璟瑋，陳 鵬，俞昊天，熊德華

(武漢理工大學(xué)硅酸鹽建筑材料國家重點(diǎn)實(shí)驗(yàn)室，武漢 430070)

0 引 言

隨著工業(yè)化與城市化的持續(xù)推進(jìn)，全球能源與環(huán)境問題日益嚴(yán)重，開發(fā)具有低能耗、低消耗、輕污染、可循環(huán)利用等特點(diǎn)的綠色建筑材料逐漸成為發(fā)展趨勢，也是我國實(shí)施可持續(xù)發(fā)展的必然要求[1]。我國每年新建建筑中85%以上仍是高能耗建筑，其運(yùn)行能耗碳排放占全國總排放21.5%[2-5]。在圍護(hù)結(jié)構(gòu)中，門窗是耗能的主要部位，門窗損失的能量約占建筑圍護(hù)部件總能耗的40%～50%。因而，高性能節(jié)能門窗對減少建筑的外圍圍護(hù)結(jié)構(gòu)能耗起著至關(guān)重要的作用。節(jié)能玻璃種類繁多，大致可分為中空玻璃、低輻射(low-E)玻璃和真空玻璃三類。中空玻璃與低輻射玻璃是目前常用的節(jié)能玻璃，通常是制成低輻射中空玻璃使用，這些玻璃雖然具有一定的節(jié)能效果，但是還不能夠完全滿足低能耗建筑的要求。

真空玻璃是一種具有優(yōu)異保溫、隔熱、降噪性能的玻璃深加工產(chǎn)品[6-9]。真空玻璃由兩片平板鋼化玻璃構(gòu)成，其中一片或兩片涂有低輻射膜；其截面與結(jié)構(gòu)示意如圖1所示，真空玻璃邊緣多采用低熔點(diǎn)玻璃粉或金屬焊料密封[10-11]；中間用微小的支撐物均勻隔開，形成0.1～0.3 mm的真空腔[12]。真空玻璃優(yōu)異的保溫隔熱性能源自其獨(dú)特的傳熱機(jī)制。如圖2所示，與中空玻璃相比，真空玻璃僅由支撐物傳熱、輻射傳熱以及殘余氣體傳熱三部分組成[13]。因此，真空玻璃具有更好的保溫、隔熱以及隔聲、降噪和防結(jié)露性能，是少有的傳熱系數(shù)低至1.0 W·m-2·K-1，可完全滿足低能耗建筑要求的節(jié)能玻璃。

圖1 真空玻璃結(jié)構(gòu)示意圖

圖2 真空玻璃和中空玻璃的傳熱過程示意圖

真空玻璃的制備與服役過程中涉及大量微小變量的計(jì)算，這使得有限元模擬在真空玻璃中有著廣泛而深入的應(yīng)用。孫敢[14]借助有限元分析軟件ANSYS，建立溫度場模型，模擬分析真空平板玻璃激光焊接的成型過程，并結(jié)合微觀組織結(jié)構(gòu)分析檢驗(yàn)激光焊接的可行性，來解決真空平板玻璃邊緣封接速度過慢的問題。李彥兵等[15]針對鋼化真空玻璃中支撐點(diǎn)的排布對玻璃力學(xué)特性的重要影響，通過建立鋼化真空玻璃力學(xué)模型，數(shù)值分析了不同支撐點(diǎn)排布方式下鋼化玻璃的力學(xué)性能。研究結(jié)果表明,隨著支撐間距增大，單位面積支撐點(diǎn)數(shù)急劇減少，使鋼化玻璃最大變形量和最大應(yīng)力迅速增大。Ghosh等[16-17]構(gòu)建了一種由普通真空玻璃和玻璃基薄膜太陽能電池膠接而成的半透明光伏真空玻璃組件，通過模擬仿真軟件，計(jì)算得出整體的傳熱系數(shù)為0.8 W·m-2·K-1，相較于雙層光伏玻璃，在保證發(fā)電效率的同時(shí)具有更佳的保溫隔熱性能。Zhu等[18]引入有限元分析法，研究邊緣封接材料厚度以及支撐柱布放方式對真空玻璃整體傳熱系數(shù)的影響，并通過現(xiàn)場實(shí)驗(yàn),輔助驗(yàn)證模擬計(jì)算的準(zhǔn)確性。結(jié)果表明,降低封接層寬度，選用導(dǎo)熱系數(shù)更小的陶瓷支撐材料以及增大支撐柱的布放間距可有效降低傳熱系數(shù)。

本文以前期制備的小尺寸真空玻璃樣品[19-21]為基礎(chǔ)，利用有限元分析法對大尺寸真空玻璃的傳熱性能進(jìn)行模擬計(jì)算，并進(jìn)一步構(gòu)建建筑模型，分析真空玻璃對節(jié)能建筑保溫隔熱性能的提升效果。

1 有限元模擬

1.1 真空玻璃物理模型的建立

1.1.1 真空玻璃三維幾何模型建立

本研究建立三種規(guī)格的真空玻璃幾何模型，如圖3所示：第一種為100 mm×100 mm，第二種為200 mm×200 mm，最后一種為500 mm×500 mm。以下是各組成部位的尺寸。基板玻璃：500 mm×500 mm×4 mm、200 mm×200 mm×4 mm或100 mm×100 mm×4 mm；邊緣封接層：寬度10 mm，厚度0.3 mm；支撐柱：直徑1 mm，高度0.3 mm，間距為20 mm。

圖3 三種規(guī)格真空玻璃的幾何模型

1.1.2 材料屬性設(shè)置

選擇COMSOL Multiphysics 5.6軟件自帶材料庫中的White plate glass[solid]為基板玻璃材料，選擇材料Air為真空腔參與氣體材料，自定義低熔點(diǎn)玻璃材料為支撐柱和邊緣封接材料各項(xiàng)屬性，具體材料屬性如表1所示(表1中低熔點(diǎn)玻璃性能參數(shù)源于課題組前期研究，其余數(shù)據(jù)源自材料庫)。

表1 組成真空玻璃的主要材料性能參數(shù)

1.1.3 物理場設(shè)置

由于存在氣體和固體，所以選擇固體和流體傳熱物理場進(jìn)行耦合，定義邊界條件和初始條件模擬真空玻璃傳熱情況。

初始條件：所有組件初始溫度為293.15 K(20 ℃)。真空腔內(nèi)殘余為流體，氣壓為0.1 Pa，導(dǎo)熱系數(shù)為0，其余組件為固體。

邊界條件[22]：內(nèi)側(cè)環(huán)境溫度為293.15 K(20 ℃)，換熱系數(shù)為3.6 W·m-2·K-1；外側(cè)環(huán)境溫度為253.15 K(-20 ℃)，換熱系數(shù)為16 W·m-2·K-1。

1.1.4 網(wǎng)格劃分及計(jì)算公式

有限元網(wǎng)絡(luò)的劃分技術(shù)是將每個(gè)域劃分成各種形狀的更小尺寸的子域，這些子域稱為有限元。基元的角點(diǎn)則稱為節(jié)點(diǎn)，在角點(diǎn)間也有節(jié)點(diǎn)，場變量就在節(jié)點(diǎn)上進(jìn)行計(jì)算。基元的節(jié)點(diǎn)越多，計(jì)算誤差越小，微小部位也可得到體現(xiàn)，同樣計(jì)算耗時(shí)也越長[23]。采用物理場控制網(wǎng)絡(luò)中的較細(xì)化網(wǎng)絡(luò)對所建立的幾何模型進(jìn)行有限元網(wǎng)絡(luò)劃分，劃分結(jié)果如圖4所示。

圖4 有限元網(wǎng)絡(luò)劃分圖

COMSOL Multiphysics軟件是基于有限元方法來求解偏微分方程的大型仿真軟件。有限元分析是通過離散法將解域離散為一系列相互關(guān)聯(lián)的有限元，最后形成離散的線性代數(shù)方程組，然后通過求解器在單元節(jié)點(diǎn)處求解方程組。固體和流體傳熱物體場模型的瞬態(tài)傳熱方程如公式(1)所示。

(1)

式中：T為溫度，K；t為時(shí)間，s;Cp為恒壓熱容，J·kg-1·K-1；ρ為密度，kg·m-3；Q為熱源，W·m-3；k為導(dǎo)熱系數(shù)，W·m-1·K-1；u為流體流速，m·s-1。

1.2 節(jié)能建筑物理模型的建立

1.2.1 節(jié)能建筑二維幾何模型建立

本研究對節(jié)能建筑使用平板玻璃、中空玻璃、真空玻璃的三種情況進(jìn)行模擬。節(jié)能建筑的二維幾何模型如圖5(a)所示，屋內(nèi)空間面積為10 m×5 m，鋁合金門為0.1 m×1.1 m。保溫墻墻體結(jié)構(gòu)為15 mm水泥砂漿+200 mm石墨聚苯板+25 mm水泥砂漿+200 mm鋼筋混凝土+5 mm混合砂漿，如圖5(b)所示。

圖5 建筑模型

1.2.2 材料屬性設(shè)置

真空玻璃各材料屬性如表1所示，節(jié)能建筑所使用的其他各類材料具體屬性如表2所示。

表2 其他材料性能參數(shù)

1.2.3 物理場設(shè)置

模擬節(jié)能建筑在冬天的傳熱情況，由于存在氣體和固體，所以采用的是固體和流體傳熱物理場結(jié)構(gòu)。真空玻璃為(4 mm基板玻璃+0.3 mm真空腔+4 mm基板玻璃)×8 m，真空腔內(nèi)空氣為流體，氣壓為0.1 Pa，導(dǎo)熱系數(shù)為0，其余組件為固體。中空玻璃為(4 mm基板玻璃+5 mm空腔+4 mm基板玻璃)×8 m，空腔內(nèi)為干燥的空氣。該空腔為不透風(fēng)的腔室，因此其傳導(dǎo)、輻射和對流引起的熱流可以用等效導(dǎo)熱系數(shù)keq來表示，如公式(2)所示。

(2)

式中：d是矩形腔在熱流速度方向的尺寸；R是矩形腔的熱阻，如公式(3)所示。

(3)

式中：ha是對流傳熱系數(shù)，hf是輻射傳熱系數(shù)，這兩個(gè)系數(shù)如公式(4)和公式(5)所示。

(4)

(5)

式中：C1=0.025 W/(m·K)；σ=5.67·10-8W/(m2·K4),是斯蒂芬-玻爾茲曼常數(shù)；Tm是矩形腔邊界上的平均溫度；E是界面輻射率，如公式(6)所示。

(6)

式中：ε1和ε2是表面發(fā)射率(在本模型中均為0.84)。

F是矩形截面的視角系數(shù)，如公式(7)所示。

(7)

式中:d是矩形腔在熱流速度方向的尺寸；b是矩形腔在垂直于熱流速度方向的尺寸。

經(jīng)過計(jì)算，中空玻璃空腔等效導(dǎo)熱系數(shù)為0.047 W·m-1·K-1。

定義模型的初始條件和邊界條件。

初始條件：所有組件初始溫度為293.13 K。

邊界條件：內(nèi)側(cè)環(huán)境溫度為293.15 K，外側(cè)環(huán)境溫度為273.15 K，對流熱通量中換熱系數(shù)為16 W·m-2·K-1。

1.2.4 網(wǎng)格劃分及計(jì)算公式

采用物理場控制網(wǎng)絡(luò)中的較細(xì)化網(wǎng)絡(luò)對所建立的幾何模型進(jìn)行有限元網(wǎng)絡(luò)劃分，劃分結(jié)果如圖6所示。

圖6 有限元網(wǎng)絡(luò)劃分圖

固體和流體傳熱物理場二維模型的瞬態(tài)傳熱方程如公式(8)所示。

(8)

由于數(shù)學(xué)物理方程描述物理場在無限空間的問題，所以在實(shí)際模型中需要對其邊界條件進(jìn)行限定。通過模擬計(jì)算在上述初始和邊界條件下7 d內(nèi)使用普通平板玻璃、中空玻璃和真空玻璃的節(jié)能建筑溫度的變化情況，來評估其保溫隔熱性能。

2 結(jié)果與討論

2.1 真空玻璃保溫隔熱性能數(shù)值模擬

當(dāng)真空玻璃規(guī)格為100 mm×100 mm時(shí)，為模擬真空玻璃的保溫性能，對模型表面施加以下邊界條件：內(nèi)側(cè)環(huán)境溫度為293.15 K，換熱系數(shù)為3.6 W·m-2·K-1；外側(cè)環(huán)境溫度為253.15 K，換熱系數(shù)為16 W·m-2·K-1。得到了0～60 min內(nèi)側(cè)玻璃和外側(cè)玻璃表面的溫度分布圖，如圖7所示。外側(cè)玻璃表面溫度從外框向內(nèi)部中心溫度逐漸降低，從265.37 K降至257.20 K，邊緣密封層區(qū)域溫度為265.37 K，中間支撐柱溫度比腔內(nèi)空氣溫度高。內(nèi)側(cè)玻璃表面溫度從外框向中心溫度逐漸上升，從266.68 K升至278.88 K，邊緣密封層區(qū)域溫度為266.68 K，中間支撐柱溫度比腔內(nèi)空氣溫度低。繪制內(nèi)側(cè)玻璃表面和外側(cè)玻璃表面平均溫度隨時(shí)間變化情況。在一定時(shí)間內(nèi)，隨著時(shí)間的增加，四個(gè)表面平均溫度均隨之降低，當(dāng)時(shí)間為30 min時(shí)，四個(gè)表面的平均溫度趨向平穩(wěn)，內(nèi)側(cè)玻璃表面平均溫度、外側(cè)玻璃表面平均溫度分別維持在271.55 K和261.98 K。

圖7 真空玻璃表面溫度及溫度隨時(shí)間變化圖(100 mm×100 mm)

當(dāng)真空玻璃的規(guī)格為200 mm×200 mm，施加相同的初始條件和邊界條件，得到0～60 min內(nèi)側(cè)玻璃和外側(cè)玻璃表面的溫度分布圖，如圖8所示。外側(cè)玻璃表面溫度從外框向內(nèi)部中心逐漸降低，從266.03 K降至255.86 K；內(nèi)側(cè)玻璃表面溫度從外框向內(nèi)部中心溫度逐漸升高，從266.82 K升至285.66 K。繪制內(nèi)側(cè)玻璃表面和外側(cè)玻璃表面平均溫度隨時(shí)間變化情況。在一定時(shí)間內(nèi)，隨著時(shí)間的增加，表面平均溫度均隨之降低，當(dāng)時(shí)間為30 min時(shí)，四個(gè)表面的平均溫度趨向平穩(wěn)，內(nèi)側(cè)玻璃表面平均溫度、外側(cè)玻璃表面平均溫度分別維持在277.43 K、259.66 K。相比于100 mm×100 mm的真空玻璃，200 mm×200 mm內(nèi)側(cè)溫度更高，外側(cè)溫度更低，這說明200 mm×200 mm的保溫性能更好。

圖8 真空玻璃表面溫度及溫度隨時(shí)間變化圖(200 mm×200 mm)

當(dāng)真空玻璃的規(guī)格為500 mm×500 mm，施加相同的初始條件和邊界條件，得到0～60 min內(nèi)側(cè)玻璃和外側(cè)玻璃表面的溫度分布圖，如圖9所示。外側(cè)玻璃表面溫度從外框向內(nèi)部中心逐漸降低，從266.08 K降至255.89 K；內(nèi)側(cè)玻璃表面溫度從外框向內(nèi)部中心溫度逐漸升高，從266.79 K升至 286.43 K。繪制內(nèi)側(cè)玻璃表面和外側(cè)玻璃表面平均溫度隨時(shí)間變化情況，如圖9所示。在一定時(shí)間內(nèi)，隨著時(shí)間的增加，表面平均溫度均隨之降低，當(dāng)時(shí)間為30 min時(shí)，表面的平均溫度趨向平穩(wěn)，內(nèi)側(cè)玻璃表面平均溫度、外側(cè)玻璃表面平均溫度分別維持在282.26 K、257.73 K。相比100 mm×100 mm及200 mm×200 mm的真空玻璃，500 mm×500 mm內(nèi)側(cè)平均溫度更高，外側(cè)平均溫度更低，這說明500 mm×500 mm的保溫性能更好。

圖9 真空玻璃表面溫度及溫度隨時(shí)間變化圖(500 mm×500 mm)

2.2 真空玻璃在節(jié)能建筑中傳熱性能數(shù)值模擬

為了模擬真空玻璃在冬天環(huán)境下的保溫性能，對模型表面施加以下邊界條件：內(nèi)側(cè)環(huán)境溫度為293.15 K，外側(cè)環(huán)境溫度為273.15 K，換熱系數(shù)為16 W·m-2·K-1時(shí)，得到了7 d時(shí)使用普通平板玻璃、中空玻璃和真空玻璃的節(jié)能建筑的溫度分布圖和等溫線分布圖。從圖10中可以發(fā)現(xiàn)：建筑中各個(gè)部位的保溫能力各不相同，保溫墻體的保溫性能最好，門窗是主要的散熱部位。當(dāng)建筑物使用平板玻璃時(shí)，從圖10(a)和(b)中可以發(fā)現(xiàn)7 d時(shí)窗戶處的溫度為277.24 K，然后通過熱輻射影響室內(nèi)溫度，以波的形式從277.24 K逐漸升至289.91 K，此時(shí)室內(nèi)平均溫度為285.04 K，由于窗戶的面積比門大很多，所以此時(shí)熱量流失的主要途徑是窗戶。當(dāng)使用中空玻璃時(shí)，從圖10(c)和(d)中可以發(fā)現(xiàn)7 d時(shí)窗戶處溫度為282.52 K，門的溫度為273.97 K，此時(shí)室內(nèi)平均溫度為286.70 K，熱量大部分通過門流失，但窗戶也損失一部分熱量。當(dāng)使用真空玻璃時(shí)，從圖10(e)和(f)中可以發(fā)現(xiàn)絕大部分熱量通過門流失，此時(shí)窗戶的保溫性能十分優(yōu)良，室內(nèi)平均溫度為288.95 K，比使用平板玻璃高出3.91 K，比中空玻璃高出2.25 K，說明真空玻璃替換平板玻璃和中空玻璃可以達(dá)到更好的節(jié)能效果，滿足冬暖夏涼的要求。

圖10 7 d時(shí)建筑物的溫度分布圖(a)、(c)、(e)和等溫線分布圖(b)、(d)、(f)及建筑室內(nèi)溫度隨時(shí)間變化圖(g)

3 結(jié) 論

利用COMSOL Multiphysics 5.6的熱傳導(dǎo)模塊進(jìn)行真空玻璃傳熱性能的數(shù)值模擬，通過建立真空玻璃的幾何模型、設(shè)置材料屬性、定義初始和邊界條件、劃分有限元網(wǎng)絡(luò)、模擬計(jì)算和后處理分析結(jié)果來模擬真空玻璃和真空玻璃在節(jié)能建筑中的傳熱性能。

(1)隨著真空玻璃規(guī)格的增大，真空腔體積變大，邊緣密封層占整個(gè)真空玻璃表面面積百分比減小，邊緣密封層的傳熱對真空玻璃的傳熱影響減少，真空玻璃的保溫性能得到增強(qiáng)。500 mm×500 mm規(guī)格的真空玻璃內(nèi)側(cè)玻璃和外側(cè)玻璃表面平均溫度分別維持在282.26 K、257.73 K,保溫性能強(qiáng)于200 mm×200 mm及100 mm×100 mm規(guī)格的真空玻璃。

(2)從節(jié)能建筑的二維模型數(shù)值模擬可以看出，使用真空玻璃7 d時(shí)，室內(nèi)平均溫度比使用平板玻璃高出3.91 K,比中空玻璃高出2.25 K，說明真空玻璃可以更好地滿足居住的舒適性以及節(jié)能環(huán)保的要求。