李煜，徐志明,陳健美,陳玲,劉冠麟,李玉強(qiáng),袁文華

(1.邵陽(yáng)學(xué)院 機(jī)械與能源工程系， 湖南 邵陽(yáng)，422000；2.湖南大學(xué) 機(jī)械與運(yùn)載工程學(xué)院，湖南 長(zhǎng)沙，410082；3.湖南涉外經(jīng)濟(jì)學(xué)院 機(jī)械工程學(xué)院，湖南 長(zhǎng)沙，410205；4.東北電力大學(xué) 能源與動(dòng)力工程學(xué)院，吉林 吉林，132012)

基于場(chǎng)協(xié)同原理的橫紋管抗垢性能數(shù)值研究

李煜1,2,3，徐志明4,陳健美3,陳玲4,劉冠麟3,李玉強(qiáng)3,袁文華1

(1.邵陽(yáng)學(xué)院 機(jī)械與能源工程系， 湖南 邵陽(yáng)，422000；2.湖南大學(xué) 機(jī)械與運(yùn)載工程學(xué)院，湖南 長(zhǎng)沙，410082；3.湖南涉外經(jīng)濟(jì)學(xué)院 機(jī)械工程學(xué)院，湖南 長(zhǎng)沙，410205；4.東北電力大學(xué) 能源與動(dòng)力工程學(xué)院，吉林 吉林，132012)

基于傳熱傳質(zhì)學(xué)控制方程和Kern-Seaton污垢形成理論，構(gòu)建了橫紋管內(nèi)基于表面化學(xué)反應(yīng)理論的CaCO3污垢熱阻模型，然后根據(jù)場(chǎng)協(xié)同原理，運(yùn)用Fluent軟件對(duì)橫紋管內(nèi)以一定濃度的CaCO3溶液為工質(zhì)的污垢生成過(guò)程進(jìn)行了數(shù)值模擬。通過(guò)改變壁溫、入口流速和CaCO3溶液濃度，模擬了不同種狀態(tài)下污垢生成過(guò)程。在此過(guò)程中實(shí)時(shí)計(jì)算平均協(xié)同角和平均場(chǎng)協(xié)同數(shù)并對(duì)比分析污垢沉積率、剝蝕率和污垢凈存速率來(lái)評(píng)價(jià)不同管型橫紋管的抗垢性能。最后將模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了驗(yàn)證，驗(yàn)證了模擬的正確性。結(jié)果表明：較其它兩管型，10-8-1型橫紋管的抗垢性能較好且其更適于管內(nèi)流速跨度較大的工況，達(dá)動(dòng)態(tài)平衡狀態(tài)時(shí)的污垢各參數(shù)值基本上不受工質(zhì)中Ca2+濃度的影響。

橫紋管；場(chǎng)協(xié)同；污垢；CaCO3

作為粗糙表面的無(wú)源強(qiáng)化換熱元件，橫紋管具有傳熱與防垢性能良好、結(jié)構(gòu)簡(jiǎn)單、制造方便、熱膨脹適應(yīng)性與防振力均強(qiáng)等優(yōu)點(diǎn)，現(xiàn)已廣泛應(yīng)用于能源動(dòng)力的諸多相關(guān)領(lǐng)域[1]。

