康茹雪，胡楊慧，宗若雯,2*，支有冉

(1.中國(guó)科學(xué)技術(shù)大學(xué)火災(zāi)科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，合肥，230026; 2.中國(guó)科學(xué)技術(shù)大學(xué)蘇州研究院蘇州城市公共安全重點(diǎn)實(shí)驗(yàn)室，蘇州，215123; 3.南京工程學(xué)院機(jī)械工程學(xué)院，南京，211167)

不同壓力環(huán)境下不透明非碳化ABS材料的熱解溫度模型研究

康茹雪1，胡楊慧1，宗若雯1,2*，支有冉3

(1.中國(guó)科學(xué)技術(shù)大學(xué)火災(zāi)科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，合肥，230026; 2.中國(guó)科學(xué)技術(shù)大學(xué)蘇州研究院蘇州城市公共安全重點(diǎn)實(shí)驗(yàn)室，蘇州，215123; 3.南京工程學(xué)院機(jī)械工程學(xué)院，南京，211167)

在不同熱流強(qiáng)度下，利用智能步進(jìn)低壓實(shí)驗(yàn)倉(cāng)對(duì)不透明非碳化聚合物材料丙烯腈-丁二烯-苯乙烯(ABS)在不同壓力情況下的熱解過(guò)程中的表面及背面溫度進(jìn)行了實(shí)時(shí)測(cè)量，建立了表面吸收一維熱解模型，預(yù)測(cè)了ABS在不同外界熱流密度下，不同壓力環(huán)境下表面溫度及背面溫度隨時(shí)間的變化；并將實(shí)驗(yàn)測(cè)量值與模型模擬值進(jìn)行對(duì)比。結(jié)果表明，非碳化聚合物的表面吸收一維熱解模型在高熱流密度下對(duì)溫度的預(yù)測(cè)情況優(yōu)于在低熱流密度下對(duì)溫度的預(yù)測(cè)結(jié)果，總體來(lái)說(shuō)，模擬計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)量值較為一致。

一維熱解模型；環(huán)境壓力；表面溫度；背面溫度

0 引言

大部分常用的聚合材料含有大量C、H原子，導(dǎo)致其有很高的易燃性和火災(zāi)危險(xiǎn)性[1]。因此，對(duì)聚合物進(jìn)行燃燒行為研究是非常必要的。熱解是聚合物氣化和其他熱化學(xué)轉(zhuǎn)化過(guò)程的第一步[2]，熱解模型的研究對(duì)預(yù)測(cè)聚合物材料熱解行為有著重要意義。

固體非碳化聚合物熱解會(huì)發(fā)生在環(huán)境壓力低于標(biāo)準(zhǔn)大氣壓的環(huán)境、高緯度地區(qū)以及飛機(jī)和航天器上。高海拔地區(qū)的環(huán)境壓力，如中國(guó)拉薩(海拔高度3650 m)，環(huán)境壓力大約為67 kPa，比標(biāo)準(zhǔn)大氣壓(101 kPa，21%O2)低34%。在飛機(jī)和航天器里，壓力也要低于常壓[3,4]。壓力在聚合物熱解氣化階段和熱表面對(duì)流熱損失方面都有影響，但目前大多數(shù)研究著重于真空或標(biāo)準(zhǔn)大氣壓下的聚合物熱分解過(guò)程[5],針對(duì)聚合物的低壓熱解行為研究少之又少。

