基于本征正交分解法的液環(huán)泵氣液兩相流場重構(gòu)

張人會 吳 昊 楊軍虎 李仁年

(1.蘭州理工大學能源與動力工程學院， 蘭州 730050； 2.蘭州理工大學甘肅省流體機械及系統(tǒng)重點實驗室， 蘭州 730050)

基于本征正交分解法的液環(huán)泵氣液兩相流場重構(gòu)

張人會1,2吳 昊1楊軍虎1,2李仁年1,2

(1.蘭州理工大學能源與動力工程學院， 蘭州 730050； 2.蘭州理工大學甘肅省流體機械及系統(tǒng)重點實驗室， 蘭州 730050)

針對液環(huán)泵內(nèi)復(fù)雜的氣液兩相流動、計算量大且優(yōu)化設(shè)計難以進行的問題，提出了采用本征正交分解法進行液環(huán)泵內(nèi)流場的重構(gòu)分析。采用泰勒多項式對二維葉片型線進行參數(shù)化控制，通過對各控制參數(shù)進行適量的擾動得到葉片型線樣本集。采用VOF模型進行液環(huán)泵內(nèi)氣液兩相流場的數(shù)值模擬，由葉片型線參數(shù)及葉輪內(nèi)流場數(shù)據(jù)構(gòu)建樣本的快照集，依照幾何相似及網(wǎng)格變形方法插值得到各相似點的流場參數(shù)，依據(jù)本征正交分解(POD)法將快照集分解為正交基的線性組合。由最小二乘法擬合目標葉型所對應(yīng)的正交基系數(shù)，實現(xiàn)了對目標葉片流場的重構(gòu)。采用POD方法對2BE-203型液環(huán)泵內(nèi)單個葉輪流道的氣液兩相流流場進行了重構(gòu)，精確地重構(gòu)了單個葉輪內(nèi)流場結(jié)構(gòu)的各個特征，除在氣液交界面附近有一定的誤差，整體預(yù)測結(jié)果具有較高的精度，大大減少了流場預(yù)估的計算量。

本征正交分解; 氣液兩相流; 液環(huán)泵; 網(wǎng)格變形

引言

液環(huán)泵是一種用來抽送氣體的流體機械, 被廣泛應(yīng)用于低真空或低壓縮比場合，其內(nèi)部為具有不穩(wěn)定自由分界面的氣液兩相流動，效率較低。由于復(fù)雜的內(nèi)部兩相流動約束導致其水力性能的優(yōu)化設(shè)計進展緩慢。RAIZMAN等[1]采用壓力探針測量了液環(huán)泵內(nèi)液環(huán)的壓力分布，并對液環(huán)的速度場進行了分析。KAKUDA 等[2]運用移動離子半隱式方法(MPS)對液環(huán)泵內(nèi)部粘性流動進行數(shù)值計算，分析了不同轉(zhuǎn)速下液環(huán)的形狀、速度矢量等，并將計算結(jié)果與實驗結(jié)果進行了對比。TETERYYUKOV等[3]提出了采用三段圓弧方法設(shè)計非圓柱形泵殼體來優(yōu)化液環(huán)泵的性能。文獻[4-5]運用CFD 軟件對液環(huán)泵內(nèi)部流動進行數(shù)值計算，分析了泵內(nèi)兩相流場的分布規(guī)律。張人會等[6]采用CFD方法對液環(huán)泵內(nèi)氣液兩相流及其自由界面的復(fù)雜運動規(guī)律進行了研究。液環(huán)泵內(nèi)復(fù)雜的氣液兩相流動導致其內(nèi)流動數(shù)值模擬的計算量大大增大。