近年來(lái)，已有不少關(guān)于橫紋管換熱等性能方面的研究。張仲彬等[2]實(shí)驗(yàn)研究了橫紋管與光管的傳熱與流阻特性，擬合了其傳熱與流阻實(shí)驗(yàn)關(guān)聯(lián)式，探析了傳熱機(jī)理并仿真對(duì)比分析了湍流強(qiáng)度和局部Nu數(shù)等關(guān)鍵強(qiáng)化傳熱參數(shù)的變化規(guī)律。李洪亮等[3]研究了以 LiBr水溶液為工質(zhì)的吸收器內(nèi)橫紋管的強(qiáng)化管外降膜吸收過(guò)程，結(jié)果表明橫紋管比光管的傳熱與傳質(zhì)效果均要好。林緯[4]以實(shí)驗(yàn)和仿真方法研究了脈動(dòng)流對(duì)橫紋管內(nèi)外傳熱流動(dòng)的影響，發(fā)現(xiàn)脈動(dòng)流對(duì)光管換熱無(wú)影響而對(duì)橫紋管換熱有強(qiáng)化作用。陳聰?shù)萚5,6]實(shí)驗(yàn)探析了熔鹽為工質(zhì)的不同結(jié)構(gòu)參數(shù)的橫紋管強(qiáng)化傳熱機(jī)理，擬合了其傳熱關(guān)聯(lián)式并對(duì)比光管發(fā)現(xiàn)9mm節(jié)距橫紋管的傳熱強(qiáng)化比達(dá)1.6。文玉良等[7]實(shí)驗(yàn)研究了橫紋管內(nèi)熔融鹽為介質(zhì)的傳熱過(guò)程，對(duì)比光管分析了熱流密度和溫度對(duì)其傳熱的影響，分析了強(qiáng)化傳熱機(jī)理，結(jié)果表明橫紋管比光管傳熱更佳。沈向陽(yáng)等[8]對(duì)比實(shí)驗(yàn)研究了熔鹽為工質(zhì)的螺旋槽管與橫紋管的傳熱特性，發(fā)現(xiàn)兩者均能有效提升傳熱系數(shù)且傳熱隨槽深增加而更佳；前者比后者的高黏度熔鹽傳熱好且兩者的強(qiáng)化傳熱倍數(shù)隨Re數(shù)增加而緩降；而前者對(duì)低黏度熔鹽傳熱強(qiáng)化倍數(shù)隨Re數(shù)增加而緩降，后者則隨Re數(shù)增加而緩升；兩者對(duì)低黏度熔鹽傳熱強(qiáng)化強(qiáng)于高黏度的，而后者更適于低黏度熔鹽高Re數(shù)傳熱。Lu等[9]研究了橫紋管內(nèi)高溫熔鹽的對(duì)流傳熱，擬合了傳熱關(guān)聯(lián)式，發(fā)現(xiàn)橫紋管能明顯強(qiáng)化高溫熔鹽的傳熱且傳熱隨槽深和Re數(shù)增加而加強(qiáng)，當(dāng)溫度很高時(shí)傳熱惡化且發(fā)生在較低的傳熱系數(shù)處。王樹(shù)濤[10]實(shí)驗(yàn)研究了混合硝酸鹽為工質(zhì)的三組不同參數(shù)的橫紋管替換光管的改造傳熱系統(tǒng)的傳熱與阻力特性，擬合了傳熱關(guān)聯(lián)式并評(píng)選了最優(yōu)強(qiáng)化管，結(jié)果表明：橫紋管與光管隨Re數(shù)的傳熱變化規(guī)律大體相似，但前者性能優(yōu)于后者；前者阻力系數(shù)比后者明顯增加且隨Re數(shù)增加而減小、隨節(jié)距P減小而增大；P=9的橫紋管傳熱最佳。Akyildiz等[11]研究了橫紋管內(nèi)的能量損失和強(qiáng)化傳熱，發(fā)現(xiàn)隨著橫紋節(jié)數(shù)的增加，其阻力系數(shù)會(huì)隨之減小，傳熱效率將會(huì)提升。閔亞光[12]通過(guò)對(duì)比光管實(shí)驗(yàn)研究了糖廠(chǎng)蒸發(fā)系統(tǒng)橫紋管強(qiáng)化傳熱和抗垢性能，結(jié)果發(fā)現(xiàn)橫紋管有很好的強(qiáng)化傳熱與抗垢性能且其積垢亦易清除。徐志明等[13]實(shí)驗(yàn)研究了橫紋管的傳熱性能，擬合出其強(qiáng)制對(duì)流換熱關(guān)聯(lián)式并對(duì)比實(shí)驗(yàn)研究了同工況下的800mg/L 人工硬水為工質(zhì)的橫紋管與光管的污垢特性，結(jié)果表明橫紋管傳熱較好且阻垢亦較佳。

