張 棣，董祥偉，李增亮，鹿勝楠

(1.海洋石油工程股份有限公司，天津 300451；2.中國石油大學(華東)，山東 青島 266580)

?

射流泵內(nèi)流場數(shù)值計算的結構化網(wǎng)格劃分方法

張 棣1，董祥偉2，李增亮2，鹿勝楠1

(1.海洋石油工程股份有限公司，天津 300451；2.中國石油大學(華東)，山東 青島 266580)

射流泵；數(shù)值模擬；結構化網(wǎng)格；疏密控制

射流泵是一種流體輸送機械及混合反應設備，其特點是本身沒有運動部件，結構簡單，且工作可靠，廣泛應用于石油開采、農(nóng)牧漁業(yè)、水利電力、交通運輸和環(huán)境保護等多個領域[1]。射流泵內(nèi)部流體流動屬于高雷諾數(shù)的強剪切湍流流動，流動不具備相似性，使得長期以來人們對內(nèi)部流動機理認識不充分，造成研究手段單一，過分依賴于工程經(jīng)驗或試驗結果。數(shù)值模擬是近幾十年逐漸發(fā)展成熟的一種研究方法，可以對復雜流動問題進行可視化研究，并且具有很高的工程應用價值。因此，學者們將數(shù)值模擬技術應用到射流泵的研究當中，研究內(nèi)容涉及射流泵內(nèi)流場特性分析、氣蝕性能評價、外特性預測、結構參數(shù)優(yōu)化等多個方面，數(shù)值模擬成為研究射流泵的重要手段[2-4]。

數(shù)值網(wǎng)格劃分屬于數(shù)值模擬技術的前處理范疇，在數(shù)值計算理論日趨完善的今天仍然被當做最關鍵的環(huán)節(jié)之一。高質(zhì)量的網(wǎng)格對于保證計算精度是至關重要的，Thompson教授于20世紀70年代發(fā)展了網(wǎng)格貼體坐標生成技術，對提高復雜幾何形狀計算域網(wǎng)格的正交性有很好的普適性[5]；自適應網(wǎng)格可以通過與計算數(shù)據(jù)的交互來實時調(diào)整網(wǎng)格，能夠解決具有復雜幾何邊界以及流場變量變化劇烈?guī)淼木W(wǎng)格劃分尺度不一的問題。但這些方法都是以精密的算法和犧牲計算時間為前提，用于流場特征明顯的射流泵流場計算，并不是有效的解決方法。因此，研究高質(zhì)量的、適合于射流泵流場計算的數(shù)值計算網(wǎng)格對于射流泵的數(shù)值模擬研究有重要意義。

雖然國內(nèi)外涉及到射流泵數(shù)值模擬的研究很多，但到目前為止還沒有一篇論文專門來探討射流泵流場數(shù)值計算網(wǎng)格的劃分方法。本文以中心式射流泵為研究對象，針對射流泵幾何結構及流場特征，較為詳細地介紹了適合于射流泵流場的結構化網(wǎng)格劃分方法，提出了射流剪切層網(wǎng)格參數(shù)的確定方法，并通過數(shù)值計算分析剪切層網(wǎng)格參數(shù)對數(shù)值計算結果的影響規(guī)律，通過某型射流泵的實際計算，驗證了網(wǎng)格的精確性及適用性。

1 網(wǎng)格劃分研究模型

網(wǎng)格劃分模型基于工程中常用的中心式射流泵，如圖1所示，除吸入管外整體呈軸對稱結構，噴嘴位于模型中軸線上，主要特征結構參數(shù)包括噴嘴直徑D0、面積比S、喉嘴距l(xiāng)1、喉管長度l2等。將模型計算域描述為一組柱錐結構體通過聯(lián)合、相交、相減等布爾運算方式形成的組合結構體，根據(jù)其幾何特征，研究其結構化網(wǎng)格劃分方法。

圖1 射流泵幾何結構示意

2 網(wǎng)格劃分方法

2.1 掃掠法及計算域分塊

射流泵的幾何結構呈現(xiàn)掃掠特征，適合于使用掃掠法進行網(wǎng)格劃分。掃掠法是一種2.5D的結構化網(wǎng)格生成方法，可操作性強，生成的網(wǎng)格正交性好，并且便于邊界層網(wǎng)格的加密。

