方奕忠 沈韓 崔新圖 黃臻成 馮饒慧 廖德駒 龐曉寧

摘 要 準穩(wěn)態(tài)法是測量不良熱導(dǎo)體熱導(dǎo)率的一種常用方法， 樣品材質(zhì)、保溫層、加熱功率等都會對測量結(jié)果造成影響。本文采用分離變量法， 以有機玻璃為例， 求解了厚度均勻的無限大不良熱導(dǎo)體兩側(cè)同時施加均勻的指向中心面的熱流密度時的熱傳導(dǎo)方程，得到樣品各處的溫度隨時間變化的嚴格解。在初始加熱時間不到十分鐘的狀態(tài)下采用截斷近似， 并在準穩(wěn)態(tài)下用常功率平面熱源加熱法測量了不良熱導(dǎo)體的熱導(dǎo)率k及比熱c。在此基礎(chǔ)上， 采用四塊有機玻璃板和兩個加熱薄膜外包泡沫和聚丙烯（PP）阻熱材料的模型， 作了COMSOL多物理場耦合仿真， 得到的理論、實驗及仿真結(jié)果三者符合得很好。

關(guān)鍵詞 準穩(wěn)態(tài); 不良熱導(dǎo)體; 熱導(dǎo)率; 虛實結(jié)合; COMSOL仿真模擬

準穩(wěn)態(tài)法是測量不良熱導(dǎo)體熱導(dǎo)率的常用方法， 測量結(jié)果受樣品材質(zhì)、保溫層性質(zhì)、熱源位置、加熱器形狀及功率等多種因素影響。本文根據(jù)虛實結(jié)合實驗[1-7]設(shè)計的理念， 采用理論、實驗和仿真相結(jié)合的方式探討不良熱導(dǎo)體內(nèi)部的溫度分布及變化過程。理論求解了一個厚度均勻無限大不良熱導(dǎo)體平板兩側(cè)同時施加均勻的指向中心面的熱流密度時板內(nèi)各處的溫度變化情況。在準穩(wěn)態(tài)情況下， 用常功率平面熱源加熱法[8]測量了四塊有機玻璃樣品板和兩塊加熱膜構(gòu)成的體系中樣品的熱導(dǎo)率及比熱[9]， 并采用COMSOL 多物理場耦合仿真軟件[10]分三種模型仿真了實驗過程，包括雙樣品模型、四樣品理想模型、四樣品實驗?zāi)Ｐ偷?，仿真模型?fù)雜程度逐步提升，逼近真實實驗。結(jié)果證明， 理論、實驗、COMSOL 仿真模擬三方面符合得很好。

1 理論模型及求解

考慮一無限大不良熱導(dǎo)體平板， 厚度為2L，初始溫度為T0，在平板兩側(cè)同時施加均勻的指向中心面x =0的熱流密度qc，如圖1所示。平板上各點x 在t 時刻的溫度值為T =T （x，t），則T滿足定解問題

2 實驗測量及結(jié)果

實驗裝置如圖2所示。實驗樣品做成長方體， 截面為正方形。為了盡可能滿足無限大平板的理論假設(shè)， 樣品的厚度小于方形截面的邊長的1/6。實驗采用的不良熱導(dǎo)體樣品為有機玻璃，樣品厚度L=0.010m，密度ρ=1196kg/m3，共有四塊，與樣品架配合使用。

為了精確測量加熱面的熱流密度qc， 實驗采用薄膜形的加熱器， 如圖2所示。加熱膜厚度l=0.123mm， 邊長與樣品相同， 由厚度為0.03mm、寬約1.50mm 的銅箔細長條片迂回近似均勻地嵌在聚酰亞胺薄膜的中間。給銅箔通電即可產(chǎn)生熱量， 其加熱功率在整個加熱面上均勻并可精確控制。四塊樣品兩兩夾持薄膜形加熱器， 保證加熱膜兩側(cè)的熱阻相同， 可認為熱流密度為功率密度的一半。測溫用的兩個T型熱電偶的熱端分別安裝在中心面和其中一個加熱面（取樣品3和樣品4之間的加熱面）的中心處。則通過樣品的熱流密度可表示為