綜上所述，關(guān)于橫紋管強(qiáng)化傳熱的研究不少，而關(guān)于其污垢方面的研究卻較少且實(shí)驗(yàn)研究較多而數(shù)值研究較少，因此筆者基于場(chǎng)協(xié)同原理[14]開(kāi)展橫紋管抗垢性能的數(shù)值研究是有創(chuàng)新意義的，文中建構(gòu)了以Fluent軟件來(lái)實(shí)現(xiàn)的基于表面化學(xué)反應(yīng)理論的CaCO3污垢熱阻模型并基于場(chǎng)協(xié)同原理分析了不同橫紋管的抗垢性能，通過(guò)對(duì)比仿真與實(shí)驗(yàn)結(jié)果且經(jīng)誤差分析論證了該模型的正確性。

1 數(shù)學(xué)物理模型

1.1 物理模型

由于橫紋管與光管均為沿中軸線(xiàn)嚴(yán)格軸對(duì)稱(chēng)的圖形，故在Fluent仿真過(guò)程中，將實(shí)際三維模型簡(jiǎn)化為二維模型。文中所采用的橫紋管和光管物理模型如圖1、2所示。圓管內(nèi)徑d=22mm，管長(zhǎng)L=2.23m，根據(jù)尼古拉茲實(shí)驗(yàn)結(jié)論，管長(zhǎng)與內(nèi)徑滿(mǎn)足如下的關(guān)系式時(shí)，便可不計(jì)湍流狀態(tài)下管入口段的影響。因此，文中忽略了管入口段對(duì)仿真的影響。文中實(shí)驗(yàn)所用橫紋管與光管的材質(zhì)與幾何參數(shù)如表1所列。其中，橫紋管和光管的不銹鋼材質(zhì)類(lèi)型均為304，幾何尺寸L1、L2、R、e參數(shù)的含義分別為每一節(jié)橫紋管的直段管長(zhǎng)、橫紋段管長(zhǎng)、管內(nèi)的半徑、橫紋槽深。

L=(25～40)d

(1)

圖1 橫紋管結(jié)構(gòu)示意圖Fig.1 Diagram of atransversally corrugated tube

圖2 光管結(jié)構(gòu)示意圖Fig.2 Geometry diagram of a plain tube

管型材質(zhì)幾何尺寸(mm)L1L2Re橫紋管不銹鋼304178112光管不銹鋼30411

1.2 數(shù)學(xué)模型

1.2.1 控制方程

為了仿真整個(gè)控制區(qū)域的傳熱和污垢生成情況，需建立如下控制方程組：

連續(xù)性方程

(2)

動(dòng)量方程

r方向

(3a)

z方向

(3b)

其中

能量方程

(4)

傳質(zhì)方程

(5)

式中，r、z分別指柱坐標(biāo)下r(管徑方向)，z(管長(zhǎng))坐標(biāo)方向；uz，ur分別是柱坐標(biāo)下z、r方向的速度分量，m/s；ρ是主流流體的密度，kg/m3；p是壓力，Pa；fr、fz是r、z方向的質(zhì)量力，m/s2；μ是主流流體的動(dòng)力粘度，Pa·s；Sr、Sz分別是控制方程沿r、z方向的源項(xiàng)，kg/m2·s2；T是溫度，K；a是熱擴(kuò)散率，m2/s；cf、cF分別是主流濃度、熱邊界層濃度，kg/m3；D是質(zhì)擴(kuò)散系數(shù)，m2/s；hm是傳質(zhì)系數(shù)，m/s；δT是熱邊界層厚度，m。

以上是一個(gè)封閉的數(shù)學(xué)模型，給定初始條件和邊界條件即可求解。本文模擬對(duì)象涉及湍流流動(dòng)，湍流流動(dòng)會(huì)引起流體動(dòng)量、能量和濃度的變化，并且導(dǎo)致數(shù)量的波動(dòng)。這種波動(dòng)具有小尺度和高頻率的特點(diǎn)，若用計(jì)算機(jī)直接模擬，會(huì)對(duì)計(jì)算機(jī)計(jì)算速度及內(nèi)存等提出很高的要求，實(shí)際情況難以滿(mǎn)足。因此，F(xiàn)luent軟件提供了湍流模型來(lái)解決這一問(wèn)題。

1.2.2 湍流模型

