基于FLUENT的水力旋流器入口結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計(jì)流場仿真

龔 俊，葉俊紅，姚明修

(1.中國石化勝利油田分公司海洋采油廠，山東 東營 257237；2.山東省油田采出水處理及環(huán)境污染治理重點(diǎn)實(shí)驗(yàn)室，中石化石油工程設(shè)計(jì)有限公司，山東 東營 257026)

水力旋流器具有分離效率高、占地面積小等優(yōu)點(diǎn)，在國內(nèi)外油田生產(chǎn)系統(tǒng)中得到廣泛使用， 其具體結(jié)構(gòu)形式多種多樣，根據(jù)構(gòu)成特點(diǎn)和工作特點(diǎn)可概括為3 種基本形式，即四段式、三段式和旋轉(zhuǎn)式。水力旋流器的油水分離效率以及油水混合物單位體積處理量所需的能耗是評價(jià)其性能的兩個(gè)重要指標(biāo)。生產(chǎn)實(shí)踐表明水力旋流器的分離效率和能耗常常是一對矛盾，如何妥善處理好這一矛盾，是水力旋流器的技術(shù)關(guān)鍵[1-5]。

水力旋流器入口其形式的優(yōu)劣是決定水力旋流器的流場分布及壓力損失的關(guān)鍵因素之一[6]。水力旋流器的入口形式多種多樣，按其數(shù)量的多少可分為單入口、雙入口及三個(gè)以上的多入口等。按入口流道形式又可分為等截面直線型及變截面( 漸縮截面) 直線型和曲線型[7]。有研究表明，調(diào)整旋流器入口形式不僅大大縮小了液-液水力水力旋流器的空間，并且能顯著提高分離效率[8-10]。

本文設(shè)計(jì)了一種雙入口形式的水力旋流器，通過CFD軟件仿真優(yōu)化了設(shè)計(jì)參數(shù)，為水力旋流器入口結(jié)構(gòu)參數(shù)的合理選擇提供新思路。

1 CFD模型建立及網(wǎng)格劃分

水力旋流器選用四段式，分為旋流腔、大錐段、小錐段、平尾段四個(gè)區(qū)域(見圖1)。其中，旋流腔和大錐段統(tǒng)稱為旋流體，也是水力旋流器的主體部分，通常是由上部的圓柱段(旋流腔)與下部的圓錐段(大錐段)組成。液體從切向入口進(jìn)入旋流腔內(nèi)產(chǎn)生高速旋轉(zhuǎn)的液流。旋流腔的直徑是水力水力旋流器的主直徑，其大小不但決定了水力水力旋流器的處理能力，而且其和旋流腔、大錐段長度的比例關(guān)系也是確定其他參數(shù)的重要依據(jù)，在很大程度上影響著水力水力旋流器油水分離的能力。

圖1 四段式水力旋流器結(jié)構(gòu)示意圖

1.1 邊界條件

(1)入口邊界條件

入口溫度：45～60℃；入口流量：1～2m3/h；入口壓力：0.5～0.7MPa。

(2)出口邊界條件出口包含收油口和出水口，均按流動(dòng)充分發(fā)展條件處理。

(3)固體壁面邊界條件

固壁邊界按照無滑移邊界條件處理，默認(rèn)壁面粗糙度為0.5。本次計(jì)算按照壁面不可滲透，不存在滑移速度考慮，并使用標(biāo)準(zhǔn)壁面函數(shù)法確定固壁附近流動(dòng)。

1.2 多相流模型的選擇

多相流模型選擇“Euler-Mixture”模型進(jìn)行計(jì)算。

1.3 模型建立與網(wǎng)格劃分

在眾多入口結(jié)構(gòu)中，阿基米德螺線形入口較其他結(jié)構(gòu)入口在零件外徑相同的前提下，流道更長，導(dǎo)流能力更強(qiáng)，更有利于混合液迅速形成較為穩(wěn)定的旋流，使旋流腔內(nèi)流場分布更加有序合理，降低了液滴剪切破碎的可能性[6]。根據(jù)阿基米德螺線標(biāo)準(zhǔn)方程，嚴(yán)格按照計(jì)算步驟完成水力旋流器入口結(jié)構(gòu)設(shè)計(jì)，并結(jié)合加工精度要求和壁厚等參數(shù)得到最優(yōu)化的入口結(jié)構(gòu)形式，見圖2。圖2中標(biāo)號1 的曲線為阿基米德螺線，標(biāo)號2 是導(dǎo)流線，標(biāo)號3 的圓為旋流腔內(nèi)孔圓，點(diǎn)A 是阿基米德螺線與旋流腔內(nèi)孔圓相交點(diǎn)。