為了研究黑色非碳化聚合材料低壓熱解過(guò)程，對(duì)其建立熱解模型是必不可少的。火災(zāi)動(dòng)力學(xué)模擬模型FDS[6]，ThermaKin[7]采用有限的熱解反應(yīng)速率來(lái)模擬材料的熱解過(guò)程，此類模型允許材料在一定溫度范圍內(nèi)發(fā)生熱解，熱解反應(yīng)用一個(gè)或多個(gè)反應(yīng)方程來(lái)表征，材料內(nèi)部熱解產(chǎn)生的揮發(fā)分對(duì)總失重速率有影響。其控制方程為有源項(xiàng)的傳熱方程，質(zhì)量方程為表征化學(xué)反應(yīng)速率的阿倫尼烏斯公式。作為模型輸入量的熱動(dòng)力學(xué)參數(shù)通過(guò)熱重分析法(TG和DTG)在低熱流強(qiáng)度范圍內(nèi)氮?dú)鈿夥罩械玫絒1,8-10]。關(guān)于聚合材料對(duì)輻射熱流的吸收方式，普遍采用兩種：表面吸收和深度吸收。前者認(rèn)為材料對(duì)紅外光譜不透明，所有外部輻射熱均被樣品表面吸收，或者認(rèn)為材料為半透明介質(zhì)，表面不吸收熱量，所有熱量均被材料內(nèi)部吸收。在實(shí)際熱解過(guò)程中兩種吸收方式同時(shí)存在[11]。Jiang等[11]、Delichatsios等[12]對(duì)PMMA涂上黑色碳層，研究其在外加熱輻射作用下的吸收方式，發(fā)現(xiàn)當(dāng)材料的厚度為3 mm時(shí)，其透射率在0.03左右，當(dāng)厚度大于4 mm時(shí)，其透射率接近為0。并且當(dāng)外加輻射小于60 kW/ m2時(shí)，材料的吸收方式接近于表面吸收。Beaulieu和Dembsey[13]也發(fā)現(xiàn)外加熱輻射熱流大于60 kW/m2后，材料的質(zhì)量損失和溫度才會(huì)因?yàn)閮?nèi)部吸收的影響出現(xiàn)明顯的非線性關(guān)系，并且對(duì)于黑色PMMA而言，表面有無(wú)黑色碳層涂料對(duì)其質(zhì)量損失并無(wú)影響。通常情況下，在工程應(yīng)用中為了簡(jiǎn)化計(jì)算，會(huì)將材料的吸收方式簡(jiǎn)化為表面吸收。為了研究低壓下不透明非碳化材料的熱解行為，在不同熱流強(qiáng)度下，在智能步進(jìn)低壓實(shí)驗(yàn)倉(cāng)內(nèi)對(duì)不透明非碳化聚合物材料丙烯腈-丁二烯-苯乙烯(ABS)在不同壓力情況下的熱解過(guò)程中的表面及背面溫度進(jìn)行了實(shí)時(shí)測(cè)量，同時(shí)在實(shí)驗(yàn)數(shù)據(jù)的基礎(chǔ)上，建立一維熱解模型進(jìn)行計(jì)算模擬，并將模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比分析。本文假設(shè)黑色非碳化聚合材料對(duì)紅外光譜為不透明，即外界輻射熱量全部被材料表面吸收，再進(jìn)一步通過(guò)熱傳導(dǎo)加熱表面以下的材料[14-16], 建立了表面吸收一維熱解模型，預(yù)測(cè)了不透明非碳化ABS在不同外界熱流密度下，不同壓力環(huán)境下表面溫度及背面溫度隨時(shí)間的變化；并將實(shí)驗(yàn)測(cè)量值與模型模擬值進(jìn)行對(duì)比。

1 實(shí)驗(yàn)方法及熱解模型

1.1 實(shí)驗(yàn)方法

低壓熱解實(shí)驗(yàn)在中國(guó)科學(xué)技術(shù)大學(xué)火災(zāi)科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室智能步進(jìn)式低壓實(shí)驗(yàn)倉(cāng)里進(jìn)行，熱解實(shí)驗(yàn)平臺(tái)自行搭建，以模擬低壓貧氧環(huán)境，氧氣體積分?jǐn)?shù)控制在4%～5%。低壓熱解試驗(yàn)臺(tái)結(jié)構(gòu)及實(shí)驗(yàn)平臺(tái)結(jié)構(gòu)示意圖如圖1。

