殷　謙，杜文靜，紀(jì)興林，程　林

?

一種回轉(zhuǎn)窯余熱回收用集熱器的實(shí)驗(yàn)研究及其結(jié)構(gòu)優(yōu)化

殷謙，杜文靜，紀(jì)興林，程林

（山東大學(xué)熱科學(xué)與工程研究中心，山東 濟(jì)南 250061）

摘要：減少回轉(zhuǎn)窯表面的散熱損失對(duì)冶金、化工等高耗能行業(yè)的節(jié)能減排具有重要意義。針對(duì)現(xiàn)有窯表面集熱器回收效率較低的現(xiàn)狀，提出了一種新型集熱器，并通過實(shí)驗(yàn)測(cè)量與數(shù)值計(jì)算分析了其換熱性能。建立了集熱器換熱量與管長(zhǎng)、管徑及管子數(shù)等結(jié)構(gòu)參數(shù)的數(shù)學(xué)關(guān)系，分別以換熱面積、進(jìn)出口壓降以及改進(jìn)的熵產(chǎn)數(shù)為優(yōu)化目標(biāo)，利用遺傳算法分別對(duì)其進(jìn)行了結(jié)構(gòu)設(shè)計(jì)，結(jié)果表明，優(yōu)化后集熱器的輻射換熱所需的換熱面積減小15%，對(duì)流換熱所需的換熱面積減少20%左右。優(yōu)化后兩種形式的換熱面所消耗的泵功明顯減少。優(yōu)化后集熱器的傳熱熵產(chǎn)數(shù)與優(yōu)化前無明顯變化，而流動(dòng)熵產(chǎn)數(shù)得到明顯降低。

關(guān)鍵詞：節(jié)能；回轉(zhuǎn)窯；集熱器；傳熱；優(yōu)化

引 言

水泥工業(yè)是國民經(jīng)濟(jì)發(fā)展的基礎(chǔ)性工業(yè)，同時(shí)也是高耗能行業(yè)。一些水泥生產(chǎn)線的能量分析表明生產(chǎn)1 t水泥需要消耗至少3～4 GJ的能量[1-2]。為降低其能耗，一些研究者[3- 4]提出利用余熱煙氣進(jìn)行發(fā)電或供熱等節(jié)能措施。另外，一些研究表明回轉(zhuǎn)窯表面散熱損失占生產(chǎn)設(shè)備表面散熱的 28.6%～56.6%，占系統(tǒng)總能耗的5%～15%[5-6]。因此，減少回轉(zhuǎn)窯表面散熱對(duì)水泥工業(yè)的節(jié)能具有重要意義。

一些學(xué)者對(duì)回轉(zhuǎn)窯表面集熱器裝置進(jìn)行了相關(guān)研究。Karamarkovi?等[7]設(shè)計(jì)了一種預(yù)熱空氣的集熱器，其他研究者[8-9]提出集熱器更多地應(yīng)用于加熱冷卻水或者鍋爐給水。為了提高集熱器的回收熱量的效率，研究人員設(shè)計(jì)了不同的集熱器布管方式，包括沿窯體周向布置[6]以及軸向布置[8]等。另外，集熱器結(jié)構(gòu)參數(shù)設(shè)計(jì)也有相關(guān)理論研究。Caputo 等[6]對(duì)集熱器回收效率和結(jié)構(gòu)參數(shù)之間的關(guān)系進(jìn)行了數(shù)學(xué)建模，通過調(diào)整管長(zhǎng)參數(shù)降低了回轉(zhuǎn)窯的能量消耗率。馬萬龍[10]建立了回轉(zhuǎn)窯表面與集熱器之間的傳熱模型，對(duì)管長(zhǎng)等結(jié)構(gòu)參量進(jìn)行了初步優(yōu)化。上述集熱器的設(shè)計(jì)均為回轉(zhuǎn)窯的低溫輻射換熱實(shí)現(xiàn)熱量的回收，不僅傳熱效率較低，而且容易引起窯體的超溫，影響設(shè)備的正常運(yùn)行?；谏鲜龇治觯菊n題組[11]提出了一種新型集熱器結(jié)構(gòu)形式，該結(jié)構(gòu)不但改進(jìn)了現(xiàn)有集熱器的管束形式，使其便于加工制造，而且增加了對(duì)流換熱面以提高換熱性能。