qc =U2/（2Fr） （單位為W/m2） （27）

其中，U 為兩并聯(lián)加熱膜的加熱電壓， F=S/A 為邊緣修正后的加熱面積。修正是因為有機玻璃樣品不是無限大平板， 且樣品的四周用聚丙烯（PP）板和中發(fā)泡泡沫板阻熱， 做不到完全絕熱。經(jīng)實驗和COMSOL模擬可知， 該修正因子可取A =0.85。S=0.09m×0.09m 為修正前的實際加熱面積， 每塊加熱膜的電阻r ≈109.74Ω， 若U =19V， 則qc=172.60W/m2。

實驗中分別測量中心面溫度和加熱面溫度隨時間的變化關(guān)系。兩個熱電偶對應(yīng)的冷端均放在同一個保溫杯里面， 保溫杯保持為室溫， 故兩個冷端的溫度保持相同。我們從室溫T0=295.70K開始給加熱膜通電， 通電電壓保持為U =19V，記錄加熱面熱電偶與中心面熱電偶之間的溫差熱電勢U1 和中心面熱電偶相對于室溫冷端的中心面熱電勢U2。大約加熱10分鐘后， 溫差熱電勢U1 將保持不變， 而中心面熱電勢U2 隨時間近似線性上升， dU2/dt 近似為常數(shù)。實驗采用的銅康銅熱電偶（T 型熱電偶）在室溫附近的電壓溫度系數(shù)λ=0.0404mV/K，據(jù)此可分別算出中心面和加熱面之間的溫度差ΔT1=U1/λ 和中心面的溫度T2=U2/λ+T0。ΔT1 和T2 隨加熱時間的變化關(guān)系如圖3所示。加熱10分鐘后， ΔT1=5.050K， dT2/dt=0.0101K/s。代入式（25）和式（26）， 得熱導(dǎo)率k=0.1709W/（m·K）， 比熱c=1429J/（kg·K）。

3 熱傳導(dǎo)過程的COMSOL仿真

3.1 雙樣品理想傳熱模型的仿真

考慮兩塊緊貼在一起的有機玻璃板，每塊厚度L=0.010m，邊長d=0.090m，則截面面積S=d2。樣品兩端以熱流密度qc 均勻地注入熱量。在兩塊板中間分界面的中心放置一域點探針1，右邊注入熱流密度的表面中心處放置另一域點探針2，所有外表面絕熱。仿真時取qc =U2A/（2rS ）=172.60W/m2，有機玻璃密度ρ=1196kg/m3，熱導(dǎo)率k=0.1709W/（m·K），恒壓比熱容c=1429J/（kg·K），網(wǎng)格化時兩塊板都選細化或較細化，如圖4所示。COMSOL 模擬結(jié)果分別見圖5、圖6和圖7?？梢姶蠹s8分鐘后，加熱面和中心面之間的溫度差基本不變，ΔT1=5.04K，與實驗值相比相對誤差在±0.10%內(nèi);而中心面溫度隨時間的變化率dT2/dt=0.0100K/s，與實驗值相比相對誤差在±0.50%內(nèi)。

3.2 四樣品理想傳熱模型的仿真