按照掃掠方式的不同，將計算域分為2塊，如圖2所示，分別是內(nèi)塊(G)和外塊(S)，其中G塊掃掠方式分別為拉伸掃掠，S塊的掃掠方式為旋轉掃掠。2塊的網(wǎng)格在分界面處共享網(wǎng)格節(jié)點。

圖2 掃掠法及計算域分塊

掃掠法的3要素包括：源面、目標面和掃掠線。首先對源面劃分面網(wǎng)格，通過將面網(wǎng)格沿掃掠線的逐層掃掠生成體網(wǎng)格。

2.2 面網(wǎng)格的劃分及網(wǎng)格加密技術

面網(wǎng)格的生成方法有映射法(Mapping)和平鋪法(Paving)等。映射法生成的四邊形網(wǎng)格正交性好，便于邊界網(wǎng)格的加密。平鋪法對復雜邊界的適應性更好[6]。

以平行于軸線的縱剖面作為S塊源面，通過面分塊減弱幾何形狀不規(guī)則對網(wǎng)格劃分帶來的影響，從而保證網(wǎng)格與邊界的正交性，采用映射法劃分網(wǎng)格；以垂直于軸向的橫斷面作為G塊源面，采用映射法與平鋪法結合的方式劃分面網(wǎng)格，如圖4的A向視圖所示，為便于邊界處網(wǎng)格加密，外圍環(huán)面使用映射法，內(nèi)部圓面使用平鋪法。

圖3展示了2種不同幾何特征的映射網(wǎng)格加密方法，基本步驟相似：利用坐標變換原理，對網(wǎng)格進行拉伸或壓縮，沿網(wǎng)格線控制網(wǎng)格密度。體網(wǎng)格的加密方案為：首先對源面面網(wǎng)格進行加密處理，然后通過掃掠法的實施，將面網(wǎng)格的加密方案傳遞到體網(wǎng)格中，充分利用掃掠法的特點——沿掃掠路徑保持面網(wǎng)格的網(wǎng)格拓撲特性。

結合面網(wǎng)格的劃分方法及映射網(wǎng)格加密技術，進行網(wǎng)格疏密控制方案研究。

圖3 映射網(wǎng)格加密技術

2.3 網(wǎng)格疏密控制方案研究及實施

網(wǎng)格疏密控制的目的是控制網(wǎng)格的數(shù)量及尺度，同時使計算域中流體變量變化劇烈的位置有足夠的分辨率[7]。射流泵流場屬于有限空間射流，有限空間射流流場內(nèi)部存在射流邊界層與固壁邊界層兩個流速梯度較大的區(qū)域，為較好地模擬該流動，應該在剪切層內(nèi)沿徑向?qū)W(wǎng)格進行加密[8-9]；相反，在非剪切層區(qū)，可使用相對較疏的網(wǎng)格密度。除此之外，由于勢流核區(qū)很短，且喉管區(qū)前部軸向流速梯度也很大，所以對相應部位的軸向網(wǎng)格進行加密；對幾何結構突變的位置(如噴嘴出口拐角)處的網(wǎng)格加密，以捕捉突變處的局部流動現(xiàn)象。

如圖4所示，為網(wǎng)格疏密控制方案示意圖，網(wǎng)格的疏密控制方案實施主要體現(xiàn)在剪切層網(wǎng)格的加密過程中，在面網(wǎng)格生成過程中，按照網(wǎng)格加密方案對分屬于不同塊(G、S)的剪切層網(wǎng)格進行加密處理。

圖4 網(wǎng)格疏密控制方案示意