本文以上述實(shí)際應(yīng)用于水泥回轉(zhuǎn)窯的集熱器為研究對(duì)象，理論分析和實(shí)驗(yàn)測(cè)定了集熱器的換熱量以及換熱特性?；趽Q熱器優(yōu)化[12-13]的理論，建立了集熱器換熱量與結(jié)構(gòu)參數(shù)之間的數(shù)學(xué)關(guān)系，分別以換熱面積、壓降和改進(jìn)的熵產(chǎn)數(shù)為目標(biāo)函數(shù)，優(yōu)化設(shè)計(jì)集熱器輻射和對(duì)流換熱面的管長(zhǎng)、管徑以及管子數(shù)等結(jié)構(gòu)參數(shù)，分析討論了不同運(yùn)行工況下優(yōu)化設(shè)計(jì)的結(jié)果。

1　集熱器的結(jié)構(gòu)及實(shí)驗(yàn)研究

1.1 集熱器整體結(jié)構(gòu)

集熱器結(jié)構(gòu)如圖1所示。集熱器安裝于直徑為4 m的回轉(zhuǎn)窯上，每個(gè)集熱器沿窯軸向長(zhǎng)度為2 m。6段冷卻水管排焊接起來構(gòu)成了輻射換熱的水冷壁管排，蛇形管束位于集熱器頂部。

1.2 集熱器的換熱量及實(shí)驗(yàn)系統(tǒng)

回轉(zhuǎn)窯內(nèi)的物料進(jìn)入回轉(zhuǎn)窯后分別經(jīng)過一系列物理化學(xué)反應(yīng)，生成水泥熟料[14]，其中各個(gè)反應(yīng)都需要一定的溫度和熱量，因此窯表面對(duì)應(yīng)的集熱器換熱量為定值[14]，通過實(shí)驗(yàn)方法與理論分析可得實(shí)際集熱器對(duì)應(yīng)窯表面的散熱量及溫度，見表 1。表中Tw1為對(duì)應(yīng)所需的窯內(nèi)壁溫，T為對(duì)應(yīng)筒體表面溫度，Q為集熱器換熱量。

圖1　回轉(zhuǎn)窯表面集熱器的結(jié)構(gòu)Fig.1 Schematic diagram of novel heat recovery exchanger

1—rotary kiln; 2—water tubes; 3—shell; 4—collecting pipes; 5—inlet header of water tubes; 6—outlet header of water tubes; 7—inlet header of coil pipes; 8—coil pipes; 9—outlet header of coil pipes

表1　集熱器對(duì)應(yīng)區(qū)域所需窯內(nèi)壁溫度Tw1、筒體外表面溫度T以及對(duì)應(yīng)集熱器所需散熱量Q[14]Table 1 Inner wall temperature, shell temperature and heat loss of measurement region

本文對(duì)集熱器進(jìn)行了傳熱性能實(shí)驗(yàn)，并在此基礎(chǔ)上進(jìn)行數(shù)值模擬計(jì)算，分析了集熱器輻射與對(duì)流換熱面的換熱特性。實(shí)驗(yàn)系統(tǒng)如圖2所示，冷卻水通過給水泵供給進(jìn)入9個(gè)支路的集熱器系統(tǒng)中，吸收回轉(zhuǎn)窯表面散熱，再匯合進(jìn)入發(fā)電系統(tǒng)中。各個(gè)支路安裝一個(gè)手動(dòng)閥門用于調(diào)節(jié)流量。

圖2　集熱器的系統(tǒng)總圖Fig.2 Schematic diagram of heat recovery system

