基于CFD UV反應(yīng)器的優(yōu)化設(shè)計與仿真

周國祥， 肖英杰， 魯 陽， 王易晨， 陸舒楠

(上海海事大學(xué) a. 商船學(xué)院； b. 航運仿真技術(shù)教育部工程研究中心， 上海 201306)

2017-04-28

周國祥(1991—)，男，重慶忠縣人，碩士，研究方向為載運工具運用工程。E-mail：zhougqmaritime@sina.com

1000-4653(2017)03-0093-05

基于CFDUV反應(yīng)器的優(yōu)化設(shè)計與仿真

周國祥a,b， 肖英杰a,b， 魯 陽a,b， 王易晨a,b， 陸舒楠a,b

(上海海事大學(xué) a. 商船學(xué)院； b. 航運仿真技術(shù)教育部工程研究中心， 上海 201306)

為提高紫外線船舶壓載水器的處理效率，在現(xiàn)有紫外線壓載水處理技術(shù)的基礎(chǔ)上，通過優(yōu)化反應(yīng)器進(jìn)出口設(shè)計和紫外線燈管的布置方案，使水流在反應(yīng)器內(nèi)部的流場得到改善，紫外線輻照強(qiáng)度在反應(yīng)器截面上的分布更加均勻，進(jìn)而提高反應(yīng)器的處理效率。通過計算流體力學(xué)(Computational Fluid Dynamics, CFD)仿真，驗證該方案的可行性。結(jié)合仿真結(jié)果，提出進(jìn)一步提高紫外線反應(yīng)器處理效率的建議。

紫外線壓載水處理； 計算流體力學(xué)； 裝置優(yōu)化； 流場改善

船舶壓載水在保證船舶安全航行方面發(fā)揮著重要作用。隨著船舶的大型化發(fā)展，單船壓載艙容量及單船壓載水排放對港口和沿海水域的影響越來越大。物種轉(zhuǎn)移很可能會造成生物入侵，帶來極大的危害。[1]為控制和防止船舶壓載水傳播有害水生物及病原體，國際海事組織(International Maritime Organization, IMO)于 2004 年制定了《船舶壓載水及沉積物控制和管理國際公約》(以下簡稱《2004 壓載水管理公約》)。該公約對壓載水的排放管理及排放標(biāo)準(zhǔn)提出了要求，要求到 2017 年所有船舶都安裝壓載水處理設(shè)備，否則公約生效后不能駛?cè)隝MO 成員國港口，違反公約將受到處分。[2]

近年來，世界各國加快了研究各類壓載水處理設(shè)備的步伐。在以往的紫外線壓載水處理器研究中，反應(yīng)器內(nèi)部紫外線強(qiáng)度容易受燈管布置的影響，且其設(shè)計不具備延展性。對此，在已有研究的基礎(chǔ)上提出具備延展性的燈管布置方案，使反應(yīng)器的處理能力和處理效率不再受燈管排列方式的限制。

1 IMO關(guān)于壓載水處理設(shè)備排放的標(biāo)準(zhǔn)

IMO船舶壓載水排放標(biāo)準(zhǔn)見表1。

IMO 有關(guān)標(biāo)準(zhǔn)[3]規(guī)定試驗認(rèn)證用壓載水處理系統(tǒng)的壓載水處理量不得低于200 m3/h。

2 紫外線殺菌的理論基礎(chǔ)

2.1紫外線殺菌原理