由于在本模擬研究過(guò)程中，流動(dòng)被假定為完全湍流，忽略了分子粘性的影響，故本文采用標(biāo)準(zhǔn)的湍動(dòng)能-湍流耗散率k-ε模型，可簡(jiǎn)化為

(6)

(7)

1.2.3 近壁面處理

標(biāo)準(zhǔn)k-ε模型僅適用于遠(yuǎn)離壁面的湍流核心區(qū)域，為使該模型適于壁面邊界層處的流動(dòng)，有必要對(duì)近壁面區(qū)做些特殊處理。Fluent中常用近壁面處理方法有兩種:壁面函數(shù)法和近壁面模型法。前者在低雷諾數(shù)狀態(tài)下模擬效果不太理想，為使模擬更加精準(zhǔn)，對(duì)于基于ε方程的模型，文中采用后者中的增強(qiáng)壁面函數(shù)法。

增強(qiáng)壁面函數(shù)法采用兩層區(qū)域模型，將湍流雷諾數(shù)Rey的值將流動(dòng)區(qū)域劃分為粘性影響區(qū)和完全湍流區(qū)，在兩區(qū)域分別用不同的模型來(lái)計(jì)算。

(8)

Rey小于200的區(qū)域?qū)儆诮谡承杂绊憛^(qū)，采用Wolfstein一方程模型對(duì)其進(jìn)行計(jì)算，Rey大于200的區(qū)域?qū)儆谕耆牧鲄^(qū)，用標(biāo)準(zhǔn)k-ε模型計(jì)算。

1.2.4 邊界條件

入口以速度為邊界條件，u=0.5～3m/s，v=0，并在湍流模型中給出湍流強(qiáng)度和當(dāng)量直徑。湍流強(qiáng)度定義為：

I=0.16(ReDH)-1/8

(9)

1.3 污垢模型

建構(gòu)污垢模型前，首先作如下假定來(lái)對(duì)模型作適當(dāng)簡(jiǎn)化：

(1)只研究析晶污垢CaCO3的生成過(guò)程，忽略其它類(lèi)型污垢，亦忽略CaCO3生成過(guò)程的逆向反應(yīng)；

(2)污垢特性參數(shù)各向同性、分布均勻；

(3)不計(jì)污垢表面粗糙度的影響；

(4)不計(jì)耗散熱和化學(xué)反應(yīng)熱；

(5)不計(jì)分子擴(kuò)散引發(fā)的能量傳遞；

(6)不計(jì)流體在污垢生成過(guò)程中的物性變化。

1.3.1 污垢沉積模型

假設(shè)該反應(yīng)僅在管壁上發(fā)生，主流體中不生成析晶污垢。從主流體到熱邊界層只有傳質(zhì)，Ca2+、CO32-不斷從主流體輸運(yùn)到熱邊界層內(nèi)，補(bǔ)充熱邊界層內(nèi)由于析晶污垢生成而失去的離子，使管壁附近的熱邊界層內(nèi)CaCO3濃度相等。

壁面處CaCO3污垢沉積率為

(10)

(11)

(12)

Δc=cf-cs

(13)

cs=98.85714-1.71071T+9.82×10-3T2

(14)式中，傳質(zhì)系數(shù)hm，m/s；kR是表面反應(yīng)速率常數(shù)；cp是比定壓比熱容，kJ/(kg·K)；ρ是污垢密度，kg/m3；μ是動(dòng)力粘度，Pa·s；KB是Boltzman常數(shù)；E是活化能，J/mol；α是熱擴(kuò)散系數(shù)，m2/s；rd是溶質(zhì)半徑，m；cf是主流濃度，kg/m3；cs是飽和濃度，kg/m3；λ是導(dǎo)熱系數(shù)，W/(m·K)；Tf、TF分別是主流平均溫度和熱邊界層溫度，K。

1.3.2 污垢剝蝕模型

由CaCO3污垢理論分析方法[15]可知，污垢層的溫差會(huì)導(dǎo)致熱應(yīng)力的產(chǎn)生，從而削弱污垢層的強(qiáng)度。當(dāng)施加在污垢層上的剪切應(yīng)力比其切變強(qiáng)度強(qiáng)時(shí)，污垢從換熱面最終剝離。

污垢剝蝕率的表達(dá)式為