實(shí)驗(yàn)測(cè)定了表1中集熱器的入口溫度、對(duì)應(yīng)窯表面溫度以及出口溫度，集熱器入口安裝有 UH50熱量表，熱量表測(cè)量流量及入口、出口工質(zhì)溫度，其量程分別為0.6～60 m3·h-1和5～130℃，精度等級(jí)為2級(jí)。窯表面溫度使用Testo835-T1紅外測(cè)溫儀測(cè)量，量程為-50～600℃，精度為±0.5℃。

1.3 實(shí)驗(yàn)誤差分析

管程傳熱系數(shù)由實(shí)驗(yàn)測(cè)量數(shù)據(jù)計(jì)算得出，因此需要對(duì)其相關(guān)量進(jìn)行不確定度分析，確定實(shí)驗(yàn)精度。假設(shè)函數(shù)Z含有n個(gè)相互獨(dú)立變量x1, x2, …, xn，Wx1, Wx2, …, Wxn為對(duì)應(yīng)n個(gè)變量的不確定度。則函數(shù)Z的不確定度為[15]

集熱器換熱量Q的不確定度WQ為

傳熱面積A的不確定度WA為

對(duì)數(shù)平均溫差ΔT的不確定度WΔT為

總傳熱系數(shù)k的不確定度Wk為

殼程傳熱系數(shù)hs的不確定度Whs為

管程傳熱系數(shù)ht的不確定度Wht為

以上求得的不確定度均為最大值，因此集熱器實(shí)驗(yàn)的總傳熱系數(shù)的不確定度在3.7%以內(nèi)，管程傳熱系數(shù)的不確定度在3.8%以內(nèi)，滿足工程要求。

1.4 實(shí)驗(yàn)結(jié)果與數(shù)值模擬分析

由1.2節(jié)中的實(shí)驗(yàn)測(cè)量結(jié)果計(jì)算得到集熱器管程的傳熱系數(shù)，將其與數(shù)值模擬結(jié)果對(duì)比，如圖 3和圖4所示。

由圖3、圖4可知，管程傳熱系數(shù)及其出口溫度的模擬值與實(shí)驗(yàn)值的偏差在 5%以內(nèi)，驗(yàn)證了數(shù)值計(jì)算的可靠性。

為探究集熱器的輻射與對(duì)流換熱面的換熱特性，對(duì)其換熱面的換熱量分別進(jìn)行了數(shù)值模擬，如圖5所示。結(jié)果表明，隨著空氣流速的增加，集熱器的水冷壁管排的輻射換熱量呈下降的趨勢(shì)，但其變化率在 5%以內(nèi)；而蛇形管側(cè)的對(duì)流換熱量隨空氣流速增加而增加，因此，集熱器內(nèi)空氣流速的變化對(duì)水冷壁側(cè)管束的換熱量無影響。輻射換熱面與對(duì)流換熱面的結(jié)構(gòu)可以分別進(jìn)行優(yōu)化。

圖3　管程傳熱系數(shù)的實(shí)驗(yàn)值與數(shù)值模擬值比較Fig.3 Comparison of experimental and numerical results of tube-side heat transfer coefficients

圖4　管程出口溫度的實(shí)驗(yàn)值與模擬值比較Fig.4 Comparison of water inlet and outlet temperatures with different m

圖5　對(duì)流與輻射換熱面的換熱量隨空氣流速的變化Fig.5 Heat transfer rates of water tubes and coiled pipes

2　集熱器的結(jié)構(gòu)優(yōu)化

2.1 集熱器傳熱的理論分析

對(duì)于前述所示的集熱器而言，分為對(duì)流換熱面?zhèn)扰c輻射換熱面?zhèn)?，分別對(duì)應(yīng)水冷壁管和蛇形管，對(duì)其進(jìn)行傳熱分析。圖6所示為輻射換熱面的結(jié)構(gòu)設(shè)計(jì)示意圖。

圖6　集熱器輻射換熱面管束結(jié)構(gòu)Fig.6 Structure diagram of water tubes

傳熱方程與能量守恒方程為

其中，輻射換熱對(duì)數(shù)平均溫差為

忽略換熱管導(dǎo)熱熱阻與污垢熱阻，則換熱器的傳熱系數(shù)為[16]