黑色聚合物材料丙烯腈-丁二烯-苯乙烯(ABS)材料詳細(xì)信息如表1所示。試驗(yàn)樣品尺寸均為100 mm100 mm3 mm，為保證實(shí)驗(yàn)的一維性，樣品四周及底部用鋁箔紙進(jìn)行包裹，為了防止材料熱解過(guò)程中的熱膨脹，紙上邊沿高出樣件表面10 mm。實(shí)驗(yàn)前分別在樣品上下表面用鋁箔紙固定直線0.5 mmK型熱電偶來(lái)測(cè)量實(shí)驗(yàn)過(guò)程中的表面及背面溫度，熱電偶布置在樣品對(duì)角線交點(diǎn)處。實(shí)驗(yàn)輻射源上方設(shè)置一根熱電偶觀察倉(cāng)內(nèi)溫度變化情況。實(shí)驗(yàn)時(shí)，先用真空泵將密閉實(shí)驗(yàn)倉(cāng)內(nèi)的壓力降至25 kPa，然后沖入氮?dú)庵钡絺}(cāng)內(nèi)壓力升至101.25 kPa，之后重啟真空泵，將倉(cāng)內(nèi)壓力降至設(shè)計(jì)工作壓力。本文中設(shè)計(jì)工作氣壓分別為35 kPa、65.2 kPa和98.1 kPa。 在實(shí)驗(yàn)前通過(guò)熱流計(jì)標(biāo)定所需輻射熱流。在熱流穩(wěn)定之前，輻射源固定在實(shí)驗(yàn)倉(cāng)導(dǎo)軌最左端，熱流穩(wěn)定后，利用滑塊及重錘將輻射源拉至樣品上方，每次實(shí)驗(yàn)保證樣品置于輻射源中央。輻射源尺寸為300 mm300 mm，實(shí)驗(yàn)樣品表面與輻射源底部之間的距離為40 mm，實(shí)驗(yàn)輻射熱流分別為18 kW/m2及36 kW/m2。實(shí)驗(yàn)過(guò)程中低壓實(shí)驗(yàn)倉(cāng)內(nèi)的壓力會(huì)隨著熱解氣的釋放而發(fā)生微小變化，變化范圍在1 kPa以內(nèi)，倉(cāng)內(nèi)非熱解區(qū)溫度變化范圍在4 ℃以內(nèi)，因?yàn)閴毫蜏囟茸兓鄬?duì)熱解區(qū)溫度和壓力很小，因此假設(shè)實(shí)驗(yàn)過(guò)程中倉(cāng)內(nèi)壓力及非熱解區(qū)溫度不變。實(shí)驗(yàn)工況如表2所示。

圖1 低壓熱解實(shí)驗(yàn)平臺(tái)Fig.1 Schematic of experimental setup

表1 樣品詳細(xì)信息

Table 1 Sample information

非碳化聚合物生產(chǎn)商商品名經(jīng)銷商ABS深圳安和達(dá)塑膠制品有限公司丙烯腈-丁二烯-苯乙烯共聚物深圳新諾實(shí)業(yè)塑膠材料有限公司

表2 實(shí)驗(yàn)工況匯總表

1.2 實(shí)驗(yàn)結(jié)果匯總

1.2.1 表面溫度

圖2為各熱流密度下，樣品表面溫度隨環(huán)境壓力變化的變化情況示意圖。由圖2中可以明顯看出，在兩種熱流密度下，表面溫度隨著環(huán)境壓力的上升而下降。在低熱流密度下，下降趨勢(shì)明顯，而在較高熱流密度下，下降趨勢(shì)則并不明顯。這是因?yàn)榈蛪簳?huì)造成熱對(duì)流系數(shù)的減小，從而減小對(duì)流熱損失，增大樣品表面的流入熱量。

圖2 不同環(huán)境壓力下表面溫度曲線圖Fig.2 Profiles of the surface temperature under different gas pressures

1.2.2 背面溫度

圖3為各熱流密度下，樣品背面溫度隨環(huán)境壓力變化的變化情況示意圖。

圖3 不同環(huán)境壓力下背面溫度曲線圖Fig.3 Profiles of the bottom temperature under different gas pressures

很明顯，壓力對(duì)背面溫度的影響較小。背面溫度對(duì)壓力的敏感度沒(méi)有表面溫度對(duì)壓力的敏感度高，這是因?yàn)樵诒趁娼^熱的條件下，材料內(nèi)部及背面對(duì)流熱損失很小，在高熱流下，甚至可以忽略不計(jì)。

1.3 基本熱解模型

在實(shí)驗(yàn)基礎(chǔ)上考慮材料內(nèi)部傳熱及熱解化學(xué)動(dòng)力學(xué)反應(yīng)，假設(shè)材料吸收方式為表面吸收，并通過(guò)控制方程建立一維熱解模型，再耦合邊界和初始條件對(duì)熱解過(guò)程進(jìn)行數(shù)值模擬，熱解模型邊界條件為背面絕熱。為了簡(jiǎn)化模型，假設(shè)材料內(nèi)部熱解氣體產(chǎn)生瞬間便逸出，且材料表面在熱解過(guò)程中接收到的熱流不變。相應(yīng)地物理模型如圖4所示。