如前所述，G塊的源面網(wǎng)格劃分采用映射法和平鋪法結合的方式，剪切層網(wǎng)格對應的部分使用了映射法，使用加密技術使得網(wǎng)格徑向尺度由內(nèi)向外逐漸變密；非剪切層區(qū)的網(wǎng)格由平鋪法生成，平鋪網(wǎng)格以映射網(wǎng)格內(nèi)邊界點為邊界元，考慮到平鋪網(wǎng)格區(qū)為非剪切層區(qū)，為得到較疏的非剪切層區(qū)網(wǎng)格，并與剪切層區(qū)網(wǎng)格的合理過渡，本文參考了網(wǎng)格細分與網(wǎng)格粗化的思想，使用了SCHNEIDERS[10]提出的網(wǎng)格過渡模板進行網(wǎng)格尺寸調(diào)節(jié)，對平鋪網(wǎng)格與映射網(wǎng)格進行過渡，如圖5所示，這一方法的優(yōu)點在于尺寸調(diào)節(jié)的能力強，調(diào)節(jié)僅對相鄰的少數(shù)單元的質(zhì)量產(chǎn)生影響。圖6a為未使用過渡方法的G塊源面網(wǎng)格圖，非剪切層區(qū)(平鋪區(qū))網(wǎng)格尺度甚至小于剪切層區(qū)(映射區(qū))，不符合疏密方案要求；圖6b為使用了過渡方法后的網(wǎng)格，非剪切層與剪切層實現(xiàn)了合理的過渡。

圖5 網(wǎng)格尺寸調(diào)節(jié)方法

圖6 網(wǎng)格比較(G塊源面)

S塊源面網(wǎng)格由映射法生成，使用網(wǎng)格加密技術對剪切層、彎折處、壁面等位置節(jié)點網(wǎng)格進行加密，隨著旋轉掃掠法的實施，加密方案傳遞到體網(wǎng)格中，如圖4中B向視圖所示。另外，為適應射流剪切層沿軸向的衰減，使用圖5b的尺寸調(diào)節(jié)技術對剪切層網(wǎng)格沿軸向進行疏化處理，使得剪切層徑向加密效果沿軸向弱化直至消失。

2.4 網(wǎng)格光順

在對S塊劃分網(wǎng)格過程中，幾何分塊能夠較好地保持網(wǎng)格與邊界的正交性，但塊與塊交接處的網(wǎng)格會產(chǎn)生彎折，彎折程度與實際幾何結構有關，這種彎折也會隨著掃掠法的實施傳遞到體網(wǎng)格中。因此有必要進行針對性的光順處理。文獻[11]提出了一種基于曲率流的四邊形主導網(wǎng)格光順方法，實現(xiàn)方法簡單，可以有效保持原網(wǎng)格的幾何特征，對邊界處網(wǎng)格的正交性影響較小。如圖7所示，為光順后與光順前網(wǎng)格對比示意圖，光順后彎折處網(wǎng)格變平滑，同時對邊界處網(wǎng)格影響較小，保持了網(wǎng)格與邊界的正交性。

圖7 彎折區(qū)光順

3 數(shù)值計算及剪切層網(wǎng)格參數(shù)的影響

對所劃分的網(wǎng)格進行數(shù)值計算，研究剪切層網(wǎng)格參數(shù)對計算結果的影響，總結剪切層網(wǎng)格參數(shù)的確定方法。

3.1 數(shù)值計算及邊界條件

對于湍流數(shù)值模擬計算，雷諾平均法是目前使用最廣泛的湍流數(shù)值模擬方法[14]，其核心是求解時均化的雷諾方程，通過湍流模型使方程封閉。湍流模型選擇基于各項異性的雷諾應力模型(RSM)，離散格式選用具有高階精度QUICK格式。

計算區(qū)域的邊界主要包括進口邊界(吸入流進口、工作流進口)、出口邊界。在進口邊界上，給定來流速度，出口邊界上給定出口壓力。數(shù)值計算采用FLUENT求解器，計算殘差值曲線如圖8所示。

圖8 數(shù)值計算收斂殘差

3.2 剪切層網(wǎng)格參數(shù)定義

剪切層網(wǎng)格參數(shù)包括底層網(wǎng)格尺寸和漸變率(δ)，如圖9～10所示?？芍羟袑泳W(wǎng)格可看作是上游壁面邊界層網(wǎng)格的延伸，因此可采用上游壁面邊界層的網(wǎng)格尺度Y+來表示剪切層底層網(wǎng)格尺度。Y+表示近壁第1層網(wǎng)格點至壁面的無量綱距離，表達式為[12]：

(1)

式中：τw為壁面剪應力；Δn為第1層網(wǎng)格點距離壁面的實際距離；ρ為流體密度；v為流體的運動黏度。