式中，hi和ho分別為管內(nèi)、管外傳熱系數(shù)，do和di分別為管外徑和管內(nèi)徑。其中do=di+2t，t為管厚。

以管外面積為基準(zhǔn)的換熱面積為

查文獻(xiàn)[16]得管內(nèi)傳熱系數(shù)公式為

輻射換熱面管外傳熱系數(shù)ho1

式中，ε為換熱系統(tǒng)黑度；C0為黑體輻射常數(shù)；F為窯表面對(duì)管束的角系數(shù)；T為窯表面溫度；Ti為管內(nèi)冷卻工質(zhì)平均溫度，Ti=(Tin+Tout)/2。

圖7　對(duì)流換熱面的結(jié)構(gòu)Fig.7 Structure diagram of coiled pipes

圖7所示為集熱器的對(duì)流換熱面結(jié)構(gòu)設(shè)計(jì)，其結(jié)構(gòu)為順排管束。其對(duì)應(yīng)能量方程和傳熱方程如下所示

其中，對(duì)流換熱的對(duì)數(shù)平均溫差為

對(duì)流換熱面的管內(nèi)傳熱系數(shù)hi2同式(13)所示，管外傳熱系數(shù)ho2如下[17]

對(duì)于輻射及對(duì)流換熱面而言，管內(nèi)流動(dòng)壓降公式相同，對(duì)流換熱面外側(cè)的空氣壓降相比管內(nèi)壓降較小，因此可以忽略。管內(nèi)Reynolds數(shù)

管內(nèi)流速

式中，Npass為管程數(shù)；μ為中心流體黏度；ρi為管內(nèi)流體密度。

管程的壓降計(jì)算公式為

其中，∑K為局部阻力系數(shù)之和，莫迪圓管摩擦系數(shù)為[16]

由于集熱器管束安裝尺寸限制，對(duì)流換熱面與輻射換熱面存在幾何結(jié)構(gòu)約束，對(duì)于輻射和對(duì)流換熱面，有① 管子數(shù):20

2.2 優(yōu)化方案

冷流體入口流量mi、入口溫度Tin和換熱量Q給定，優(yōu)化問題可描述為:（1）在給定集熱器進(jìn)出口壓降Δp的情況下，使其換熱面積A最小；（2）優(yōu)化進(jìn)出口壓降 Δp使其達(dá)到最小；（3）集熱器傳熱與流動(dòng)的改進(jìn)熵產(chǎn)數(shù) Ns最小。本文采用 Matlab遺傳算法工具箱求解上述問題，遺傳算法能夠自適應(yīng)地控制搜索過程，在潛在的解集中求得全局最優(yōu)解[18]。

2.2.1 壓降給定下的集熱器優(yōu)化 在給定進(jìn)出口壓降的情況下，即

優(yōu)化目標(biāo)為集熱器所需換熱面積A最小，即式(12)最小，設(shè)計(jì)參數(shù)為管子數(shù)n、管長(zhǎng)l以及管徑d。約束條件的方程有式(8)～式(11)、式(15)～式(18)及對(duì)應(yīng)傳熱系數(shù)方程。

2.2.2 壓降為目標(biāo)函數(shù)的集熱器優(yōu)化 集熱器管束不同流動(dòng)狀態(tài)下，優(yōu)化目標(biāo)為集熱器進(jìn)出口壓降Δp最小，即式(22)最小，設(shè)計(jì)參數(shù)與約束條件與2.2.1節(jié)相同。

2.2.3 以改進(jìn)熵產(chǎn)數(shù)為目標(biāo)的集熱器優(yōu)化 由于使用 Bejan[19]提出的熵產(chǎn)數(shù)進(jìn)行換熱器優(yōu)化會(huì)產(chǎn)生熵產(chǎn)悖論，因此本文采用Hesselgreaves[20]提出的改進(jìn)熵產(chǎn)數(shù)進(jìn)行優(yōu)化來避免上述問題。其定義式為

式中，Ns,ΔT、Ns,Δp分別為傳熱和流動(dòng)的改進(jìn)熵產(chǎn)數(shù)；Sgen,ΔT、Sgen,Δp分別為傳熱和流動(dòng)引起的熵產(chǎn)，其定義式分別為