圖4 一維熱解物理模型：背面絕熱Fig.4 One-dimension pyrolysis model: Adiabatic bottom condition

1.3.1 控制方程

一維熱解模型控制方程主要包括聚合物內(nèi)能量守恒方程。假設(shè)黑色非碳化聚合材料熱流吸收方式為表面吸收。則聚合物內(nèi)部傳熱控制方程如下所示[17]。

ρsSv[ΔH+(TS-T

(1)

(2)

(3)

(4)

樣品熔融后密度會(huì)發(fā)生變化[18]，

ρs=ρ0(1-βsΔt)

(5)

βs(T)=0.1995×10-6T2-

1.6178×10-5T-1.2271×10-2

(6)

其中，βs為ABS的熱膨脹系數(shù)。

1.3.2 初始條件和表面控制體

初始條件為：

x=δpol:Ts(x,0)=T

(7)

x=0:T1-T2=0

(8)

其中，下標(biāo)pol，分別指聚合物和環(huán)境值。

非碳化聚合物樣品上表面控制體與空氣接觸且有傳熱傳質(zhì)過(guò)程，對(duì)上表面控制體應(yīng)用能量平衡關(guān)系，可得能量平衡方程為：

h(Ts-T

(9)

Qout=

ρsSv[ΔH+(TS-T)(Cs-Cg)]

(10)

樣品下表面邊界絕熱，因此下表面控制體能量平衡方程為：

ρsSv[ΔH+(TS-T

(11)

其中，h可通過(guò)自然對(duì)流經(jīng)驗(yàn)公式算得，即熱平板頂表面的自然對(duì)流公式為[17, 19]，由公式很明顯看出，壓力主要通過(guò)影響氣體的熱擴(kuò)散系數(shù)從而影響熱對(duì)流系數(shù)來(lái)影響流入樣品熱量：

(12)

其中：

(13)

(14)

(15)

(16)

其中，Ra為瑞利數(shù)，L為特征長(zhǎng)度m，g為重力加速度m·s-2，比熱容為常數(shù)，Cp,g=1.063 J·g-1K-1，導(dǎo)熱系數(shù)(kg)、氣體密度(ρg)、氣體熱膨脹系數(shù)(β)和動(dòng)力粘度(νg)是溫度的函數(shù)[20]：

kg(T)=1.5207×10-11T3-4.8574×10-8T2

+1.0184×10-4T-3.9333×10-4

(17)

(18)

(19)

νg(T)=-1.1555×10-14T3+9.5728×10-11T2

+3.7604×10-8T-3.4484×10-6

(20)

2 模擬結(jié)果與討論

根據(jù)上述熱解模型來(lái)計(jì)算ABS不同壓力情況下的表面溫度與背面溫度，并與實(shí)驗(yàn)測(cè)量值進(jìn)行對(duì)比。模型中的網(wǎng)格劃分為0.1mm大小的均勻網(wǎng)格，時(shí)間步長(zhǎng)為0.005s。在計(jì)算表面溫度及背面溫度之前，先對(duì)模型中的每個(gè)偏微分方程進(jìn)行離散化處理，離散方法為顯示差分法。根據(jù)離散化得到的離散化方程進(jìn)行Matlab編程并求解計(jì)算。

2.1 輸入?yún)?shù)

計(jì)算中用到的一些熱物理參數(shù)如表3所示。

表3 模擬輸入的ABS熱物理參數(shù)匯總

2.2 模擬結(jié)果與對(duì)比

2.2.1 表面溫度

表面溫度是ABS熱解過(guò)程中最直觀的參數(shù)之一。在熱解階段的初期，熱解反應(yīng)主要發(fā)生在材料表面附近的一個(gè)薄層，因此熱解速率的強(qiáng)烈程度在于表面溫度的高低。圖5～圖6為不同壓力情況，不同熱流密度情況下，ABS表面溫度實(shí)驗(yàn)值與模擬值的對(duì)比。其中Exp(Experiment)代表實(shí)驗(yàn)測(cè)得的表面溫度，Sim(Simulation)代表模擬計(jì)算得到的表面溫度。

圖5 18 kW/m2輻射強(qiáng)度下ABS表面溫度實(shí)驗(yàn)值與模擬值對(duì)比Fig.5 Comparison of experimental theoretical results of surface temperature of ABS under heat flux of 18 kW/m2

圖6 36 kW/m2輻射強(qiáng)度下ABS表面溫度實(shí)驗(yàn)值與模擬值對(duì)比Fig.6 Comparison of experimental and theoretical results of surface temperature of ABS under heat flux of 36 kW/m2

從圖5、圖6中可以看出，在高熱流密度下，模型計(jì)算得到的表面溫度與實(shí)驗(yàn)值較為接近。但是，在低熱流密度下，模型計(jì)算的表面溫度與實(shí)驗(yàn)值之間的差距較大，但在低壓情況下，模擬值與實(shí)驗(yàn)值的一致性更好。在兩種熱流密度下，實(shí)驗(yàn)值均低于模型計(jì)算值。造成此種情況的原因可能是：1)模型假設(shè)為表面吸收，實(shí)際上還會(huì)有部分熱量被材料內(nèi)部吸收，導(dǎo)致實(shí)驗(yàn)值較低；2)模型內(nèi)輸入的材料熱物理參數(shù)取自文獻(xiàn)，與實(shí)際樣品會(huì)有偏差；3)模型的邊界條件與實(shí)際熱解過(guò)程會(huì)有偏差；4)ABS受熱會(huì)發(fā)生熔融，模型中僅考慮了材料熱膨脹對(duì)熱解造成的影響，有些因素未被完全考慮。低熱流密度下，材料的熱解速率較慢，因此相變的影響比較大，而高熱流密度下，相變的影響較小。5)壓力影響材料的熱對(duì)流系數(shù)，導(dǎo)致低壓下材料表面的流入熱量較多，同時(shí)加快材料的熱解速率，進(jìn)一步減小了材料相變對(duì)模擬情況的影響。