(15)

式中，u是主流平均速度，m/s；mf是污垢質(zhì)量，kg；Tw為壁面溫度，K；β為線(xiàn)性膨脹系數(shù)，1/K；dp為晶體粒徑，m。

1.3.3 污垢熱阻模型

污垢生成是一個(gè)動(dòng)態(tài)的過(guò)程，既有污垢不斷沉積到換熱面，同時(shí)亦不斷有污垢從換熱面上剝離，此兩現(xiàn)象均影響熱阻值，故最終監(jiān)測(cè)到的污垢熱阻值Rf是兩者的迭加[16]。Kern和Seaton提出下述常微分方程來(lái)描述該過(guò)程：

(16)

式中，ρf是污垢層平均密度，kg/m3；λf是污垢層導(dǎo)熱系數(shù)，W/(m·K)。

1.3.4 層流有限速率模型

污垢沉積模型中已提及污垢生成的化學(xué)反應(yīng)主要在壁面及壁面附近的熱邊界層內(nèi)進(jìn)行，此區(qū)域內(nèi)可基本忽略湍流脈動(dòng)對(duì)流體流動(dòng)的影響，故模擬污垢生成時(shí)需模擬壁面反應(yīng)與容積反應(yīng)且在模擬后者時(shí)選用層流有限速率模型，用其計(jì)算時(shí)不計(jì)湍流脈動(dòng)對(duì)流動(dòng)的影響，用Arrhenius公式計(jì)算每種物質(zhì)的化學(xué)源項(xiàng)。化學(xué)物質(zhì)i的化學(xué)反應(yīng)凈源項(xiàng)可通過(guò)對(duì)有其參加的NR個(gè)化學(xué)反應(yīng)的反應(yīng)源求和得到：

(17)

式中，物質(zhì)i的產(chǎn)生摩爾速率

(18)式中，Γ為第三體對(duì)反應(yīng)速率的凈影響，此處取0。

物質(zhì)i第r個(gè)反應(yīng)可表示為

(19)

式中，kf,r為反應(yīng)r的前向速率常數(shù)，用Arrhenius公式計(jì)算。

kf,r=ArTβre-Er/RT

(20)

式中，指前因子Ar取5656m4/(kg·s)，反應(yīng)活化能Er取21830J/mol，溫度指數(shù)βr取1。

2 橫紋管抗垢性能比較

眾多影響污垢沉積的因素中，最主要的為流體速度和流體中CaCO3濃度，下面以單變量的方式分別從這兩方面著手分析橫紋管的抗垢性能。此處流速分別為：1.0m/s，1.5m/s，2 m/s和2.5 m/s；CaCO3濃度分別為：0.2kg/m3，0.4kg/m3和0.8kg/m3。

2.1 污垢沉積、剝蝕的數(shù)值模擬結(jié)果

圖3是同工況下三種管型的橫紋管內(nèi)沉積率、剝蝕率及凈存速率的時(shí)變規(guī)律。圖中顯示，10-8-1(即表1中L1、L2、e的尺寸分別為10、8、1mm，其它管型依此定義)型橫紋管在整個(gè)污垢生成過(guò)程中的凈存速率遠(yuǎn)低于17-8-2型與10-12-2型的。故10-8-1型橫紋管的抗垢性能較好。接下來(lái)依據(jù)場(chǎng)協(xié)同原理具體論證這一結(jié)論。

(a)

(b)

(c)圖3 橫紋管內(nèi)沉積率、剝蝕率及凈存速率的時(shí)變規(guī)律：(a)10-8-1型；(b)10-12-2型；(c)17-8-2型Fig.3 Traces of deposition rate,removal rate and net deposition rate of fouling in transversally corrugated tubes:(a)10-8-1 type; (b)10-12-2 type; (c)17-8-2 type

2.2 污垢生成過(guò)程中各參數(shù)變化規(guī)律

(a)

(b)

(c)圖4 各參數(shù)的時(shí)變規(guī)律：(a)平均協(xié)同角；(b)平均熱邊界層厚度；(c)平均場(chǎng)協(xié)同數(shù)Fig.4 Traces of each parameter:(a)average synergy angle; (b)average thermal boundary layer thickness; (c)average field synergy number