流體機械水力優(yōu)化設(shè)計的流場計算量隨著其設(shè)計變量的增加呈幾何級數(shù)增加[7-9]，因此其水力性能優(yōu)化問題的重點在于減少流場預(yù)估的計算量。為了減少流場的計算量，ANTONY[10]提出伴隨方法(Adjoint method)，并將該方法應(yīng)用于航空翼型的氣動優(yōu)化設(shè)計，該方法最大的困難在于要求解復(fù)雜的伴隨系統(tǒng)。MOHAMMADI等[11]提出不完全敏感性方法來減少流場約束下幾何邊界的優(yōu)化問題，但該方法有苛刻的適用條件。本征正交分解(Proper orthogonal decomposition，POD)法最早由HOTELLING[12]提出，是一種功能強大的數(shù)據(jù)特征提取方法，被廣泛應(yīng)用于低維湍流結(jié)構(gòu)的描述[13]、結(jié)構(gòu)振動分析[14]、圖像處理[15]、PIV技術(shù)[16]、以及翼型流場的重構(gòu)分析及其氣動優(yōu)化[17-19]。本文提出基于本征正交分解(POD)方法進行液環(huán)泵內(nèi)氣液兩相流場的重構(gòu)，避免其優(yōu)化過程敏感性分析中的多次流場計算，減少優(yōu)化設(shè)計的計算量。

1 液環(huán)泵流場分析

以2BE-203型液環(huán)泵為研究對象，對其內(nèi)部復(fù)雜的氣液兩相流場進行重構(gòu)分析，其結(jié)構(gòu)如圖1所示。

圖1 液環(huán)泵結(jié)構(gòu)圖Fig.1 Structures of liquid-ring pump

其流動求解區(qū)域包括：葉輪、腔體、進氣段、出氣段和補液管?？紤]其為雙吸式泵，兩側(cè)流道完全對稱，在計算過程中只模擬其單側(cè)流動，對整個流動區(qū)域使用六面體結(jié)構(gòu)化網(wǎng)格進行劃分，網(wǎng)格總數(shù)1 091 166，其計算域網(wǎng)格如圖2所示。

圖2 液環(huán)泵求解區(qū)域網(wǎng)格Fig.2 Mesh of computational domain for liquid-ring pump

其中任一葉輪流道有39×30個網(wǎng)格單元，如圖3所示。

圖3 葉輪流道網(wǎng)格Fig.3 Blades mesh domain

液環(huán)泵額定工況為：進氣壓力p=40 kPa，進水流量Q=0.43 m3/h，轉(zhuǎn)速n=1 170 r/min。液環(huán)泵內(nèi)流動為具有非穩(wěn)態(tài)氣液自由分界面的氣液兩相流，氣液分解面的形狀對泵的性能有重要影響，為精確捕捉界面流動，采用VOF氣液兩相流模型[20]，采用RNGk-ε兩方程湍流模型、理想氣體定律、滑移網(wǎng)格、PISO算法，時間步長Δt=3×10-5s，由32核64 G內(nèi)存小型工作站120 h計算后，得到較為穩(wěn)定的壓力、相態(tài)、速度場分布，如圖4所示。

由圖4可以看出壓力由輪轂到殼體內(nèi)壁、由吸氣區(qū)到排氣區(qū)逐漸增大；在泵殼體內(nèi)壁形成近似等厚度的液環(huán)，氣液交界面呈不穩(wěn)定的鋸齒形；流道內(nèi)速度基本上從輪轂到殼體內(nèi)壁逐漸增大，但從吸氣區(qū)到壓氣區(qū)逐漸減小。

圖4 流道截面的流場分布Fig.4 Flow distributions of channel section

液環(huán)泵性能曲線的模擬結(jié)果與實驗結(jié)果對比如圖5所示。實驗結(jié)果由合作生產(chǎn)廠家的液環(huán)泵實驗系統(tǒng)測試得到。

圖5 模擬結(jié)果與實驗結(jié)果對比Fig.5 Comparison of hydraulic performance between CFD and experiment