優(yōu)化目標(biāo)為式(19)最小，設(shè)計(jì)參數(shù)與約束條件與2.2.1節(jié)相同。

3　優(yōu)化結(jié)果及討論

3.1 集熱器輻射換熱面的優(yōu)化

對(duì)于前述的集熱器的輻射換熱面，取入口流量m1=0.85 kg·s-1，入口水溫Tw,i1=30℃，比定壓熱容cp=4174 J·kg-1·K-1，輻射換熱量Q1=53433 W，輻射換熱面的實(shí)際結(jié)構(gòu)如文獻(xiàn)[14]中實(shí)驗(yàn)確定得出，見表2。優(yōu)化后的輻射換熱面中的管程N(yùn)pass,1=6，局部阻力系數(shù)見表3[21]，輻射換熱面的局部阻力系數(shù)∑K=4.2。

表2　集熱器輻射換熱面的實(shí)際結(jié)構(gòu)形式Table 2 Structural parameters of actual heat recovery exchanger

表3　局部阻力系數(shù)值Table 3 Values of minor loss coefficients

按實(shí)際輻射換熱面結(jié)構(gòu)形式及運(yùn)行工況求得其進(jìn)出口壓降 Δp為 13.52 Pa。令換熱系統(tǒng)黑度ε=0.8，窯表面對(duì)管束的角系數(shù)F=0.82，黑體輻射常數(shù)C0=5.67 W·m-2·K-4，將上述已知量、熱源溫度和工質(zhì)進(jìn)出口溫度代入式(14)中得管外傳熱系數(shù)h0=7.41 W·m-2·K-1。運(yùn)用遺傳算法工具箱對(duì)2.2.1節(jié)中的問題進(jìn)行求解，優(yōu)化后的集熱器結(jié)構(gòu)見表4。

比較表2、表4的集熱器結(jié)構(gòu)可知，優(yōu)化后集熱器所需換熱面積減小15%，提高了其換熱性能。

表4　優(yōu)化后的集熱器結(jié)構(gòu)參數(shù)Table 4 Structural parameters after optimization

對(duì)于輻射換熱面，求解2.2.2節(jié)中的優(yōu)化問題，優(yōu)化結(jié)果見表 5。優(yōu)化后進(jìn)出口壓降相比優(yōu)化前有所降低。

表5　優(yōu)化前后的壓降及結(jié)構(gòu)參數(shù)Table 5 Pressure drop and structural parameters of heat recovery exchanger

對(duì) 2.2.3節(jié)中的問題進(jìn)行求解，優(yōu)化輻射換熱面的結(jié)構(gòu)，降低其傳熱與流動(dòng)不可逆損失。結(jié)果見表6。

表6　優(yōu)化前后的改進(jìn)熵產(chǎn)數(shù)Table 6 Modified entropy numbers of heat recovery exchanger

表中 1、2分別代表優(yōu)化前和優(yōu)化后的結(jié)構(gòu)參數(shù)及熵產(chǎn)數(shù)，優(yōu)化后壓降引起的改進(jìn)熵產(chǎn)數(shù)Ns,Δp相比優(yōu)化前明顯降低，而傳熱引起的熵產(chǎn)數(shù) Ns,ΔT無變化，因此優(yōu)化后降低了輻射換熱面的流動(dòng)不可逆損失。

3.2 集熱器對(duì)流換熱面的優(yōu)化

對(duì)于集熱器的對(duì)流換熱面，取管內(nèi)入口流量m2=0.85 kg·s-1，熱空氣流速vair∈[0.5,3]，冷卻水溫度 Tw,i2=40℃，熱空氣出口溫度 Ta,o=50℃，對(duì)流換熱量 Q2=17811 W，優(yōu)化中換熱面的管程N(yùn)pass,2=14，其局部阻力系數(shù)見表3。對(duì)2.2.1節(jié)的問題進(jìn)行優(yōu)化，如圖8所示，圖中A1和A2分別表示優(yōu)化前后對(duì)流換熱面所需面積。結(jié)果表明，隨熱空氣流速 vair的增加，優(yōu)化前后所需換熱面積減小，優(yōu)化前相比優(yōu)化后的換熱面積減少20%左右。