根據(jù)相關(guān)研究[4-5]，在不同溫度下，DNA和RNA吸收紫外線光譜的范圍為240～280 nm，其中：波長為253.7 nm的紫外線吸收效果最好；波長為185 nm的紫外線能產(chǎn)生臭氧，臭氧的強(qiáng)氧化性也能起到殺菌消毒的作用。紫外線的殺菌原理與其對核酸、蛋白質(zhì)及酶的作用有關(guān)。當(dāng)水中的細(xì)菌、病毒和藻類等生物受到一定劑量的C波段紫外線照射時，胞體內(nèi)的DNA和RNA結(jié)構(gòu)遭到破壞，基因的復(fù)制、轉(zhuǎn)錄和表達(dá)受阻，蛋白質(zhì)和酶合成受阻，從而使病毒的遺傳物質(zhì)失去活性，細(xì)菌無法進(jìn)行二分裂，藻類因無法進(jìn)行正常的新陳代謝而失去活性，從而達(dá)到消毒和凈化的目的。根據(jù)研究，當(dāng)海洋生物在紫外線的照射下接受的劑量達(dá)到一定值時，微生物不再具有活性。部分微生物在接受紫外線照射之后，在可見光的環(huán)境中，在光復(fù)活酶的作用下，DNA結(jié)構(gòu)會恢復(fù)正常。但是，當(dāng)其接受的總劑量增加到一定值時，將不再具備光復(fù)活性。部分微生物由于缺乏光復(fù)活酶，不具備光復(fù)活性。在所有微生物中，殺滅藍(lán)藻所需劑量最高，達(dá)到300 000 μJ/cm2。

表1 IMO船舶壓載水排放標(biāo)準(zhǔn)

2.2紫外線燈管的選擇

對比各種紫外線光源可知，低壓汞齊燈單根燈管的輸出功率約為300 W，輻射出的紫外線功率可達(dá)到普通低壓汞燈的10倍。汞齊燈獨特的鍍膜涂層保證其不會發(fā)生石英玻璃的透射損失。此外，可在工作16 000 h后依然保持90%以上的紫外線輸出功率。因此,在燈的整個使用壽命周期內(nèi)都會有持續(xù)穩(wěn)定的殺菌效果。汞齊燈對溫度波動不敏感，甚至可在 90 ℃的環(huán)境溫度下使用。[6]由于具有壽命長和功率高的特點，汞齊燈有節(jié)約燈管和系統(tǒng)元件使用數(shù)量、降低能源消耗和維護(hù)時間的巨大潛力，因此可明顯地降低消毒系統(tǒng)的運營成本。

該設(shè)計選擇253.7 nm單波長的低壓汞齊燈。為保證反應(yīng)器的殺菌效果，必須保證壓載水在處理過程中接受紫外線達(dá)到一定的劑量。

2.3紫外線接受總量計算

在計算壓載水接受紫外線的劑量時，將壓載水看作是由微小的點元組成的，通過計算點元從進(jìn)入反應(yīng)器到流出反應(yīng)器這段時間內(nèi)的接受總劑量來確定壓載水接受紫外線的劑量，其接受總劑量Dose的計算式為

(1)

式(1)中：Dose為點元接受總劑量，μJ/cm2；I為反應(yīng)器內(nèi)某單點紫外線照射強(qiáng)度，μW/cm2；t為照射時間，也是點元在反應(yīng)器內(nèi)的停留時間，s。

點元在反應(yīng)器內(nèi)的處理時間t可由路徑長度和點元速度求得，即

t=L/v

(2)

式(2) 中：L為反應(yīng)器縱向長度，m；v為點元縱向分速度，m/s。

由式(1)和式(2)可知，反應(yīng)器設(shè)計的關(guān)鍵在于確定輻照總劑量，輻照總劑量的大小與紫外線燈管的有效輸出功率、輻照強(qiáng)度和總照射時間有關(guān)。燈管的輸出功率為定值。因此，反應(yīng)器設(shè)計的關(guān)鍵在于保證合適的輻照強(qiáng)度、確定反應(yīng)器縱向距離和配置合理的內(nèi)部流場，從而保證壓載水在處理過程中能接受足夠的輻照劑量。

紫外線強(qiáng)度計算采用徑向輻射模型，將紫外線燈看作是一個能量均勻分布的直線光源，并作以下假設(shè)[7]：

1) 紫外線燈輻射能量沿軸向均勻分布。

2) 紫外線垂直于燈管軸線及石英套管外壁，以柱面的形式沿徑向向外傳播能量。

3) 只考慮水中污染物對紫外線的吸收，不計空氣、石英和水界面的折射及水中污染物的散射和反射。

2.4朗伯定律