從圖5可看到，隨著進氣口真空度的增大，其質(zhì)量流量逐漸下降，當進氣口吸氣壓力大于40 kPa時，其模擬值與實驗值吻合較好，當進氣口壓力低于40 kPa時，模擬值與實驗值之間有一定的偏差，隨著進氣口真空度的提高，液環(huán)泵內(nèi)氣液交界面變得越來越不穩(wěn)定，VOF模型對自由界面的捕捉精度下降。本研究對液環(huán)泵流場重構(gòu)分析中選定40 kPa工況作為預(yù)測工況。

2 Gappy POD本征正交分解方法

POD方法利用降階技術(shù)對復(fù)雜系統(tǒng)進行特征分析，對樣本集構(gòu)成的矢量集進行本征正交分解，樣本集矢量可表示為正交基的線性組合。POD方法常被用于缺失數(shù)據(jù)的填補。

樣本快照集U為

U=[U1,U2]

(1)

在矢量集U1中，所有的元素都是已知的，而在矢量集U2中，部分元素是缺失的，所缺失的元素需要補充完整。矢量集U1為

(2)

式中，Φ是POD基礎(chǔ)，可通過快照集U的奇異值分解計算得到，即

Φ=[Φ1,Φ2]

(3)

式中，Φ1、Φ2對應(yīng)于快照集U1和不完整矢量U2的POD基礎(chǔ)。

矢量U2可表示為

(4)

系數(shù)βj可通過最小二乘逼近法得到，即

Mβ=f

(5)

其中Mij=(Φi,Φj)，fi=(U1,Φi)，通過求解式(5)得到系數(shù)βj，所以不完整矢量U2可通過式(4)進行補充。

3 葉片型線控制方法和網(wǎng)格變形

3.1 葉片型線參數(shù)控制

由于液環(huán)泵葉片都屬于二維圓柱形葉片，圓柱葉片型線上的點可以表示為：f(r,θ)=0。通過采用泰勒公式來表示葉片型線，任意葉片型線均可以表示為其初始型線基礎(chǔ)上的泰勒展開，即

θ(r)=θ0+a1r+a2r2+…+anrn

(6)

在初始葉型的基礎(chǔ)上對各控制參數(shù)進行擾動可得流場重構(gòu)所需的樣本集，本算例中以圖6所示目標葉型附近的8個葉型作為初始樣本集。

圖6 樣本葉片型線Fig.6 Blade shape curves of samples

3.2 相似網(wǎng)格及變形方法

POD方法是根據(jù)各樣本葉型的內(nèi)部流場對目標葉型內(nèi)流場進行重構(gòu)分析，由于各葉輪內(nèi)流道幾何相似，進行流場的重構(gòu)時需要根據(jù)相似坐標進行流場樣本集的構(gòu)建，并進行流場的預(yù)測。因此對各葉輪樣本的單個流道均進行39×30的結(jié)構(gòu)化網(wǎng)格劃分。而對部分由于計算過程收斂性要求對葉輪結(jié)構(gòu)化網(wǎng)格進行細化的葉輪，在CFD流場分析完成后，采用網(wǎng)格變形技術(shù)重構(gòu)其39×30的結(jié)構(gòu)化網(wǎng)格，然后按照相似的網(wǎng)格節(jié)點進行流場的插值，從而得到各相似坐標點對應(yīng)的樣本矢量集。

定義沿著圓周角θ的無量綱坐標εi,j為

εi,j=(θi,j-θi,1)/(θi,30-θi,1)

(7)

式中θi,j為(i,j)節(jié)點所對應(yīng)的圓周角坐標，類似地可定義r方向的無量綱坐標ηi,j。對于根據(jù)目標葉輪所生成的39×30的結(jié)構(gòu)化網(wǎng)格可重構(gòu)其他樣本葉輪的結(jié)構(gòu)化網(wǎng)格，如圖7所示。

圖7 網(wǎng)格變形對比Fig.7 Comparison of mesh deformation

4 葉輪內(nèi)流場的重構(gòu)及其誤差分析

由此可得其樣本矢量集為

(8)