給定對(duì)流換熱面的換熱面積，在變化 vair的情況下，求解2.2.2節(jié)的優(yōu)化問題，如圖9所示。圖中Δp1和Δp2分別代表優(yōu)化前后的對(duì)流換熱面的壓降。結(jié)果表明，優(yōu)化后進(jìn)出口壓降相比優(yōu)化前明顯減少，有效降低了泵功的功耗。另外，隨著空氣流速的增加，優(yōu)化后的壓降呈降低的趨勢(shì)。

圖8　優(yōu)化前后所需換熱面積對(duì)比Fig.8 Required heat transfer area with different vair

圖9　優(yōu)化前后壓降隨空氣流速的變化Fig.9 Pressure drop with different vair

圖10　優(yōu)化前后傳熱熵產(chǎn)數(shù)與壓降熵產(chǎn)數(shù)的對(duì)比Fig.10 Modified entropy numbers with different vair

針對(duì)對(duì)流換熱面，以改進(jìn)熵產(chǎn)數(shù)為優(yōu)化目標(biāo)，求解2.2.3節(jié)中的問題，結(jié)果如圖10所示。圖中Ns,ΔT1和 Ns,ΔT2分別代表優(yōu)化前后傳熱的改進(jìn)熵產(chǎn)數(shù)；Ns,Δp1和Ns,Δp2分別代表優(yōu)化前后壓降的改進(jìn)熵產(chǎn)數(shù)。優(yōu)化后的傳熱熵產(chǎn)數(shù)Ns,ΔT2相比優(yōu)化前無明顯變化，而壓降熵產(chǎn)數(shù) Ns,Δp2比優(yōu)化前明顯減少。因此，優(yōu)化后的結(jié)構(gòu)降低了流動(dòng)的不可逆性。

4　結(jié) 論

（1）以回轉(zhuǎn)窯余熱回收用新型集熱器為研究對(duì)象，通過實(shí)驗(yàn)測(cè)量與數(shù)值模擬的方法確定了其傳熱性能。數(shù)值模擬計(jì)算其輻射和對(duì)流換熱面的換熱量，發(fā)現(xiàn)空氣流速在一定范圍內(nèi)變化僅對(duì)對(duì)流換熱面的傳熱性能有影響，不會(huì)影響輻射換熱面的傳熱性能，因而上述兩種換熱面可以分別進(jìn)行結(jié)構(gòu)優(yōu)化。

（2）建立了集熱器的換熱量與結(jié)構(gòu)參數(shù)之間的數(shù)學(xué)模型，分別以換熱面積、進(jìn)出口壓降和改進(jìn)熵產(chǎn)數(shù)為目標(biāo)函數(shù)，運(yùn)用遺傳算法對(duì)集熱器的輻射換熱面和對(duì)流換熱面分別進(jìn)行了結(jié)構(gòu)優(yōu)化。結(jié)果表明，優(yōu)化后集熱器的輻射換熱所需的換熱面積減小15%，對(duì)流換熱所需的換熱面積減少20%左右。優(yōu)化后兩種形式的換熱面所消耗的泵功明顯減少。優(yōu)化后集熱器的傳熱不可逆損失與優(yōu)化前無明顯變化，而流動(dòng)不可逆損失得到明顯降低。

符 號(hào) 說 明

A ——換熱器面積，m2

C0——黑體輻射常數(shù)，W·m-2·K-4

cp——比定壓熱容，J·kg-1·K-1

d ——管徑，m

F ——窯表面對(duì)管束的角系數(shù)

f ——莫迪圓管摩擦系數(shù)

h ——管內(nèi)外傳熱系數(shù)，W·m-2·K-1

K ——局部阻力系數(shù)

k ——換熱器傳熱系數(shù)，W·m-2·K-1

l ——管長(zhǎng)，m

M, m ——流量，kg·s-1