由式(1)可知，Y+值是由流體物性參數(shù)(ρ，V)和壁面剪切力(τw)共同決定的，Y+值與近壁底層網(wǎng)格尺度Δn是對應的。在計算過程中，應通過不斷調(diào)整底層網(wǎng)格尺度Δn，使得Y+值滿足湍流模型的要求[13]。

圖9 剪切層網(wǎng)格參數(shù)示意

圖10 網(wǎng)格數(shù)量與剪切層網(wǎng)格尺度關系

圖11為數(shù)值計算得到的速度等值線圖，圖中標識了(a)、(b)、(c)、(d)4個區(qū)域，分別對應了剪切層形成區(qū)(射流初始區(qū))、射流勢流核區(qū)(末端)、喉管回流區(qū)、擴散段充分發(fā)展區(qū)，通過對比各區(qū)域計算結果，分析剪切層網(wǎng)格參數(shù)對數(shù)值計算的影響。

圖11 速度等值線圖

圖12 應變率梯度預報結果比較((a)區(qū)，δ=1.1)

a (b)區(qū)速度剖面

b (c)區(qū)速度剖面

3.4 剪切層網(wǎng)格參數(shù)的確定方法總結

4 算例驗證及適應性分析

網(wǎng)格1與網(wǎng)格4具有相同的網(wǎng)格節(jié)點數(shù)量，后者以前者為基礎進行加密。由表2可知，在同樣的網(wǎng)格數(shù)量下，網(wǎng)格1的計算結果誤差較大，誤差高于4.5%。

圖14 算例網(wǎng)格

D0/mmS/mm2l1/mml2/mm12．110．924280

表2 計算結果與試驗結果比較

關于不同湍流模型的預報結果，在小流量比(0.81)下，RSM模型的計算誤差最小，雙方程模型的計算誤差略高于RSM模型；在大流量比(1.71)下，RSM模型保持了預報的精確性，而雙方程模型的計算誤差明顯上升，這與文獻[16]得到的結論一致。

如圖15所示，喉管段的速度分布呈現(xiàn)不對稱性，吸入管偏置造成了這種不對稱特性。但這種不對稱性在應變率(變量梯度)分布中卻不明顯。

圖15 喉管段速度剖面(流量比1.71)

如圖16所示，應變率分布的對稱性要好于速度分布，另外，各截面應變率的峰值均出現(xiàn)在在G、S塊交界面處或附近，且隨著流場沿軸向的發(fā)展，各峰值迅速下降。該種分布規(guī)律一方面說明軸對稱式的剪切層網(wǎng)格配置方案是合適的，另一方面也說明沒有必要將射流剪切層加密方案沿軸向延伸至太長，網(wǎng)格劃分方案中針對剪切層的網(wǎng)格的軸向過渡處理是必要的。

圖16 應變率分布(流量比1.71)

5 結論

1) 根據(jù)中心射流泵計算域幾何結構特點，采用掃掠法劃分六面體結構化網(wǎng)格。通過幾何分塊保證網(wǎng)格與邊界的正交性，對整體網(wǎng)格進行光順處理以消除分塊帶來的網(wǎng)格彎折效應。

2) 結合映射網(wǎng)格的拉伸方法提出邊界及剪切層網(wǎng)格加密技術，對射流剪切層、壁面邊界層以及幾何彎折部位進行網(wǎng)格加密，以適應流場梯度的變化。

3) 利用結構化網(wǎng)格關聯(lián)性強的特點，提出使用噴嘴內(nèi)壁無量綱距離Y+來代表剪切層底層網(wǎng)格尺度，并通過數(shù)值計算分析總結了剪切層網(wǎng)格參數(shù)的確定方法。

4) 應用網(wǎng)格劃分方法對某型射流泵進行建模和數(shù)值計算，選擇不同的湍流模型，得到的計算結果與試驗值吻合良好，從而驗證了該網(wǎng)格劃分方法的正確性。該網(wǎng)格劃分方法對有限空間射流流場的網(wǎng)格劃分也具有一定的參考價值。

[1] 金博.射流排砂泵動態(tài)特性分析[J].石油礦場機械，2012，41(7)：55-58.

[2]FanJ，EvesJ，ThompsonHM，etal.Computationalfluiddynamicanalysisanddesignoptimizationofjetpumps[J].Computers&Fluids，2011(46)：212-217.