2.2.2 背面溫度

背面溫度是ABS熱解過(guò)程中另外一個(gè)非常直觀的參數(shù)。圖7～圖8是不同壓力情況，不同熱流密度情況下，ABS背面溫度實(shí)驗(yàn)值與模擬值的對(duì)比。

圖7 18 kW/m2輻射強(qiáng)度下ABS背面溫度實(shí)驗(yàn)值與模擬值對(duì)比Fig.7 Comparison of experimental and theoretical results of bottom temperature of ABS under heat flux of 18 kW/m2

低熱流密度下ABS的大規(guī)模熱解時(shí)刻可以從背面溫度曲線反映出來(lái)，也就是材料經(jīng)過(guò)初期溫度迅速升高后到達(dá)拐點(diǎn)的時(shí)刻，該時(shí)刻之后材料背面溫度緩慢升高，材料的凈入能量通過(guò)熱解損失了大部分。從圖7和圖8中可以很明顯的看出，模擬值中拐點(diǎn)的出現(xiàn)較明顯于實(shí)驗(yàn)值，這是因?yàn)樵趯?shí)際熱解過(guò)程中，除過(guò)熱傳導(dǎo)吸收的熱量，樣品內(nèi)部還會(huì)通過(guò)深度吸收得到一部分熱量，這導(dǎo)致實(shí)際流入樣品內(nèi)部的熱量要稍高于模型中流入的熱量。熱解初期模擬值的背面溫度要稍高于實(shí)驗(yàn)值，這是因?yàn)槟Ｐ图僭O(shè)為背面完全絕熱狀態(tài)，但是在實(shí)際的低熱流熱解情況中，背面的石棉板仍然會(huì)吸收一部分熱量，導(dǎo)致材料背面的溫度低于模擬值。高熱流密度作用下的背面溫度模擬結(jié)果同樣好于低熱流密度，這是因?yàn)榈蜔崃髅芏认?，材料的熱解不受熱滯后效?yīng)的影響，所以熱解的逐層熱解效應(yīng)稍弱。高熱流密度下模擬值與實(shí)驗(yàn)值拐點(diǎn)的出現(xiàn)相比低熱流密度下較為不明顯，這是因?yàn)樵诟呙芏葻崃饕约盁釡笮?yīng)的雙重作用下，樣品底部通過(guò)對(duì)流，傳導(dǎo)損失的熱量微乎其微，因此材料的熱解更接近與模型假設(shè)的逐層熱解。此外，從圖7、圖8中可以看到，在兩種熱流密度下，常壓下拐點(diǎn)的出現(xiàn)相較于低壓較為明顯，這也說(shuō)明，壓力導(dǎo)致的熱對(duì)流損失對(duì)常壓情況下的熱解影響更大。