2.3 抗垢性能隨雷諾數(shù)Re的變化規(guī)律

(a)

(b)

(c)圖5 抗垢性能隨雷諾數(shù)Re的變化規(guī)律：(a)平均協(xié)同角；(b)平均熱邊界層厚度；(c)平均場(chǎng)協(xié)同數(shù)Fig.5 Anti-fouling performance vs.Re:(a)average synergy angle; (b)average thermal boundary layer thickness; (c)average field synergy number

2.4 抗垢性能隨Ca2+濃度的變化規(guī)律

(a)

(b)

(c)圖6 抗垢性能隨Ca2+濃度的變化規(guī)律：(a)平均協(xié)同角；(b)平均熱邊界層厚度；(c)平均場(chǎng)協(xié)同數(shù)Fig.6 Anti-fouling performance vs.Ca2+ concentration:(a)average synergy angle; (b)average thermal boundary layer thickness; (c)average field synergy number

3 實(shí)驗(yàn)驗(yàn)證

3.1 實(shí)驗(yàn)原理

采用文獻(xiàn)[13]中的方法實(shí)時(shí)在線(xiàn)監(jiān)測(cè)本實(shí)驗(yàn)中的CaCO3污垢生成過(guò)程。實(shí)驗(yàn)系統(tǒng)構(gòu)成如下：冷卻水系統(tǒng)、循環(huán)水泵、管道調(diào)節(jié)閥、水浴水箱、高位水箱、低位水箱、實(shí)驗(yàn)管、電加熱器、主機(jī)、溫控器、數(shù)據(jù)采集器、溫度傳感器、數(shù)據(jù)通信卡等，實(shí)驗(yàn)系統(tǒng)如圖7所示。

(a)

(b)圖7 污垢熱阻動(dòng)態(tài)測(cè)量系統(tǒng)：(a)原理示意圖；(b)實(shí)物圖Fig.7 Dynamic measurement system of thermal resistance for fouling:(a)schematic diagram;(b)actual photograph

用循環(huán)水泵將實(shí)驗(yàn)工質(zhì)從低位水箱輸送至可控制水位高度而對(duì)整套實(shí)驗(yàn)系統(tǒng)起定壓作用的高位水箱中，當(dāng)高位水箱的水位高出預(yù)設(shè)水位線(xiàn)時(shí)，工質(zhì)便經(jīng)溢流管回流至低位水箱中，這便維持了實(shí)驗(yàn)系統(tǒng)中流速的恒定。低位水箱中安設(shè)的換熱器可維持流經(jīng)圓管的工質(zhì)的入口溫度恒定。該實(shí)驗(yàn)用流量計(jì)與閥門(mén)來(lái)調(diào)節(jié)工質(zhì)的流速以此來(lái)調(diào)節(jié)冷卻水流量，從而調(diào)控工質(zhì)溫度；用熱電阻和溫控器來(lái)調(diào)節(jié)電加熱器來(lái)維持水浴溫度恒定。數(shù)據(jù)采集器采集實(shí)驗(yàn)數(shù)據(jù)后輸入計(jì)算機(jī)并經(jīng)計(jì)算機(jī)中的軟件進(jìn)行數(shù)據(jù)處理，最終實(shí)時(shí)輸出換熱管內(nèi)的污垢熱阻值。

該實(shí)驗(yàn)?zāi)Ｐ偷奶攸c(diǎn)如下：

(1)完全從污垢熱阻的原始定義出發(fā)，沒(méi)有任何假設(shè)，應(yīng)用范圍廣；

(2)只能獲取換熱面的平均污垢熱阻而無(wú)法計(jì)算其局部污垢熱阻；

(3)需先以清潔水為工質(zhì)，測(cè)得其傳熱系數(shù)值，再以污垢水為工質(zhì)，測(cè)得其結(jié)垢狀態(tài)的傳熱系數(shù)值，期間的轉(zhuǎn)換過(guò)程要人工監(jiān)控。

(4)測(cè)得的污垢熱阻值波動(dòng)較小，故其穩(wěn)定性、可靠性均較好。

3.2 模擬結(jié)果與實(shí)驗(yàn)結(jié)果的比較