Npass——管程數(shù)

Ns——改進(jìn)的熵產(chǎn)數(shù)

n ——管子數(shù)

Pr ——Prandtl數(shù)

Q ——換熱量，W

R ——熱阻，K·W-1

Re ——Reynolds數(shù)

r ——半徑，m

S ——熵產(chǎn)，W·K-1

T ——溫度，℃

ΔT ——對(duì)數(shù)平均溫差，℃

t ——管壁厚度，m

v——冷卻水流速，m·s-1

vair——空氣流速，m·s-1

W——不確定度

ε——換熱系統(tǒng)黑度

λ——熱導(dǎo)率，W·m-1·K-1

μ——黏度，kg·m-1·s-1

ρ——密度，kg·m-3

下角標(biāo)

a——空氣側(cè)

gen——熵產(chǎn)

i——管內(nèi)

i1——輻射換熱面?zhèn)热肟?/p>

i2——對(duì)流換熱面?zhèn)热肟?/p>

o——管外

o1——輻射換熱面?zhèn)瘸隹?/p>

o2——對(duì)流換熱面?zhèn)瘸隹?/p>

s——?dú)?cè)

w——水側(cè)

wl——管壁面

References

[1] ATMACA A, YUMRUTA? R. Analysis of the parameters affecting energy consumption of a rotary kiln in cement industry [J]. Appl. Therm. Eng., 2014, 66(1/2): 435-444.

[2] MUJUMDAR K S, GANESH K V, KULKARNI S B, et al. Rotary cement kiln simulator (rocks): integrated modeling of pre-heater, calciner, kiln and clinker cooler [J]. Chem. Eng. Sci., 2007, 62(9): 2590-2607.

[3] MADLOOL N A, SAIDUR R, HOSSAIN M S, et al. A critical review on energy use and savings in the cement industries [J]. Renewable and Sustainable Energy Rev., 2011, 15(4): 2042-2060.

[4] SUI X, ZHANG Y, SHAO S, et al. Exergetic life cycle assessment of cement production process with waste heat power generation [J]. Energy Convers.Manage., 2014, 88(4): 684-692.

[5] ATMACA A, YUMRUTA? R. Thermodynamic and exergoeconomic analysis of a cement plant (Ⅱ): Application [J]. Energy Conver. Manage., 2014, 79(3): 799-808.

[6] CAPUTO A C, PELAGAGGE P M, SALINI P. Performance modeling of radiant heat recovery exchangers for rotary kilns [J]. Appl. Therm. Eng., 2011, 31(14/15): 2578-2589.

[7] KARAMARKOVI? V, MARA?EVI? M, KARAMARKOVI? R, et al. Recuperator for waste heat recovery from rotary kilns [J]. Appl. Therm. Eng., 2013, 54(2): 470-480.

[8] S??üT Z, OKTAY Z, KARAKO H. Mathematical modeling of heat recovery from a rotary kiln [J]. Appl. Therm. Eng., 2010, 30(8/9): 817-825.

[9] ENGIN T, ARI V. Energy auditing and recovery for dry type cement rotary kiln systems –– a case study [J]. Energy Conver. Manage., 2005, 46(4): 551-562.

[10] 馬萬龍. 水泥回轉(zhuǎn)窯筒體表面余熱回收利用技術(shù)研究 [D]. 大連:大連理工大學(xué), 2012.

MA W L. Cement rotary kiln surface waste heat recovery andutilization technology research [D]. Dalian: Dalian University of Technology, 2012.

[11] 程林, 杜文靜. 正六邊形水泥回轉(zhuǎn)窯表面熱量回收裝置: CN103123219A [P]. 2013-05-29.

CHENG L, DU W J. The orthohexagonal heat recovery devices on rotary kilns: CN103123219A [P]. 2013-05-29.

[12] YANG J, FAN A, LIU W, et al. Optimization of shell-and-tube heat exchangers conforming to TEMA standards with designs motivated by constructal theory [J]. Energy Conver. Manage., 2014, 78(2): 468-476.