圖8 36 kW/m2強(qiáng)度下ABS背面溫度實(shí)驗(yàn)值與模擬值對(duì)比Fig.8 Comparison of experimental and theoretical results of bottom temperature of ABS under heat flux of 36 kW/m2

3 結(jié)論

開(kāi)展了低壓下不透明非碳化聚合物ABS的熱解實(shí)驗(yàn)研究，對(duì)其熱解過(guò)程中的表面溫度和背面溫度進(jìn)行了實(shí)時(shí)測(cè)量，基于能量守恒，建立了表面吸收一維熱解模型，預(yù)測(cè)了不透明非碳化ABS在不同外界熱流密度下，不同壓力環(huán)境下表面溫度及背面溫度隨時(shí)間的變化，并將實(shí)驗(yàn)測(cè)量值與模型模擬值進(jìn)行對(duì)比。主要結(jié)論如下：

1)實(shí)驗(yàn)過(guò)程中除了熱傳導(dǎo)，材料內(nèi)部還會(huì)直接吸收部分輻射熱，因此模擬結(jié)果會(huì)普遍高于實(shí)驗(yàn)值。此外，ABS的熱物理參數(shù)和實(shí)驗(yàn)初始條件均與實(shí)驗(yàn)樣品有誤差，所以會(huì)影響模擬結(jié)果。高熱流密度和低壓情況下，因?yàn)閷?duì)流熱損失小，樣品熱解速率較快，所以壓力和相變對(duì)樣品的模擬情況影響較小，模擬值與實(shí)驗(yàn)值也較為接近。

2)總體來(lái)說(shuō)，兩種熱流密度下，表面溫度的實(shí)驗(yàn)值和模擬值較為一致。高熱流密度下的模擬情況要好于低熱流密度下，低壓下的模擬情況要優(yōu)于高壓。

3)樣品背面溫度的實(shí)驗(yàn)值與模擬值也較為一致。在兩種熱流密度下，常壓下拐點(diǎn)的出現(xiàn)相較于低壓較為明顯，表明壓力導(dǎo)致的熱對(duì)流損失對(duì)常壓情況下的熱解影響更大。

[1] Li J, Stoliarov SI. Measurement of kinetics and thermodynamics of the thermal degradation for non-charring polymers[J]. Combustion and Flame, 2013, 160(7): 1287-1297.

[2] Basile L, et al. Influence of pressure on the heat of biomass pyrolysis[J]. Fuel, 2014, 137: 277-284.

[3] Lange KE, et al.Bounding the spacecraft atmosphere design space for future exploration missions[R]. National Aeronautics and Space Administration, 2005.

[4] McAllister S, et al. Piloted ignition delay of PMMA in space exploration atmospheres[J]. Proceedings of the Combustion Institution. 2009, 32(2): 2453-2459.

[5] Koptelov AA, Karyazov SV. Effect of pressure on the thermal decomposition of polymers[J]. Doklady Physical Chemistry, 2003, 389(4): 101-105.

[6] Lautenberger C, et al. The application of a genetic algorithm to estimate material properties for fire modeling from bench-scale fire test data[J]. Fire Safety Journal, 2006, 41(3): 204-214.

[7] Stoliarov SI, Lyon RE. Thermo-kinetic model of burning for pyrolyzing materials[J]. Fire Safety Science, 2008, 9: 1141-1152.

[8] 蘇誠(chéng)偉, 等. 聚甲醛的熱重分析研究[J]. 功能高分子學(xué)報(bào), 1994, 7(2): 165-171.

[9] 曾文茹, 等. 非等溫動(dòng)力學(xué)積分-微分法研究聚苯乙烯熱解反應(yīng)機(jī)理及動(dòng)力學(xué)參數(shù)[J]. 高分子材料科學(xué)與工程, 2006, 22(5): 162-165.