[3]Xue-GuanSong，Joon-HongPark，Seung-GyuKim，etal.Performancecomparisonanderosionpredictionofjetpumpsbyusinganumericalmethod[J].Mathematicalandcomputermodelling，2013(57)：245-253.

[4] 陳鳳官，綦耀光，王渭，等.環(huán)空射流泵排采工藝在煤層氣井中的適用性研究[J].石油礦場機械，2012，41(6)：46-49.

[5]ThompsonJF，ThamesFC，MastinCW.Automaticnumericalgenerationofbody-fittedcurvilinearcoordinatesystemfieldcontaininganynumberofarbitrarytwo-dimensionbodies[J].JournalofComputationalPhysics，1974(15)：299-319.

[6]TedDBlacker，MichaelBStephenson.Paving：Anewapproachtoautomatedquadrilateralmeshgeneration[J].InternationalJournalforNumericalMethodsinEngineering，1991(32)：811-847.

[7] 沈海龍，蘇玉民.肥大型船伴流場數(shù)值模擬的網(wǎng)格劃分方法研究[J].哈爾濱工程大學學報，2008，29(11)：1191-1198.

[8] 龍新品，朱勁木.射流泵內(nèi)部流動的數(shù)值模擬[J].武漢大學學報(工學版)，2002，35(6)：1-6.

[9] 龍新平，陸宏圻，李煒.有限空間射流的數(shù)值模擬[J].應用基礎與工程科學學報，1997，5(1)：76-83.

[10]SCHNEIDERSR.Refiningquadrilateralandhexahedralelementmeshes[C]//Proceedingof5thInternationalConferenceonNumericalGridGenerationinComuputationalFieldSimulationsCinton，USA，April1996.Clinton：MississippiStateUniversity，1996：679-688.

[11] 胡事民，來煜坤，楊永亮.基于曲率流的四邊形主導網(wǎng)格的光順方法[J].計算機學報，2008，31(9)：1622-1628.

[12]GeorgiKalitzin，GorazdMedic，Gianlucalaccarino，etal.Near-wallbehaviorofRANSturbulencemodelsandimplicationsforwallfunctions[J].Journalofcomputationalphysics，2005(204)：265-291.

[13] 李曉俊，袁壽其，潘中永，等.離心泵邊界層網(wǎng)格的實現(xiàn)及應用評價[J].農(nóng)業(yè)工程學報，2012，20(28)：67-72.

[14] 覃文潔，胡春光，郭良平，等.近壁面網(wǎng)格尺寸對湍流計算的影響[J].北京理工大學學報，2006，26(5)：388-392.

[15] 蔡標華.射流泵初生空化及其試驗研究[D].武漢：武漢大學動力與機械學院，2005.

[16] 周凌九，袁玲麗.射流泵內(nèi)部流動計算中不同湍流模擬方法的比較[J].排灌機械工程學報，2013，31(1)：25-30.

Research and Application of Mesh Partition Method for Numerical Simulation of Jet Pump

ZHANG Di1，DONG Xiangwei2，LI Zengliang2，LU Shengnan1

(1.OffshoreOilEngineeringCo.，Ltd.，Tianjin300451，China；2.ChinaUniversityofPetroleum(Huadong)，Qingdao266580，China)

A structured partition method for computational domain of jet pump is researched，in which block division and smooth techniques are applied to control the mesh orthogonality.Mesh density scheme is established with jet shear layer at the core and the influence of mesh parameter of shear layer on numerical calculation.Mesh partition methods appropriate for numerical calculation of jet pump were obtained by comparing the results between numerical and experimental results.The results show that，sweeping method is appropriate for jet pump simulation and it benefits the implementation process of mesh density scheme.It also benefits the orthogonality and relevance，which is conducive to parameter research.Grid size of first layer and transition ratio determine quality and quantity of the mesh.Simulation results of external characteristic have little relationship with size of first layer.The results indicated that this mesh partition method is effective for flow field simulation of jet pump.

jet pump；mathematical simulation；structured grid；spacing control

2016-05-17

張 棣(1986-)，男，山東濱州人，工程師，博士，主要從事海洋石油裝備研究，E-mail：zd739@126.com。

1001-3482(2016)11-0001-07

TE

A

10.3969/j.issn.1001-3482.2016.11.001