[13] KHALIFEH SOLTAN B, SAFFAR-AVVAL M, DAMANGIR E. Minimizing capital and operating costs of shell and tube condensers using optimum baffle spacing [J]. Appl. Therm. Eng., 2004, 24(17/18): 2801-2810.

[14] YIN Q, CHEN Q, DU W J, et al. Design requirements and performance optimization of waste heat recovery systems for rotary kilns [J]. Int. J. Heat Mass Transfer, 2016, 93(2): 1-8.

[15] 曹興. 連續(xù)與搭接螺旋折流板換熱器理論分析與實(shí)驗(yàn)研究[D]. 濟(jì)南: 山東大學(xué), 2012.

CAO X. Theoretical analysis and experimental study of shell-and-tube heat exchanger with continuous helical baffles and overlapped helical baffles [D]. Jinan: Shandong University, 2012.

[16] SHAH R K. Fundamentals of Heat Exchanger Design [M]. New Jersey: John Wiley & Sons, Inc., 2003: 88-182.

[17] 楊世銘, 陶文銓. 傳熱學(xué)[M]. 北京: 高等教育出版社, 2006: 260-261.

YANG S M, TAO W Q. Heat Transfer [M]. Beijing: Higher Education Press, 2006: 260-261.

[18] 雷英杰, 張善文, 李續(xù)武. MATLAB遺傳算法工具箱及應(yīng)用[M].西安: 西安電子科技大學(xué)出版社, 2005: 1-2.

LEI Y J, ZHANG S W, LI X W. MATLAB Algorithm Toolbox and Application [M]. Xi’an: Xi’an University of Electric Science and Technology Press, 2005: 1-2.

[19] BEJAN A. The concept of irreversibility in heat exchanger design: counterflow heat exchangers for gas-to-gas application [J]. J. Heat Transfer Trans. ASME, 1977, 99(3): 374-380.

[20] HESSELGREAVES J E. Rationalisation of second law analysis of heat exchangers [J]. Int. J. Heat Mass Transfer, 2000, 43(22): 4189-4204.

[21] WHITE F M. Fluid Mechanics [M]. Rhode Island: McGraw-Hill, 1998: 325-404.

2016-01-20收到初稿，2016-03-28收到修改稿。

聯(lián)系人:程林。第一作者:殷謙（1989—），男，博士研究生。

Received date: 2016-01-20.

中圖分類號(hào):TK 124

文獻(xiàn)標(biāo)志碼:A

文章編號(hào):0438—1157（2016）07—2740—08

DOI:10.11949/j.issn.0438-1157.20160093

基金項(xiàng)目:國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目（2013CB228305）。

Corresponding author:Prof. CHENG Lin, Cheng@sdu.edu.cn supported by the National Basic Research Program of China (2013CB228305).

Experimental measurement and structure optimization of heat recovery exchangers on rotary kilns

YIN Qian, DU Wenjing, JI Xinglin, CHENG Lin
(Institute of Thermal Science and Technology, Shandong University, Jinan 250061, Shandong, China)

Abstract:A novel heat recovery exchanger installed on the rotary kiln shell surface is proposed in this paper. The heat recovery exchanger contains water tubes and coiled pipes, which function as the radiative and the convective heat transfer surfaces, respectively. Numerical studies and experimental measurements are carried out to investigate the heat transfer characteristics of the heat exchanger. Several optimization models are proposed with the heat transfer area, the pressure drop and the modified entropy generation numbers set as the objective functions. The models describe the relation between the heat transfer rates and structural parameters of the heat recovery exchanger, i.e. the tube length, the tube numbers and the tube diameter. The optimized design parameters are obtained by applying the genetic algorithm toolbox in Matlab. The results indicate that the optimized heat transfer areas of water tubes and coiled pipes are decreased by 15% and 20%, respectively. The corresponding pressure drop is significantly decreased after optimization. In the optimization process, the modified entropy generation numbers are decreased due to fluid friction, while the modified entropy generation numbers remain unchanged due to heat transfer.

Key words:energy conservation; rotary kilns; heat recovery exchanger; heat transfer; optimization