式中vo,i——目標葉型(αo,i)所對應(yīng)的流場參數(shù)

z——樣本數(shù)

通過式(2)對剔除目標樣本矢量的矢量集U進行本征正交分解；通過式(5)進行二乘擬合，最后可由式(4)對式(8)中的目標葉型流場數(shù)據(jù)vo,i進行重構(gòu)。預(yù)測流場與CFD分析流場對比及其誤差分別如圖8～13所示。

圖8 CFD與POD重構(gòu)的壓力場分布Fig.8 Pressure distributions of CFD and POD

圖9 壓力場誤差分布Fig.9 Distribution of pressure errors

圖8為單個葉輪內(nèi)流道的壓力場預(yù)測結(jié)果與CFD對比，由圖可以看出POD預(yù)測出的壓力場與其CFD分析結(jié)果基本一致，從葉輪的輪轂到輪緣壓力逐漸增大，在葉片工作面出口附近壓力最高。圖9為壓力預(yù)測誤差(pCFD-pPOD)的分布情況，可以看出除在氣液交界面的葉片工作面附近出現(xiàn)一個較高的局部誤差(約10 000 Pa左右)外，整個葉輪內(nèi)流道的預(yù)測誤差基本在2 000～4 000 Pa左右,該值相對于葉輪出口平均壓力的相對誤差為1.7%～3.4%。

圖10 CFD與POD重構(gòu)的體積分數(shù)分布Fig.10 Volume fraction distributions of POD and CFD

圖11 體積分數(shù)誤差分布Fig.11 Distribution of volume fraction errors

圖10為單個葉輪內(nèi)流道的氣液相態(tài)的預(yù)測結(jié)果與CFD結(jié)果的對比(α=1為液相，α=0為氣相)，可以看出POD方法能較為精確地重構(gòu)其氣液分解面的形狀，而且還重構(gòu)了葉輪中心的一個微小氣泡。從圖11的誤差分布可以看出，除在氣液交界面附近的相態(tài)有一定誤差外，其余整個葉輪流道的相態(tài)體積分數(shù)誤差平均約為0.013。

圖12為單個葉輪內(nèi)流道相對速度場的重構(gòu)結(jié)果與CFD結(jié)果的對比，可以看出POD方法精確地重構(gòu)出了原葉輪內(nèi)相對速度場的分布特征，在葉片背面靠近出口處、氣液交界面靠近葉片的工作面附近等幾個局部的高速區(qū)域都被POD方法重構(gòu)出來。圖13為重構(gòu)速度場與CFD速度場間的誤差，重構(gòu)速度場誤差的分布與壓力場重構(gòu)誤差分布相似，除在氣液交界面及葉片出口的工作面附近存在一定的誤差外，其余葉輪內(nèi)大部分區(qū)域的相對速度誤差在0.2 m/s左右。

圖12 CFD與POD重構(gòu)的速度場分布Fig.12 Velocity distributions of POD and CFD

圖13 速度場誤差分布Fig.13 Distribution of velocity errors

所提出的POD方法可以精確地實現(xiàn)對液環(huán)泵內(nèi)單個葉輪流道的復(fù)雜氣液兩相流場的重構(gòu)，對其他葉輪流道流場及殼體流場的重構(gòu)方法完全一樣。該方法可以作為代理模型用于對液環(huán)泵優(yōu)化過程中流場的預(yù)估，減少對復(fù)雜流場的CFD模擬次數(shù)，大大提高了復(fù)雜流動優(yōu)化的速度。

5 結(jié)論

(1)提出采用本征正交分解法，對液環(huán)泵內(nèi)的多相流場進行重構(gòu)分析。對由葉型控制參數(shù)及流場參數(shù)構(gòu)成的樣本矢量集分解為正交基的線性表示，由最小二乘法擬合目標葉型在該正交基上的各個系數(shù)，可以得到目標葉型所對應(yīng)的各流場參數(shù)。