(a)

(b)圖8 污垢熱阻的實(shí)驗(yàn)與模擬對(duì)比：(a)污垢熱阻曲線(xiàn)；(b)相對(duì)誤差曲線(xiàn)Fig.8 Comparison of experiment and simulation results of thermal resistance for fouling:(a)trace of thermal resistance for fouling;(b)relative error trace

本實(shí)驗(yàn)為了驗(yàn)證數(shù)值模擬的正確性，故模擬了以上三種管在入口流速為0.37m/s的CaCO3污垢生成情況，其中17-8-2型橫紋管的如8(a)所示。將熱阻值的實(shí)驗(yàn)值與模擬值進(jìn)行對(duì)比分析并計(jì)算其相對(duì)誤差，如8(b)所示。其它兩種管的情況類(lèi)似，在此不再另外給出其污垢熱阻的實(shí)驗(yàn)與模擬對(duì)比曲線(xiàn)。由于本模型沒(méi)有考慮污垢生成過(guò)程的誘導(dǎo)期，故開(kāi)始階段誤差較大。兩者間相對(duì)誤差隨著污垢的不斷增長(zhǎng)而逐漸減小，達(dá)到誤差允許范圍之內(nèi)，最終穩(wěn)定在約20%處，這對(duì)污垢實(shí)驗(yàn)來(lái)說(shuō)是足夠的，因此驗(yàn)證了該模型的正確性。

4 結(jié)論

本文模擬了CaCO3污垢在不同管型橫紋管內(nèi)的生長(zhǎng)過(guò)程，分析評(píng)價(jià)了其污垢特性并用實(shí)驗(yàn)數(shù)據(jù)對(duì)模擬結(jié)果進(jìn)行了驗(yàn)證，驗(yàn)證了模擬的正確性，得到了如下結(jié)論：

(1)10-8-1型橫紋管的平均場(chǎng)協(xié)同數(shù)比其它兩種管型的大，故其抗垢性能較好。

(2)橫紋管的抗垢性能隨雷諾數(shù)Re的增加而減弱，較其它兩種管型，10-8-1型橫紋管的平均場(chǎng)協(xié)同數(shù)隨雷諾數(shù)Re增加而減小的速率較小，故其更適于管內(nèi)流速跨度較大的工況。

(3)污垢達(dá)動(dòng)態(tài)平衡狀態(tài)時(shí)的各參數(shù)值基本上不受工質(zhì)中Ca2+濃度的影響。

[1]黃德彬,鄧先和,邢華偉.異型管內(nèi)的水力、傳熱性能的數(shù)值模擬 [J].高?；瘜W(xué)工程學(xué)報(bào),2003,17(2):146-150.

[2]張仲彬,陳玲,徐志明，等.橫紋管強(qiáng)化傳熱實(shí)驗(yàn)研究與機(jī)制分析 [J].中國(guó)電機(jī)工程學(xué)報(bào),2010,30(增刊):101-105.

[3]李洪亮,柳坤.橫紋管強(qiáng)化吸收器的傳熱傳質(zhì)實(shí)驗(yàn)研究 [J].化學(xué)工程,2011,39(1):14-18.

[4]林緯.橫紋管脈沖流動(dòng)與壁面振動(dòng)對(duì)流傳熱特性研究 [D].武漢:武漢工程大學(xué),2011.

[5]陳聰,吳玉庭,王樹(shù)濤，等.橫紋管內(nèi)混合硝酸鹽傳熱特性的實(shí)驗(yàn)研究 [J].工程熱物理學(xué)報(bào),2013,34(5):927-930.

[6]Chen C, Wu Y, Wang S, et al. Experimental investigation on enhanced heat transfer in transversally corrugated tube with molten salt [J]. Experimental Thermal and Fluid Science, 2013, 47: 108-116.

[7]文玉良,丁靜,楊曉西，等.高溫熔鹽橫紋管傳熱特性與強(qiáng)化機(jī)理研究 [J].工程熱物理學(xué)報(bào),2010,31(1):113-115.

[8]沈向陽(yáng),陸建峰,丁靜，等.熔鹽在螺旋槽管和橫紋管內(nèi)強(qiáng)化傳熱特性 [J].工程熱物理學(xué)報(bào),2013,34(6):1149-1152.