有色溶液對光的吸收程度與溶液的濃度、液層的厚度及入射光的波長有關(guān)。當(dāng)入射光的波長、溶液的濃度及溫度一定時，溶液的吸光度與液層的厚度成正比，該關(guān)系稱為朗伯定律。[8]

由朗伯定律可知，燈管發(fā)出的紫外線的能量在從管壁傳遞到流域中某一點的過程中，除了會被外層石英玻璃管吸收，還會被管外的水層吸收。因此，在計算水層中的紫外強(qiáng)度分布情況時，要綜合考慮上述2種衰減。《紫外線消毒設(shè)計規(guī)范》[9]中規(guī)定結(jié)晶石英套管在波長為253.7 nm的紫外線透過率≥90%。這里假設(shè)石英玻璃管對紫外線的透過率ηQ=85%。假設(shè)不計紫外線能量在燈管與石英套管間空氣中的衰減，則由朗伯定律可知，距離燈管管壁為d的某點的紫外輻照強(qiáng)度I可表示為

I=I0ηQe-αd

(3)

式(3)中：I為水層厚度為d的點輻照強(qiáng)度，μJ/cm2；I0為石英燈管管壁的紫外線照射強(qiáng)度，μW/cm2；ηQ為石英玻璃管對紫外線的透過率，ηQ=85%；α為水的光吸收系數(shù)，取定值-0.232 6 cm-1。

由采用的輻照模型的3個假設(shè)可知，每只紫外線燈的輻照能量在以燈管為軸線的圓心柱面上均勻分布，且在同一個柱面上任意點的紫外線輻照強(qiáng)度都是相等的。輻照強(qiáng)度的大小隨距離燈管管壁的距離的增大而減小。設(shè)紫外線燈的額定輸入功率為P，波長為253.7 nm 的紫外線的輸出效率為η，燈管的有效長度(弧長)為L，可得到半徑為r的柱面上任意一點的初始紫外線強(qiáng)度I0為

(4)

將式(4)代入式(3)得

(5)

2.5紫外線殺菌動力學(xué)

紫外線殺菌的實質(zhì)是光化學(xué)反應(yīng)，符合一級反應(yīng)動力學(xué)模型。根據(jù)已有的試驗研究結(jié)果[10]，一級反應(yīng)公式對病毒殺滅率的適用性可達(dá)99.999%。

3 紫外線壓載水處理器反應(yīng)器的設(shè)計

3.1設(shè)計要求

船舶在世界各地航行時，會根據(jù)自身需求在不同的水域加載壓載水，因此其壓載水中的細(xì)菌、病毒等生物存在一定的廣泛性，這對處理程序提出了更高的要求。為能殺滅大部分微生物，要求處理過程中的最低總劑量達(dá)到300 000 μJ/cm2。

采用雙入口、單出口設(shè)計。2個入口設(shè)置在反應(yīng)器底部，與桶壁相切，進(jìn)水速度沿著相反的方向；出水口設(shè)置在反應(yīng)器頂部，沿反應(yīng)器中心線向上。用SolidWorks進(jìn)行三維建模，得到反應(yīng)器三維立體圖見圖1。

圖1 反應(yīng)器三維立體圖

3.2燈管布置

反應(yīng)器截面布置圖見圖2，其中小圓為石英玻璃管。每3個石英玻璃管的中心點組成1個等邊三角形，每7塊石英玻璃組成1個正六邊形，形狀類似于蜂窩，故命名為蜂窩排列。[11-14]

圖2 反應(yīng)器截面布置圖

蜂窩排列的創(chuàng)新點：

1) 截面照射強(qiáng)度分布均勻，無照射盲區(qū)。截面可看作是由若干個等邊三角形組成的，最外層等邊三角形的強(qiáng)度分布相同;內(nèi)部所有三角形的強(qiáng)度分布均相同，且略大于外層等邊三角形的強(qiáng)度。

2) 處理能力具備拓展性。良好的截面拓展性使得反應(yīng)器的處理能力具備拓展性，這亦使得設(shè)計大流量反應(yīng)器不再受燈管布置導(dǎo)致截面強(qiáng)度分布不均的限制，僅與燈管管壁之間的間距和反應(yīng)器的高度有關(guān)。