圖2 阿基米德螺線進(jìn)口結(jié)構(gòu)示意圖

利用了ICEM-CFD軟件對水力旋流器的物理模型劃分了全六面體網(wǎng)格，網(wǎng)格數(shù)量控制在100萬左右。在進(jìn)水口交界面及細(xì)微區(qū)域都進(jìn)行了一定的加密，確保得到更加精確的計(jì)算結(jié)果，見圖3。

圖3 水力旋流器模型

圖4 模型進(jìn)水口放大圖

在網(wǎng)格劃分過程中，對網(wǎng)格的主要尺寸都進(jìn)行了最大與最小值網(wǎng)格質(zhì)量評測工作，網(wǎng)格質(zhì)量在0.4以上，滿足計(jì)算要求。

圖5 網(wǎng)格劃分模型

2 仿真結(jié)果及分析

通過仿真計(jì)算得到了水力旋流器分離效率、內(nèi)部的壓力分布、內(nèi)部的氣體流線分布以及各截面面積等相應(yīng)的數(shù)值。

2.1 流線分布

圖6 水力旋流器內(nèi)部的流線分布

由圖6水力旋流器內(nèi)部截面速度分布圖可以發(fā)現(xiàn)，進(jìn)入腔體內(nèi)的液體大多數(shù)都是基于流入的速度沿著壁面在高速轉(zhuǎn)動(dòng)，受重力以及壓力的影響，水流向下運(yùn)動(dòng)，在中心筒附近形成一個(gè)流動(dòng)較密集的區(qū)域。而中心區(qū)域液體受壓力影響向頂部出口流動(dòng)。

2.2 速度分布

圖7 水力旋流器內(nèi)部的速度分布

由圖7截面的速度分布矢量圖可以發(fā)現(xiàn)明顯顯示環(huán)狀的速度分布，越向內(nèi)部速度逐漸變低，中心點(diǎn)的速度則變大，其方向也由環(huán)狀變成了向垂直方向變化。

2.3 速度矢量分布

圖8 旋流器中間截面的速度分布矢量圖

經(jīng)過大量的模擬運(yùn)算得到水力旋流器各個(gè)結(jié)構(gòu)尺寸的優(yōu)化調(diào)節(jié)范圍，最終獲得圖8較理想的的內(nèi)部速度矢量分布結(jié)果。

2.4 壓力分布

圖9為沿著頂部橫截面上的壓力云圖，可以發(fā)現(xiàn)，水力旋流器中心壓力較低，兩邊的壓力較高，壓力成環(huán)狀向外分布，壓力分布由內(nèi)向外均勻分布。

圖10為縱截面壓力云圖，可以發(fā)現(xiàn)，水力旋流器自入口至出水口壓力逐漸降低。

圖9 中間截面的壓力分布云圖

圖10 縱截面上的壓力云圖分布

2.5 油水分離效果分析

圖11水力旋流器縱截面含油體積分布圖和圖12水力旋流器橫截面含油體積分布圖結(jié)果顯示，液體中的油相集中的水力旋流器的中心區(qū)域，體積分?jǐn)?shù)相對較高，油相部分通過中心出油管排出水力旋流器，起到了較好的油水分離效果。

圖11 水力旋流器縱截面含油體積分?jǐn)?shù)分布圖

圖12 水力旋流器頂部出口含油體積分?jǐn)?shù)分布圖

3 結(jié)論

入口結(jié)構(gòu)對水力旋流器內(nèi)部流場分布影響較大，設(shè)計(jì)了阿基米德螺旋線入口形式的水力旋流器，仿真內(nèi)部速度場具有良好的對稱性，內(nèi)部流動(dòng)穩(wěn)定性好、渦流區(qū)域少、離心力場強(qiáng)，具有較高的分離效率。