[10] 曾文茹, 等. 非線性等轉(zhuǎn)化率法研究聚苯乙烯熱解反應(yīng)活化能與轉(zhuǎn)化率的關(guān)系[J]. 高分子材料科學(xué)與工程, 2008, 24(8): 128-131.

[11] Jiang FH, et al. Absorption of thermal energy in PMMA by in-depth radiation[J]. Fire Safety Journal, 2009, 44(1): 106-112.

[12] Delichatsios MA, et al. Upward fire spread: key flammability properties, similarity solutions and flammability indices[J]. Fire Safety Science, 1991, 3: 217-226.

[13] Beaulieu PA, Dembsey NA. Flammability characteristics at applied heat flux levels up to 200 kW/m2[J]. Fire and Materials, 2008, 32(2): 61-86.

[14] Whiting P, et al. A one-dimensional mathematical model of laser induced thermal ablation of biological tissue[J]. Lasers in Medical Science, 1992, 7(1): 357-368.

[15] Staggs JEJ. Modelling thermal degradation of polymers using single-step first-order kinetics[J]. Fire Safety Journal, 1999, 32(1): 17-34.

[16] Delichatsios MA, Chen Y. Asymptotic, approximate, and numerical solutions for the heatup and pyrolysis of materials including reradiation losses[J]. Combustion and Flame, 1993, 92(3): 292-307.

[17] Incropera FP, et al. Fundamentals of heats and mass transfer[M]. John Wiley & Sons. Incorporation, 2007, 6.

[18] 黃永華, 等. 熱膨脹系數(shù)簡(jiǎn)易測(cè)量裝置研制及若干材料測(cè)量[J]. 化工學(xué)報(bào), 2016, 67(S2): 38-45.

[19] Atreya A. Convection heat transfer, in SFPE handbook of fire protection[M]. Society of Fire Protection Engineers, 2016, 5.

[20] Gong JH, et al. A numerical study of thermal degradation of polymers: Surface and in-depth absorption[J]. Applied Thermal Engineering, 2010, 106:1366-1379.

[21] Li J, Stoliarov SI. Measurement of kinetics and thermodynamics of the thermal degradation for charring polymers[J]. Polymer Degradation and Stability, 2014, 106: 2-15.

[22] Li J, et al. Gasication experiments for pyrolysis model parameterization and validation[J]. International Journal of Heat and Mass Transfer, 2014,77: 738-744.

Pyrolysis temperature model of opaque non-charring ABS under different pressures

KANG Ruxue1， HU Yanghui1， ZONG Ruowen1, 2，ZHI Youran3

(1.State Key Laboratory of Fire Science, University of Science and Technology of China, Hefei 230026, China; 2. Suzhou Key Laboratory of Urban Public Safety, Suzhou Institute for Advanced Study, University of Science and Technology of China, Suzhou 215123, China; 3.Institute of Mechanical Engineering，Nanjing Institute of Technology Nanjing，Nanjing 211167, China)

A series of experiments on the opaque Acrylonitrile Butadiene Styrene (ABS) in reduced low-pressure chamber have been conducted in this work to study the pyrolysis of this polymer under different external heat fluxes. The variation of temperature at the surface and bottom of the sample during the whole process has been measured. A surface absorption one-dimension numerical model is developed to predict the surface and bottom temperature of ABS in a series of reduced pressures under different external heat fluxes. The model is verified by the measured experimental data. The results indicates that the theoretical result of temperature under higher heat flux is better than that under lower heat flux. The experimental results fit well with the theoretical results.

One-dimension pyrolysis model; Ambient pressure; Surface temperature; Bottom temperature

2016-12-22；修改日期：2017-02-27

國(guó)家自然科學(xué)基金面上項(xiàng)目(51576185，51606092)；國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFC0802101)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(WK2320000034)

康茹雪(1992-)，女，陜西咸陽(yáng)人。中國(guó)科學(xué)技術(shù)大學(xué)火災(zāi)科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室碩士研究生，安全科學(xué)與工程專業(yè)，主要從事聚合物熱物理及煙氣毒性等方面的研究。

宗若雯，E-mail：zongrw@ustc.edu.cn

1004-5309(2017)-00079-08

10.3969/j.issn.1004-5309.2017.02.03

X915.5

A