3.3燈管布置

截面內(nèi)任意一點的照射強(qiáng)度是截面所有燈管在該點的照射強(qiáng)度的疊加。有研究表明，當(dāng)水層厚度達(dá)到10 cm時，紫外線穿透率僅為10%。因此,選擇燈管間距的臨界值為10 cm。

為便于計算，對計算區(qū)域進(jìn)行分區(qū)和標(biāo)記(見圖3)。數(shù)字代表區(qū)域，字母代表燈管。將圖3中B,E,I,F所圍區(qū)域作為截面輻照強(qiáng)度計算單元。忽略相隔較遠(yuǎn)的燈管對各區(qū)域的貢獻(xiàn)，將區(qū)域1接受的輻照強(qiáng)度簡化為燈管A，B，D，E，F(xiàn)在區(qū)域1的疊加，區(qū)域2所接受的輻照強(qiáng)度簡化為燈管A，B，C，E，F(xiàn)在區(qū)域2的疊加。由于外層三角形的劑量分布強(qiáng)度略低于內(nèi)部三角形，因此若區(qū)域1的輻照強(qiáng)度滿足設(shè)計要求，則整體設(shè)計方案在截面上的輻照強(qiáng)度滿足設(shè)計要求。

等邊三角形的邊長為L，以燈管D的中心點為原點建立直角坐標(biāo)系。設(shè)點Q為區(qū)域1中的任意一點，則Q點與各燈管之間的距離分別為

(6)

(7)

圖3 計算區(qū)域標(biāo)記圖

(8)

(9)

(10)

(11)

(12)

(13)

各燈管在Q點的劑量為

(14)

代入計算點與計算域各燈管之間的距離得Q點的輻照強(qiáng)度為

(15)

由式(15)可知，Q點的輻照強(qiáng)度是關(guān)于等邊三角形邊長的函數(shù)，當(dāng)水層厚度達(dá)到10 cm時，紫外線穿透率僅為10%。在設(shè)計中，選擇燈管中心距離為11 cm，即每相鄰2根燈管管壁之間的距離為5 cm。

利用MATLAB對紫外線輻照強(qiáng)度進(jìn)行計算，結(jié)果見圖4和圖5。

由圖4可知，在菱形顯示區(qū)域內(nèi)，其紫外線強(qiáng)度所在區(qū)間為150 000～160 000 μW/cm2，平均輻照強(qiáng)度約為155 000 μW/cm2，不存在大的跳躍，紫外線強(qiáng)度分布均勻。整個反應(yīng)器截面由多個相同的三角形區(qū)域組成，因此整個反應(yīng)器截面的輻照強(qiáng)度分布均勻。在以往的設(shè)計方案中，有的采用同心圓的布置方案，有的采用環(huán)形燈管圍繞在處理水流外側(cè)的布置方案，有的通過計算采用不規(guī)則排列的布置方案。對于這些布置方案，紫外線在橫截面的輻照強(qiáng)度都存在輻照暗點，導(dǎo)致截面輻照強(qiáng)度存在大的跳躍，從而導(dǎo)致輻照強(qiáng)度分布不均勻，燈管的截面布置方式不具備可延伸性。與以往的設(shè)計方案相比，該方案在截面輻照強(qiáng)度分布方面有了很大改善。

圖4 計算域輻照強(qiáng)度分布圖

圖5 計算域強(qiáng)度數(shù)值分布

考慮邊緣三角形的強(qiáng)度分布，在計算域內(nèi)(如圖5所示)，大部分區(qū)域的輻照強(qiáng)度分布在148 000～168 000 μW/cm2，最小輻照強(qiáng)度約為137 500 μW/cm2，平均輻照強(qiáng)度約為150 900 μW/cm2，其分布滿足高斯分布。輻照強(qiáng)度遠(yuǎn)大于設(shè)計平均輻照強(qiáng)度100 000 μW/cm2，能滿足設(shè)計要求。

3.4反應(yīng)器內(nèi)部流暢模擬