[9]Lu J,Shen X，Ding J, et al. Convective heat transfer of high temperature molten salt in transversely grooved tube [J]. Applied Thermal Engineering, 2013, 61: 157-162.

[10]王樹(shù)濤.混合硝酸鹽橫紋管強(qiáng)化傳熱實(shí)驗(yàn)研究 [D].北京:北京工業(yè)大學(xué),2012.

[11]Akyildiz F T, Siginer D A, Khezzar L. Energy losses and heat transfer enhancement in transversally corrugated pipes [J]. International Journal of Heat and Mass Transfer, 2011, 54: 3801-3806.

[12]閔亞光.糖廠(chǎng)蒸發(fā)系統(tǒng)橫紋管強(qiáng)化傳熱和抗垢性能的研究 [J].齊齊哈爾輕工學(xué)院學(xué)報(bào),1994,10(1):19-25.

[13]徐志明,楊善讓,甘云華.橫紋管污垢性能的實(shí)驗(yàn)研究 [J].中國(guó)電機(jī)工程學(xué)報(bào),2005,25(5):159-163.

[14]過(guò)增元.場(chǎng)協(xié)同原理與強(qiáng)化傳熱新技術(shù)[M].北京:中國(guó)電力出版社,2004.

[15]Brahim F,Augustin W,Bohnet M.Numerical simulation of the fouling process [J].International Journal of Thermal Sciences,2003,42(3):323- 334.

[16]Muller-Steinhagen H,Malayeri M R,Watkinson A P.Fouling of heat exchangers-new approaches to solve an old problem [J].Heat Transfer Engineering,2005,26(1):1-4.

Numerical study on anti-fouling performance of a transversally corrugated tube based on field synergy principle

LI Yu1,2,3,XU Zhiming4,CHEN Jianmei3,CHEN Ling4,LIU Guanlin3,LI Yuqiang3,YUAN Wenhua1

(1.Department of Mechanical and Energy Engineering,Shaoyang University,Shaoyang 422000,China;2.College of Mechanical & Vehicle Engineering,Hunan University,Changsha 410082,China;3.Department of Mechanical Engineering,Hunan International Economics University,Changsha 410205,China;4.School of Energy and Power Engineering,Northeast Dianli University,Jilin 132012,China)

Based on governing equations of heat and mass transfer as well as Kern-Seaton’s Fouling Theory,a thermal resistance model of CaCO3fouling based on the Surface Chemical Reaction Theory was set up.Then according to the Field Synergy Principle,the fouling process occurring in a transversally corrugated tube where the working fluid was certain concentration of CaCO3solution was numerically simulated.The fouling processes in various states were simulated via changing the wall temperature,inlet flow velocity and concentration of CaCO3solution.During the process,in order to evaluate the anti-fouling performance of different transversally corrugated tube patterns in the fouling state,the average field synergy angle and average field synergy number were calculated as well as the fouling deposition rate,removal rate and net fouling deposition rate of all tubes were compared and analyzed.Finally,simulation results were verified by experiments,which verified the correctness of simulation.The results show that the anti-fouling performance of the 10-8-1 type transversally corrugated tube is better than the other two tubes as well as is more suitable for operation conditions with large flow velocity range,besides,each parameter is hardly affected by the concentration of Ca2+in the working fluid when fouling is up to a dynamic equilibrium state.

transversally corrugated tube; field synergy; fouling; CaCO3

1672-7010(2017)02-0038-10

2017-02-26

國(guó)家自然科學(xué)基金資助項(xiàng)目(51276056,51176045,91541121); 湖南省教育廳科研項(xiàng)目(13C492)

李煜(1984-)，男，江西萍鄉(xiāng)人，講師，在讀博士研究生，從事工程傳熱傳質(zhì)和燃燒學(xué)研究，E-mail:augustus168@163.com

徐志明(1959-)，男，吉林九臺(tái)人，教授，博士，博士生導(dǎo)師，從事強(qiáng)化傳熱與換熱設(shè)備污垢研究，E-mail:xuzm@mail.nedu.edu.cn

TK124

A