(2)從對2BE-203型液環(huán)泵單個葉輪流道內(nèi)的氣液兩相流的重構(gòu)結(jié)果看，POD方法能精確地重構(gòu)出其壓力場、速度場及相態(tài)場分布的結(jié)構(gòu)特征，除在氣液交界面處有一定誤差外，整個葉輪內(nèi)壓力的平均相對誤差為1.7%～3.4%，相態(tài)體積分數(shù)平均誤差約為0.013，相對速度誤差平均值為0.2 m/s。

(3)算例計算結(jié)果表明所提出的POD方法能夠精確地對液環(huán)泵內(nèi)多相流場進行精確的重構(gòu)，可以作為代理模型用于流體機械優(yōu)化設(shè)計中對流場的預(yù)估，大大減少了優(yōu)化設(shè)計過程中流場的計算量。

1 RAIZMAN I A, MATS é B. Experimental investigation of the velocity field in a liquid ring vacuum-pump[J]. Chemical and Petroleum Engineering, 1972, 8(2): 134-137.

2 KAKUDA K, USHIYAMA Y, OBARA S, et al. Flow simulations in a liquid ring pump using a particle method[J]. CMES, 2010, 66(3): 215-226.

3 TETERYUKOV V I. Effect of body shape on performance of water-ring vacuum pumps[J]. Chemical and Petroleum Engineering, 1966, 2(8): 511-513.

4 黃苗苗, 李國君, 匡曉峰. 水環(huán)真空泵內(nèi)部氣液兩相流動的數(shù)值分析[J].船舶力學,2011,15(7):722-729. HUANG Miaomiao, LI Guojun, KUANG Xiaofeng. Numerical simulation of the gas-liquid tow-phase flow in water-ring vacuum pump[J]. Journal of Ship Mechanics,2011,15(7): 722-729. (in Chinese)

5 黃思, 阮志勇, 鄧慶健,等. 液環(huán)真空泵內(nèi)氣液兩相流動的數(shù)值分析[J].真空,2009,46(2):49-52. HUANG Si,RUAN Zhiyong,DENG Qingjian, et al. Numerical analysis of gas-liquid tow-phase flow in liquid-ring vacuum pump[J]. Vacuum, 2009, 46(2): 49-52. (in Chinese)