將模型導(dǎo)入到Gambit內(nèi)進(jìn)行邊界定義及面網(wǎng)格和體網(wǎng)格劃分，將反應(yīng)器模型劃分為入流區(qū)、壁面區(qū)和出流區(qū)等3部分進(jìn)行網(wǎng)格劃分。將劃分好的網(wǎng)格文件導(dǎo)入到Fluent中進(jìn)行流場模擬。流場的模擬以質(zhì)量守恒、動量守恒和能量守恒為基礎(chǔ)，建立連續(xù)方程、動量方程和能量方程對流場內(nèi)各計算點的位置參數(shù)進(jìn)行求解。利用 Fluent 軟件中的離散相[15]湍流模型對裝置內(nèi)的流場進(jìn)行數(shù)值模擬分析，設(shè)置進(jìn)口速度為2.5 m/s。

圖6為反應(yīng)器速度云圖，圖7為質(zhì)量停留時圖，圖8為速度分布，圖9為Z軸方向速度分布。由圖6可知:在整個反應(yīng)器內(nèi)部，速度分布比較均勻，速度較大的位置位于進(jìn)口及出口與反應(yīng)器頂部連接位置；反應(yīng)器內(nèi)的最大流速約為3.8 m/s，出現(xiàn)在出口與反應(yīng)器連接處。由圖8可知：最大速度所占比例非常小，不到1%，不會影響反應(yīng)器的整體處理效果；無明顯的高流速區(qū)域，最小流速為0.000 1 m/s，主要位于反應(yīng)器頂部流場靜區(qū)；在反應(yīng)器內(nèi)部，速度0.0～0.5 m/s所占比例最大(主要分布在反應(yīng)器中上部)，其次為0.5～1.0 m/s(主要分布在反應(yīng)器入口附近)。由圖9可知：沿Z軸方向的分流速主要集中在0.0～1.0 m/s，其中0.0～0.5 m/s分布在反應(yīng)器中上部，速度分布均勻，不存在較大的速度跨越。整體分布比較均勻，這樣的速度分布有利于改善反應(yīng)器的處理效果。結(jié)合圖4和圖7可知：反應(yīng)器截面平均輻照強(qiáng)度約為155 000 μW/cm2，反應(yīng)器平均處理時間約為3.97 s，當(dāng)壓載水流經(jīng)反應(yīng)器時，所接受的紫外線劑量約為615 350 μJ/cm2，能殺死對紫外線劑量要求較高的藍(lán)藻屬生物。

圖6 反應(yīng)器速度云圖圖7 質(zhì)點停留時圖

圖8 速度分布圖9 Z軸方向速度分布

4 結(jié)束語

將紫外線消毒技術(shù)應(yīng)用到壓載水處理領(lǐng)域中，通過對紫外線壓載水處理器的燈管布置進(jìn)行優(yōu)化，并用MATLAB對反應(yīng)器截面的紫外線強(qiáng)度進(jìn)行計算分析，用Fluent軟件對反應(yīng)器內(nèi)部的流場和輻照劑量進(jìn)行仿真分析。

1) MATLAB計算結(jié)果表明：蜂窩排列的截面強(qiáng)度分布均勻，不存在輻照死角，其最小光照強(qiáng)度和平均輻照強(qiáng)度均比以往的設(shè)計方案高。

2) Fluent的模擬結(jié)果表明：反應(yīng)器內(nèi)部流場分布均勻，在處理過程中，壓載水接受的輻射總劑量能滿足殺菌要求。仿真結(jié)果表明：反應(yīng)器頂部存在流場靜區(qū)，在反應(yīng)器的結(jié)構(gòu)上存在繼續(xù)優(yōu)化的空間。反應(yīng)器頂部可設(shè)計為弧形曲面，以消除流場靜區(qū)。若在反應(yīng)器內(nèi)壁增設(shè)螺旋曲線，不僅可改善內(nèi)部流場，還可將螺旋線作為流體導(dǎo)軌，從而使水流在反應(yīng)器內(nèi)形成螺旋上升的流場；水流在流經(jīng)反應(yīng)器內(nèi)部時，所流經(jīng)的路徑得以延長，從而延長接受紫外線照射的時間，使接受的總劑量增加，提高反應(yīng)器的處理效率。由于當(dāng)前三維建模技術(shù)存在一定的局限性，因此在模型內(nèi)部添加螺旋導(dǎo)軌的設(shè)想無法實現(xiàn)，通過仿真的方式驗證模型的有效性存在一定的困難。