考慮四塊同材料同尺寸的勻質(zhì)有機玻璃板緊靠在一起，每塊板厚度L=0.010m，邊長d=0.090m。如圖2所示，在第1塊、第2塊樣品之間和第3塊、第4塊樣品之間分別放置一加熱薄膜，故需考慮6塊長方體的熱傳導(dǎo)。加熱薄膜由銅箔及聚酰亞胺薄膜組成，銅箔密度ρ1 =8960kg/m3，恒壓比熱容c1 =385J/（kg·K），熱導(dǎo)率k1=400W/（m·K）;聚酰亞胺密度ρ2=1420kg/m3，恒壓比熱容c2=1090J/（kg·K），熱導(dǎo)率k2=0.5W/（m·K），銅箔厚度l2=0.03mm，寬1.5mm，按迂回的長條片均勻嵌在聚酰亞胺薄膜的中間薄層。由于整個加熱膜厚度l=0.123mm，銅箔的厚度占了約1/4，銅箔的面積約等于整個加熱膜面積的2/3，故銅箔體積約占整個加熱膜的1/6，聚酰亞胺約占5/6?？扇〖訜崮ふw密度ρ=16ρ1+56ρ2=2676.8kg/m3，恒壓比熱容c=16c1+56c2=972.5J/（kg·K），熱導(dǎo)率k=16k1+56k2=67.09W/（m·K）。四塊有機玻璃板和兩塊加熱膜構(gòu)成一個整體，且整體結(jié)構(gòu)的外表面均設(shè)為絕熱。兩個加熱膜的加熱熱流均取qc =U2A/rS，其中的修正因子A =0.85。網(wǎng)格化時加熱膜的單元尺寸需選較細化或極細化，有機玻璃板可選細化，如圖8 所示。COMSOL模擬結(jié)果分別示于圖9、圖10和圖11，可見大概加熱8分鐘后，加熱面與中心面之間的溫度差基本恒定，ΔT1=4.99K，與實驗值相比相對誤差落在±0.60%內(nèi);中心面域點探針1的溫度對時間的變化率dT2/dt=0.0100K/s，與實驗值相比相對誤差落在±0.50%內(nèi)。表明四樣品理想模型仿真結(jié)果與實驗結(jié)果符合得很好。

3.3 四樣品真實傳熱模型的仿真

在四樣品理想模型基礎(chǔ)上，在左右兩端面用中發(fā)泡厚泡沫板阻熱，每塊泡沫板邊長122mm，厚47mm，恒壓比熱容c3=4200J/（kg·K），熱導(dǎo)率k3 =0.046W/（m·K），密度ρ3 =400kg/m3。四塊樣品及兩塊加熱膜共六塊板的四周用PP板進行阻熱，所有阻熱板和泡沫板外表面均與周圍環(huán)境的空氣有熱交換，滿足第三類邊條件，且取熱交換系數(shù)h=150W/（m2·K）。PP 板的底部高為45mm，左右邊寬均為63mm，頂部高為40mm，模型中四周PP板的總厚度均為40.246mm。PP板密度ρ4 =920kg/m3，恒壓比熱容c4 =1900J/（kg·K），熱導(dǎo)率k4 =0.90W/（m·K），室溫取T0=295.70K，加熱膜加熱功率Q=U2/r，不含修正因子，U 和r 取值與四樣品理想模型一致，網(wǎng)格化時加熱膜的單元尺寸可選較細化（或超細化），有機玻璃板可選常規(guī)（或細化），見圖12。用COMSOL模擬的結(jié)果分別見圖13、圖14和圖15?？梢?，當加熱時間達8分鐘后，加熱面和中心面的溫度差穩(wěn)定在ΔT1=5.15K，與實驗值相比相對誤差落在±0.99%內(nèi);中心面的升溫速率dT2/dt=0.0105K/s，與實驗值相比相對誤差落在±2.0%內(nèi)。仿真結(jié)果與實驗及前兩種模型的仿真結(jié)果基本符合。圖16為中心截面的等溫線分布，對稱性好。

4 結(jié)語

本文在無限大不良熱導(dǎo)體導(dǎo)熱時溫度變化的嚴格解的基礎(chǔ)上，采用COMSOL 多物理場耦合仿真軟件分別仿真了雙樣品理想模型、四樣品理想模型和四樣品真實模型傳熱過程中的溫度分布及溫度變化情況，仿真結(jié)果與理論、實驗基本符合。清楚說明了在實驗設(shè)計中由于作了簡化，需要在熱流密度中引入修正因子A ，且取A =0.85是合理的。

參 考 文 獻