6 張人會，郭廣強，楊軍虎,等. 液環(huán)泵內(nèi)部氣液兩相流動及其性能分析[J/OL].農(nóng)業(yè)機械學報, 2014, 45(12):99-103. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20141216&flag=1.DOI:10.6041/j.issn.1000-1298.2014.12.016. ZHANG Renhui, GUO Guangqiang, YANG Junhu, et al. Investigation on the inner gas-liquid flow in liquid-ring pump and it’s performance[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(12):99-103. (in Chinese)

7 MOHAMMADI B, PIRONNEAU O. Applied shape optimization for fluids[M]. Oxford: Oxford University Press, 2010.

8 張人會，鄭凱，楊軍虎,等. 基于不完全敏感性方法的低比轉(zhuǎn)速離心葉輪優(yōu)化研究[J].機械工程學報, 2014, 50(4):162-166. ZHANG Renhui, ZHENG Kai, YANG Junhu, et al. The optimization of low specific speed centrifugal impeller based on incomplete sensitivities[J]. Journal of Mechanical Engineering, 2014, 50(4):162-166. (in Chinese)

9 IULIANO E， PEREZ E A. Application of surrogate-based global optimization to aerodynamic design[M]. Berlin: Springer International Publishing, 2016.

10 ANTONY J. Aerodynamic design via control theory[J]. Journal of Scientific Computing, 1988(3): 402-407.

11 MOHAMMADI B, MOLHO J I, SANTIAGO J G. Incomplete sensitives for design of minimal dispersion fluidic channels[J]. Computer Methods in Applied Mechanics and Engineering, 2003, 192: 4131-4145

12 HOTELLING H. Analysis of a complex of statistical variables into principal components[J]. Journal of Education Psychology, 1933, 24: 417-441, 498-520.

13 HOLMES P, LUMLEY J L, BERKOOZ G. Turbulence, coherent structures, dynamical systems and symmetry[M]. Cambridge: Cambridge University Press, 1998.

14 FEENY B F, KAPPAGANTU R. On the physical interpretation of proper orthogonal modes in vibrations[J]. Journal of Sound and Vibration, 1998, 211(4): 607-616.

15 SIROVICH L, KIRBY M. Low-dimensional procedure for the characterization of human faces[J]. Journal of the Optical Society of America, 1987, 4(3): 519-524.

16 RABEN S G, CHARONKO J J, VlACHOS P P. Adaptive gappy proper orthogonal decomposition for particle image velocimetry data reconstruction[J]. Measurement Science and Technology, 2012, 23(2): 25303-25318.

17 BUI-THANH T, DAMODARAN M, WILLCOX K E. Aerodynamic data reconstruction and inverse design using proper orthogonal decomposition[J]. AIAA Journal, 2004, 42(8): 1505-1516.

18 OYAMA A, NONOMURA T, FUJII K. Data mining of Pareto-optimal transonic airfoil shapes using proper orthogonal decomposition[J]. Journal of Aircraft, 2010, 47(5): 1756-1762.

19 白俊強, 邱亞松, 華俊. 改進型Gappy POD翼型反問題方法[J].航空學報,2013,34(4):762-771. BAI Junqiang, QIU Yasong,HUA Jun. Inverse design of airfoil based on Gappy POD method[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(4): 762-771. (in Chinese)

20 HIRT C W, NICHOLS B D. Volume of fluid(VOF) method for the dynamics of free boundaries[J]. Journal of Computational Physics,1981,39(1): 201-225.

Reconstruction for Gas-Liquid Flow of Liquid-ring Pump Based on Proper Orthogonal Decomposition

ZHANG Renhui1,2WU Hao1YANG Junhu1,2LI Rennian1,2

(1.CollegeofEnergyandPowerEngineering,LanzhouUniversityofTechnology,Lanzhou730050,China2.KeyLaboratoryofFluidMachineryandSystems,GansuProvince,LanzhouUniversityofTechnology,Lanzhou730050,China)

According to the complicated gas-liquid flow in liquid-ring pump, the large amount of calculation, and the difficulty for its optimization design, the flow field reconstitution for the gas-liquid flow of liquid-ring pump based on proper orthogonal decomposition were proposed. The blade was parameterized by Taylor polynomial, and the experiment samples can be designed by introducing small perturbation of the control parameter. The transient gas-liquid flow in liquid-ring pump was simulated by using the VOF model. The snapshot set consisted of the control parameter for the blade shape and the flow field data. According to the geometric similarity of the impeller flow passage and the mesh deformation technology, the flow field data of the similar position of each point were interpolated. The snapshot set can be decomposed as linear combination of orthogonal basis by using proper orthogonal decomposition. The coefficients of the objective blade orthogonal basis were fitted by the least square method. The flow field data of the objective blade were reconstructed. In the calculation case for the type of 2BE-203 liquid-ring pump, the gas-liquid flows of single blade passage were reconstructed, and almost all the flow structures were accurately predicted. The prediction had high accuracy except near the gas-liquid interface. The calculation amount for the multiphase flow field was greatly reduced.

proper orthogonal decomposition; gas-liquid flow; liquid-ring pump; mesh deformation

10.6041/j.issn.1000-1298.2017.06.050

2016-09-22

2016-11-10

國家重點研發(fā)計劃項目(2016YFB0200901)和國家自然科學基金項目(51469014)

張人會(1977—)，男，副教授，博士生導師，主要從事流體機械內(nèi)部流動及性能優(yōu)化研究，E-mail： zhangrh@lut.cn

TH311

A

1000-1298(2017)06-0381-06