[1] 丁尋韜,李智.風(fēng)險預(yù)防原則在海洋生物入侵中的適用研究[J].海峽法學(xué),2014,9(3):72-77.

[2] 張碩慧,劉乒,張爽,等.《船舶壓載水及沉積物控制和管理國際公約》履約面臨的問題及對策[J].水運管理,2009,31(1):29-37.

[3] 黃加亮,李品芳.IMO國際船舶壓載水和沉積物控制與管理公約介紹[J].航海技術(shù),2006(1):46-49.

[4] 劉佳,黃翔峰,沈捷,等.污水紫外消毒微生物光復(fù)活原理及其控制技術(shù)[J].環(huán)境污染與防治,2007,29(11):841-843.

[5] LUIS A D， MELO G D e A， CANESM C A A.UV Dose Investigation for Imersed Lamp Purifier for Electronic Ballast UV Lamps Design[C]∥216th International Power Electronics and Motion Control Conference and Exposition,2014.

[6] 丁有生.高強(qiáng)度低壓汞齊紫外線燈特性探討[J] .燈與照明,2014(4):53-56.

[7] 唐玲.基于CFD的紫外消毒模型的開發(fā)[D] .哈爾濱:哈爾濱工業(yè)大學(xué),2008.

[8] 王建鳴.朗伯—比爾定律的物理意義及其計算方法探討[J].高等函授報,2000,13(3):32-33.

[9] 紫外線消毒設(shè)計規(guī)范:GB 1987—2005[S] .

[10] 譚輝.船舶UV壓載水處理系統(tǒng)反應(yīng)器的設(shè)計研究[D] .大連:大連海事大學(xué),2008.

[11] ZHANG Nanxia, LI Na, XIAO M W.Ballast Water Treatment by Using Porous Ceramic Filter and UV Radition Combination[C].International Conference on Biomedical Engineering and Biotechnology,2012.

[12] ZHANG Manxia, LIU Shengjie,TANG Xiaojia.Evaluation of Micro-Pore Ceramic Filtration and UV Radiation Combination on[C].International Conference on Digital Manufacturing & Automation,2010.

[13] WANG S C, LIU Y H.High Power Factor Electronic Ballast With Intelligent Energy-Saving Control for Ultraviolet Drinking Water Treatment Systems[J].IEEE Transcations on Industrial Electronics,2008,55(1):142-153.

[14] ALONSO J M, LPEZ E, RIBAS J. Design and Implementation of an Electronic Ballast for UV-Based Ozone[C]. Annual Conference of the Industrial Electronics,2002.

[15] 孫國祥,汪小品,丁為民,等.基于CFD離散相模型霧滴沉積特性的模擬分析[J].農(nóng)業(yè)工程學(xué)報,2012,28(6):13-19.

CFDBasedOptimizationDesignandSimulationofUVReactorforBallastWaterTreatment

ZHOUGuoxianga,b,XIAOYingjiea,b,LUYanga,b,WANGYichena,b,LUShu’nana,b

(a. Merchant Marine College; b. Engineering Research Center of Shipping Simulation, Ministry of Education, Shanghai Maritime University, Shanghai 201306, China)

In order to improve the treatment efficiency of the ship’s ballast water with UV radiation method, on the basis of the existing UV radiation method of ballast water treatment, an optimized design of the inlet and outlet of the reactor and the layout of the ultraviolet lamp tube are proposed which improves the flow field in the reactor and the ultraviolet radiation distribution. The effectiveness of these design improvements is verified through simulation with Fluent. A few suggestions about further improvement are put forward.

ultraviolet ballast water treatment; CFD; device optimization; flow field improvement

U664

A