[1] 廖德駒， 沈韓， 崔新圖， 等. 基于通用設(shè)備和Multisim 的虛實結(jié)合電學(xué)實驗教學(xué)模式[J]. 物理與工程， 2021，31（1）：96-103.

LIAO D J， SHEN H， CUI X T， et al. Virtual-reality combinedteaching mode for electric experiments of college physicslaboratory base on common used instruments and Multisim[J]. Physics and Engineering， 2021，31（1）： 96-103. （inChinese）

[2] 陳曉東， 高梓軒， 何辛濤， 等. 基于近場電磁掃描系統(tǒng)的矩形金屬波導(dǎo)模式成像實驗[J]. 物理與工程， 2021，31（4）：160-165.

CHEN X D， GAO Z X， HE X T， et al. Imaging experimentof rectangular metal waveguide modes based on nearfieldelectromagnetic waves scanning systems[J]. Physicsand Engineering， 2021，31（4）：160-165. （in Chinese）

[3] 趙福利， 陳敏， 董建文， 等. 在過程性學(xué)習(xí)評價基礎(chǔ)上開展線上+線下《光學(xué)》課程教學(xué)[J]. 物理與工程， 2021， 31（4）：154-160.

ZHAO F L， CHEN M， DONG J W， et al. Online and offlineteaching implementation of optics with process evaluation[J]. Physics and Engineering， 2021， 31（4）： 154-160.（in Chinese）

[4] 馮饒慧， 沈韓， 高鳳彬， 等. 以衍射實驗為例談虛實結(jié)合實驗教學(xué)模式[J]. 物理與工程， 2020， 30（4）： 102-107.

FENG R H， SHEN H， GAO F B， et al. Virtual-actual experimentalteaching mode for diffraction experiment[J].Physics and Engineering， 2020， 30（4）： 102-107. （in Chinese）

[5] 方奕忠， 沈韓， 崔新圖， 等. 基于CompactDAQ 的多通道溫度測量實驗及COMSOL模擬[J]. 物理與工程， 2021， 31（1）：109-113.

FANG Y Z， SHEN H， CUI X T， et al. Experiment ofmulti-channel temperature measurement based on Compact-DAQ and COMSOL simulation[J]. Physics and Engineering，2021， 31（1）： 109-113. （in Chinese）

[6] 崔新圖，方奕忠，沈韓，等. 有限長通電圓柱導(dǎo)體溫度分布的嚴格解及COMSOL 模擬[J]. 大學(xué)物理，2018，37（8）：13-18.

CUI X T， FANG Y Z， SHEN H， et al. The exact solutionand COMSOL simulation of temperature inside a finitelength energized cylinder[J]. College Physics， 2018， 37（8）：13-18. （in Chinese）

[7] 陳天賜， 黃臻成， 丁喜冬， 等. 常功率平面熱源法仿真及其在熱學(xué)教學(xué)的應(yīng)用[J]. 物理與工程， 2021， 31（Z1）： 61-65.

CHEN T C， HUANG Z C， DING X D， et al. Simulation ofplane heat source method with constant heat rate and its applicationin thermodynamic course[J]. Physics and Engineering，2021， 31[Z1]： 61-65. （in Chinese）

[8] 王補宣， 韓禮鐘， 王維城， 等. 同時測定熱絕緣材料α 和λ的常功率平面熱源法[J]. 工程熱物理學(xué)報， 1980， 1（1）：80-87.

WANG B X， HAN L Z， WANG W C， et al. A plane heatsource method for simultaneous measurement of the thermaldiffusivity α and conductivity λ of insulating materials withconstant heat rate[J]. Journal of Engineering Thermophysics，1980， 1（1）： 80-87. （in Chinese）

[9] 沈韓. 基礎(chǔ)物理實驗[M]. 北京： 科學(xué)出版社， 2015： 258-264.

[10] 王剛，安琳. COMSOL Multiphysics 工程實踐與理論仿真：多物理場數(shù)值分析技術(shù)[M]. 北京： 電子工業(yè)出版社， 2012：72-83.