邊界層內(nèi)水介質(zhì)在壓差驅(qū)動(dòng)下沿小孔向氣腔射流問(wèn)題研究

鮑文春,劉元清,李 巖,武龍龍

(北京宇航系統(tǒng)工程研究所,北京100076)

邊界層內(nèi)水介質(zhì)在壓差驅(qū)動(dòng)下沿小孔向氣腔射流問(wèn)題研究

鮑文春,劉元清,李 巖,武龍龍

(北京宇航系統(tǒng)工程研究所,北京100076)

對(duì)邊界層內(nèi)小孔氣水多相流場(chǎng)下射流問(wèn)題開(kāi)展數(shù)值仿真及定常水洞試驗(yàn)研究,建立了適用于邊界層內(nèi)壓差驅(qū)動(dòng)下小孔向氣腔射流多相流場(chǎng)問(wèn)題研究的數(shù)值仿真計(jì)算模型,針對(duì)典型孔參數(shù)及氣水流場(chǎng)條件,對(duì)比分析了仿真試驗(yàn)數(shù)據(jù),驗(yàn)證了數(shù)值仿真模型的正確性及模型計(jì)算精度。結(jié)合流體質(zhì)點(diǎn)受力及運(yùn)動(dòng)模型及平板邊界層理論,分析了氣水域壓力場(chǎng)特征及水域流動(dòng)規(guī)律對(duì)小孔射流過(guò)程的作用機(jī)理及影響規(guī)律,開(kāi)展了孔參數(shù)對(duì)射流多相流場(chǎng)特征及射流量的影響研究。獲得了小孔射流量估算方法,為航行體上防水裝置設(shè)計(jì)提供數(shù)據(jù)支撐。

小孔;射流;多相流;邊界層;水洞試驗(yàn)

0 引言

在水下高速運(yùn)動(dòng)航行體工程設(shè)計(jì)和研究中,經(jīng)常采用通氣技術(shù)獲得良好流體動(dòng)力特性[1]。通氣技術(shù)可以分為主動(dòng)通氣和被動(dòng)通氣兩類,被動(dòng)通氣又稱為排氣技術(shù)。采用排氣技術(shù)的航行體內(nèi)部留有儲(chǔ)氣腔,并在表面特定位置開(kāi)設(shè)一定數(shù)量的氣縫/孔 (見(jiàn)圖1)。發(fā)射前在氣腔內(nèi)預(yù)置一定壓力的氣體,在航行體向水面運(yùn)動(dòng)過(guò)程中,外界環(huán)境壓力逐漸減小,預(yù)置氣體在內(nèi)外壓差作用下從排氣孔/縫排出,形成氣泡 (氣膜)附著在航行體表面。通過(guò)氣泡或氣膜的形態(tài)進(jìn)行航行體運(yùn)動(dòng)過(guò)程中流體動(dòng)力調(diào)節(jié)[2]。

圖1 航行體表面排氣孔結(jié)構(gòu)示意圖Fig.1 Sketch map small exhaust holes on the under voyage

航行體水下彈射過(guò)程屬于非定常、非線性、跨介質(zhì)的多相流動(dòng)問(wèn)題,航行體采用高壓燃?xì)鈴椛涑鐾?燃?xì)饨橘|(zhì)在發(fā)射筒口發(fā)展演化,導(dǎo)致筒口區(qū)域的介質(zhì)組成及力學(xué)環(huán)境復(fù)雜多變。當(dāng)航行體上小孔結(jié)構(gòu)隨航行體運(yùn)動(dòng)到筒口附近時(shí),復(fù)雜的多相流場(chǎng)環(huán)境將使小孔面臨進(jìn)水問(wèn)題,或者稱為水介質(zhì)沿小孔向氣腔內(nèi)射流問(wèn)題,進(jìn)水過(guò)程將給航行體上儀器設(shè)備正常工作帶來(lái)不可預(yù)估的風(fēng)險(xiǎn)。國(guó)內(nèi)外學(xué)者針對(duì)小孔射流問(wèn)題進(jìn)行了一定研究,但研究成果主要集中在氣膜冷卻問(wèn)題及壓差驅(qū)動(dòng)下小孔向氣腔內(nèi)射流問(wèn)題,針對(duì)氣腔進(jìn)水問(wèn)題的研究資料較少。李廣超等[4]等對(duì)發(fā)動(dòng)機(jī)葉片的氣膜冷卻問(wèn)題進(jìn)行了數(shù)值仿真研究,指出了圓柱孔射流冷卻的有害渦流動(dòng)結(jié)構(gòu),論述了幾何結(jié)構(gòu)和氣動(dòng)參數(shù)對(duì)氣膜冷卻特性的影響;戴萍等[5]研究了不同吹風(fēng)比下入射角度對(duì)縮放槽縫孔氣膜冷卻效果的影響;劉海涌等[6]對(duì)渦輪葉片錯(cuò)排出流孔流場(chǎng)進(jìn)行了試驗(yàn)研究,獲得通道和射流孔側(cè)壁靜壓的變化規(guī)律;李芳等[7]對(duì)流體在多孔仿生射流表面上的流動(dòng)特性進(jìn)行了數(shù)值模擬,探討了仿生射流表面減阻機(jī)理;施瑤等[8]采用LES與VOF相結(jié)合的方法對(duì)多載荷AUV后平衡艙空腔的進(jìn)水過(guò)程進(jìn)行了數(shù)值模擬,考察了進(jìn)水孔直徑、進(jìn)水孔位置以及航行深度等因素對(duì)進(jìn)水過(guò)程的影響,相關(guān)研究成果主要針對(duì)大孔徑進(jìn)水孔。本文對(duì)運(yùn)動(dòng)航行體邊界層內(nèi)水介質(zhì)沿小孔徑向航行體內(nèi)部氣域射流問(wèn)題進(jìn)行研究,并分析孔參數(shù)及內(nèi)外壓條件等影響因素對(duì)射流流體體積的影響規(guī)律,為小孔射流總量估算及航行體上防水裝置設(shè)計(jì)提供數(shù)據(jù)支撐。

1 數(shù)值計(jì)算模型

1.1 多相流模型

在研究氣水多相流場(chǎng)射流問(wèn)題時(shí),氣水流體介質(zhì)具有清晰的相界面,相界面的分布情況直接影響流場(chǎng)中的物性分布及氣水介質(zhì)間的相互作用力。為精確構(gòu)造非定常流場(chǎng)中的相界面,采用PLIC[9]技術(shù)進(jìn)行界面捕捉。在多相流模型方面, VOF[10-11]類方法因?yàn)榭梢暂^好保持流體質(zhì)量守恒,且能輕易處理界面拓?fù)浣Y(jié)構(gòu)的變化,受到很多學(xué)者的青睞。本文采用VOF方法來(lái)構(gòu)造相界面。VOF多相流模型中的控制方程形式如下:

混合物連續(xù)方程:

其中,αg為空氣相的體積分?jǐn)?shù);u i、u j為速度分量;ρm為混合物密度,ρm=(1-αg)ρl+αgρg,ρl、ρg分別為水的密度及空氣的密度;μm為混合物的動(dòng)力黏性系數(shù),μm=(1-αg)μl+αgμg,μl、μg分別為水和空氣的動(dòng)力黏性系數(shù),μt為湍流黏性系數(shù)。

氣相體積分?jǐn)?shù)方程:

計(jì)算過(guò)程中先求解氣相的體積分?jǐn)?shù),然后根據(jù)單元內(nèi)體積分?jǐn)?shù)守恒條件來(lái)獲得主相的體積分?jǐn)?shù)及混合物密度等。

1.2 湍流模型

采用標(biāo)準(zhǔn)κ-ε湍流模型[12],湍流動(dòng)能κ及耗散率ε的控制方程分別為:

其中,σt=ρmCμκ2/ε為湍流黏性系數(shù);經(jīng)驗(yàn)常數(shù)Cμ=0.09,Cε1=1.44,Cε2=1.42;σε=1.3,σk=1.0表征基于湍流動(dòng)能及耗散率的湍流普朗特?cái)?shù);為由平均速度梯度引起的湍流動(dòng)能生成項(xiàng)。

1.3 計(jì)算域及邊界條件設(shè)置

本文計(jì)算采用三維數(shù)值仿真計(jì)算模型,計(jì)算域尺度及邊界條件定義方式見(jiàn)圖2。數(shù)值離散方法采用有限體積法,湍流輸運(yùn)方程采用1階差分格式,時(shí)間離散為1階顯式。壓力速度耦合方式采用SIMPLEC算法,壓力及動(dòng)量方程離散采用2階迎風(fēng)格式。

圖2 計(jì)算域及邊界條件定義Fig.2 Computation region and boundary conditions

2 數(shù)值計(jì)算有效性驗(yàn)證

為驗(yàn)證該方法計(jì)算結(jié)果的準(zhǔn)確程度,在中尺度水洞中開(kāi)展小孔氣水射流問(wèn)題試驗(yàn),測(cè)量獲得一定壓差條件下小孔的射流質(zhì)量流率。試驗(yàn)設(shè)施及測(cè)量設(shè)備情況見(jiàn)圖3。

針對(duì)直徑D為6mm,孔深H為10mm小孔,水域流速為3.7m/s,水流與壁面之間攻角為0°,在不同氣水流場(chǎng)外壓條件下進(jìn)行小孔射流過(guò)程仿真計(jì)算及水洞試驗(yàn)。

圖3 小孔氣水射流水洞試驗(yàn)整體方案Fig.3 The plane of experimental of small holes jet

數(shù)值仿真計(jì)算獲得的小孔射流體積流率與相應(yīng)工況試驗(yàn)測(cè)量工況的對(duì)比情況見(jiàn)表1。不同射流條件下,數(shù)值仿真獲得的小孔射流體積流率隨氣水流域壓比變化規(guī)律與試驗(yàn)測(cè)量結(jié)果基本一致,數(shù)值仿真獲得的進(jìn)水量大于試驗(yàn)測(cè)量值,二者相對(duì)誤差約為29%。產(chǎn)生上述差異的主要原因有兩點(diǎn):第一,因孔面積及孔內(nèi)外壓差均較小,試驗(yàn)中對(duì)小進(jìn)水量的測(cè)量存在一定的差異;第二,基于RANS方程的數(shù)值模型對(duì)小孔邊界層內(nèi)的流動(dòng)不能完全模擬,但上述差異滿足工程設(shè)計(jì)的需求。故初步驗(yàn)證了本文建立的小孔氣水射流數(shù)值仿真計(jì)算模型的正確性及適用性。

表1 射流體積流率值仿真與試驗(yàn)數(shù)據(jù)對(duì)比Tab.1 Comparisons of jet volume flow rate in numerical results and experimentals

3 流場(chǎng)特征對(duì)小孔氣水射流流動(dòng)機(jī)理分析

3.1 射流驅(qū)動(dòng)力影響分析

為分析流體沿小孔向氣域射流機(jī)理問(wèn)題,對(duì)不同氣水域外壓條件作用下小孔射流問(wèn)題進(jìn)行仿真計(jì)算。

表2給出了水域無(wú)運(yùn)動(dòng)時(shí)小孔射流控制參數(shù)對(duì)比結(jié)果,圖4給出了相應(yīng)的小孔射流多相流場(chǎng)情況。表明水域無(wú)運(yùn)動(dòng)時(shí),射流多相流場(chǎng)特征由小孔兩側(cè)氣水介質(zhì)壓比決定。

表2 水域無(wú)運(yùn)動(dòng)工況射流控制參數(shù)對(duì)比Tab.2 Comparisons of jet parameters with no movement in water region

圖4 水域無(wú)運(yùn)動(dòng)時(shí)氣水射流速度場(chǎng)分布情況對(duì)比Fig.4 Comparisons of velocity distribution of gas-water jet with no movement in water region

表3給出了水域有運(yùn)動(dòng)時(shí)孔射流控制參數(shù)對(duì)比結(jié)果,圖5給出了相應(yīng)的小孔射流多相流場(chǎng)情況。表明水域有運(yùn)動(dòng)時(shí),射流多相流場(chǎng)特征由小孔兩側(cè)氣水介質(zhì)壓差決定。

表3 水域有運(yùn)動(dòng)工況射流控制參數(shù)對(duì)比(U0=3.7m/s)Tab.3 Comparisons of jet parameters with movement in water region

圖5 水域運(yùn)動(dòng)速度U0=3.7m/s氣水射流速度場(chǎng)分布情況對(duì)比Fig.5 Comparisons of velocity distribution of gas-water jet with movement in water region is U0=3.7m/s

數(shù)值仿真結(jié)果表明,有無(wú)水域運(yùn)動(dòng)速度下小孔氣水射流問(wèn)題具有不同的流動(dòng)機(jī)理及多相流場(chǎng)特征,如圖6所示。無(wú)來(lái)流情況下,射流充滿整個(gè)射流孔;有來(lái)流情況下,水流僅沿小孔部分區(qū)域進(jìn)入氣域。分析認(rèn)為流動(dòng)邊界層的有無(wú)是上述射流景象差異的主要原因,外流場(chǎng)有來(lái)流情況下,孔附近存在發(fā)展較為充分的流動(dòng)邊界層。

圖6 小孔氣水射流景象(壓比=0.995,壓差0.5kPa)Fig.6 Sketch of the gas-water jet in small holes

3.2 水域運(yùn)動(dòng)參數(shù)對(duì)小孔射流影響分析

針對(duì)直徑D為6mm,孔深H為10mm小孔,開(kāi)展不同水域來(lái)流條件下,小孔射流過(guò)程仿真計(jì)算及水洞試驗(yàn),分析水域運(yùn)動(dòng)參數(shù)的影響機(jī)理及作用規(guī)律分析。

綜合分析所有進(jìn)水的工況,得到了腔內(nèi)外壓差與進(jìn)水量之間的初步關(guān)系。定義考慮動(dòng)壓的表面壓力系數(shù):

其中,Δp為氣水域壓強(qiáng)差,U為水域流速,α為水域流動(dòng)相對(duì)壁面攻角,ρ為水的密度。

從水流速度為6.5m/s和9m/s的工況可以看到,進(jìn)水量隨壓力系數(shù)Cp線性增加。而來(lái)流速度等于3.7m/s時(shí)進(jìn)水量隨壓力系數(shù)增加而增加,最后趨于一個(gè)定值。壓力系數(shù)相同時(shí),攻角為0°和6°進(jìn)水量差別不大,而攻角等于9°時(shí)進(jìn)水量明顯增加。

圖7 表面壓力系數(shù)與進(jìn)水量的關(guān)系Fig.7 Relationship between surface pressure coefficient and water inflow

從圖7中可以看到速度3.7m/s與6.5m/s和9m/s有明顯的分界,下面將從流體邊界層理論來(lái)定性分析該現(xiàn)象。物體在黏性流體運(yùn)動(dòng)時(shí)表面有一層很薄的邊界層。現(xiàn)定義坐標(biāo)系原點(diǎn)在翼型頭部與孔中心高度一致,沿水流方向?yàn)閄正方向,沿孔深方向指向外部為Y正 (如圖8所示)。雷諾數(shù)定義為:

其中,V x為x處速度(文中取水流流動(dòng)速度U);L為孔中心距原點(diǎn)的x方向坐標(biāo),文中取值為77.2mm;v為流體運(yùn)動(dòng)黏度,文中取值為1.007×10-4m2/s。

圖8 坐標(biāo)系定義Fig.8 Definition of the coordinate system

來(lái)流速度為3.7m/s時(shí),雷諾數(shù)Re=2.84× 105,小孔處于層流邊界層內(nèi);來(lái)流速度等于6.5m/s時(shí),Re=4.98×105。計(jì)算表明來(lái)流速度v＞6.5m/s時(shí),小孔附近為湍流邊界層。由平板邊界層理論可知[12]:由層流轉(zhuǎn)變?yōu)槲闪鞯呐R界雷諾數(shù)Re=5×105~3×106。此時(shí)從前駐點(diǎn)開(kāi)始的一段距離為層流邊界層,隨著離前駐點(diǎn)距離增大,邊界層逐漸變厚。

圖9為壁面附近邊界層內(nèi)單個(gè)微元流體質(zhì)點(diǎn)c流經(jīng)孔口附近的運(yùn)動(dòng)示意圖,以及忽略流體質(zhì)點(diǎn)剪切變形的受力分析。假設(shè)流體質(zhì)點(diǎn)在Δx和Δy上無(wú)能量耗散,所以

式(10)表明流體質(zhì)點(diǎn)c主要受到Y(jié)方向的作用力,其大小約為

其中,ΔA表示流體質(zhì)點(diǎn)Y方向的微元面積。流體質(zhì)點(diǎn)流經(jīng)孔口的時(shí)間Δt為:

其中,d y為流線位置孔寬,m c為微元流體質(zhì)點(diǎn)質(zhì)量,E x為質(zhì)點(diǎn)動(dòng)能。

圖9 流體質(zhì)點(diǎn)運(yùn)動(dòng)及受力示意圖Fig.9 Sketch of movement and force for fluid particle

由牛頓第二定律知流體質(zhì)點(diǎn)向孔口偏移距離Δd為

單位時(shí)間內(nèi)的射流量為

其中C為修正常系數(shù),由坐標(biāo)定義可知

其中Mc表示流體質(zhì)點(diǎn)的動(dòng)量,對(duì)式(16)積分得

其中v(y)表示壁面速度分布規(guī)律。

圖10是孔附近壁面邊界層流場(chǎng)示意圖,圖中3條流速拋面線,從左到右依次表示來(lái)流速度等于3.7m/s、6.5m/s和9m/s。當(dāng)來(lái)流速度為3.7m/s時(shí),空口附近雷諾數(shù)沒(méi)有達(dá)到臨界雷諾數(shù),壁面附近邊界層為層流邊界層。邊界層內(nèi) (δ)的流體質(zhì)點(diǎn)與邊界層外 (δ+)的流體質(zhì)點(diǎn)幾乎沒(méi)有相互交換,且貼近水翼壁面的流體質(zhì)點(diǎn)摩擦能量耗散,此時(shí)邊界層內(nèi)流體質(zhì)點(diǎn)動(dòng)能較小,在相同壓差阻力作用下較易改變運(yùn)動(dòng)狀態(tài)。而來(lái)流速度U大于6.5m/s時(shí),雷諾數(shù)接近臨界雷諾數(shù),而一般轉(zhuǎn)捩點(diǎn)位置都非常靠前,故此時(shí)只有孔附近邊界層發(fā)展為劇烈的紊流邊界層,邊界層內(nèi)的流體質(zhì)點(diǎn)與邊界層外的流體質(zhì)點(diǎn)交換劇烈,即邊界層內(nèi)的流體質(zhì)點(diǎn)保持較大動(dòng)能,此時(shí)在相同壓差阻力作用下較難改變運(yùn)動(dòng)狀態(tài)。另一方面,孔口孔徑相同,流速小的流體質(zhì)點(diǎn)經(jīng)過(guò)孔口所用的時(shí)間Δt長(zhǎng),反之時(shí)間Δt短。

圖10 邊界層流場(chǎng)示意圖Fig.10 Sketch of the boundary layer flow

綜上分析可知,壓差相同時(shí)來(lái)流速度小于臨界雷諾數(shù)需要的速度時(shí),壁面附近邊界層為層流邊界層,流體質(zhì)點(diǎn)水平動(dòng)能E x小,即單個(gè)微元流體質(zhì)點(diǎn)流過(guò)孔口偏移距離Δd越大。壁面的水介質(zhì)由無(wú)數(shù)微元流體質(zhì)點(diǎn)組成,從概率上來(lái)講E x越小,進(jìn)入孔內(nèi)的流體質(zhì)點(diǎn)越多,宏觀表現(xiàn)是流體射流量越大。同理可知,水域流動(dòng)速度相同時(shí),壓差Δp越大,流體射流量越多。

4 孔參數(shù)對(duì)氣水射流流動(dòng)機(jī)理分析

4.1 小孔徑深比影響分析

在上述研究基礎(chǔ)上,進(jìn)一步分析小孔徑深比H/D(徑深比定義為射流孔深度與小孔直徑之比)對(duì)射流過(guò)程的影響。對(duì)以下徑深比的小孔射流過(guò)程進(jìn)行仿真研究,見(jiàn)表4。

表4 小孔徑深比影響數(shù)值仿真工況Tab.4 Effect of diameter-depth ratio on numerical simulation

圖11給出了不同徑深比下小孔射流速度隨壓比變化規(guī)律,結(jié)果表明,薄壁孔的射流面平均流速相對(duì)較大。

圖11 徑深比對(duì)射流速度影響Fig.11 Effect of diameter-depth ratio on the velocity of jet flow

不同徑深比下小孔射流質(zhì)量流率隨氣水介質(zhì)壓比變化規(guī)律見(jiàn)圖12,研究結(jié)果表明,小孔射流質(zhì)量流率不受徑深比影響。

圖13給出了不同徑深比下多相流場(chǎng)特征。可見(jiàn):徑深比為3.00及1.67時(shí),小孔射流的流動(dòng)特征相近,在小孔內(nèi)壁存在局部回流渦,氣流沿小孔下緣逆向水流方向進(jìn)入氣域;當(dāng)徑深比為0.33時(shí),小孔內(nèi)部的局部流動(dòng)渦消失,氣流沿小孔下緣順?biāo)蛄鲃?dòng)方向進(jìn)入氣域。

4.2 小孔面積影響分析

保持小孔徑深比不變,開(kāi)展射流孔面積影響分析,設(shè)置如表5所示仿真研究工況。

圖12 小孔徑深比對(duì)射流質(zhì)量流率影響規(guī)律Fig.12 Effect of diameter-depth ratio on the mass flow rate of jet

圖13 徑深比對(duì)射流影響分析(P氣/P水=0.995)Fig.13 Analysis of the diameter-depth ratio on the jet flow

表5 小孔面積影響研究數(shù)值仿真工況_Tab.5 Effect of the small holes on numerical simulations_

圖14給出了射流體積流率隨氣水介質(zhì)壓比變化規(guī)律。不同排氣面積下,射流體積流率隨壓比變化規(guī)律一致,射流體積流率隨小孔面積增加而增大。

圖14 小孔面積對(duì)射流質(zhì)量流率影響Fig.14 Effect of the area of small holes on the velocity of jet flow

圖15給出了射流體積流率與排氣面積之比。不同孔徑下,射流體積流率與孔面積之比為常值,即射流體積流率隨孔面積成正比。

圖15 射流體積流率與排氣面積之比Fig.15 The ratio between jet volume flow rate and ventilation area

5 結(jié)論

本文針對(duì)邊界層內(nèi)小孔氣水多相流場(chǎng)下射流問(wèn)題開(kāi)展數(shù)值仿真及定常水洞試驗(yàn)研究,對(duì)比相同工況下數(shù)值仿真及水洞試驗(yàn)結(jié)果,驗(yàn)證了數(shù)值模型的正確性及模型計(jì)算精度。并針對(duì)不同氣水兩側(cè)流場(chǎng)外壓條件、不同水域運(yùn)動(dòng)狀態(tài)、不同小孔類型及小孔面積下小孔射流問(wèn)題開(kāi)展數(shù)值仿真研究,獲得了小孔進(jìn)水量工程估算公式,為航行體上防進(jìn)水結(jié)構(gòu)的設(shè)計(jì)提供數(shù)據(jù)支撐,本文獲得如下結(jié)論:

1)水域無(wú)運(yùn)動(dòng)時(shí),氣水射流問(wèn)題主要受小孔兩側(cè)壓比驅(qū)動(dòng),水域有運(yùn)動(dòng)時(shí),氣水射流問(wèn)題由孔兩側(cè)氣水介質(zhì)壓差驅(qū)動(dòng),小孔射流量隨著氣水壓差的增加而增大;

2)小孔射流過(guò)程可歸納為層流邊界層內(nèi)射流狀態(tài)及湍流邊界層內(nèi)射流狀態(tài),相同氣水流場(chǎng)外壓條件下,處于湍流邊界層內(nèi)的小孔射流體積流率更小;

3)從小孔徑深比角度來(lái)看,壓差驅(qū)動(dòng)下小孔氣水多相流場(chǎng)射流問(wèn)題,射流流量不受徑深比影響。但射流孔局部流動(dòng)特征可分為兩種典型射流,即薄壁孔(H/D＜0.5)射流及深孔(H/D＞2)射流。對(duì)于深孔射流,小孔內(nèi)壁存在局部回流渦,隨著孔徑深比的減小,回流渦消失;

4)相同來(lái)流條件下,射流體積流率隨小孔面積增加而增大,且射流體積流率與孔面積之比為常值。

[1] 李杰,魯傳敬,陳鑫,等,附著空泡對(duì)潛射導(dǎo)彈彈道影響分析[J].彈道學(xué)報(bào),2014,26(3):54-58.

[2] 鮑文春,權(quán)曉波,魏海鵬.航行體排氣水下發(fā)射流體動(dòng)力數(shù)值仿真研究 [J].導(dǎo)彈與航天運(yùn)載技術(shù), 2014(5):14-18.

[3] 陳浮,馬貴輝,程少華,等.直、斜孔排氣對(duì)航行體繞流流動(dòng)影響:Part 1——流場(chǎng)結(jié)構(gòu)[J].工程熱物理學(xué)報(bào),2016,37(3):507-513.

[4] 李廣超,柏樹(shù)生,吳冬,等.氣膜孔形狀對(duì)渦輪葉片氣膜冷卻影響的研究進(jìn)展 [J].熱能動(dòng)力工程, 2010,25(6):581-585.

[5] 戴萍,林楓.入射角度對(duì)縮放槽縫孔氣膜冷卻效果的影響[J].汽輪機(jī)技術(shù),2010,52(4):263-267.

[6] 劉海涌,孔滿昭,劉松齡,等.有錯(cuò)排射流沖擊的受限長(zhǎng)通道及出流孔流場(chǎng)的實(shí)驗(yàn)研究 [J].推進(jìn)技術(shù), 2006,27(4):307-311.

[7] 李芳,趙剛,劉維新,等.多孔仿生射流表面減阻特性數(shù)值模擬[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版),2014, 44(6):1698-1703.

[8] 施瑤,潘光.多載荷AUV進(jìn)水過(guò)程的LES/VOF數(shù)值模擬[J].西北工業(yè)大學(xué)學(xué)報(bào),2015(3):438-443.

[9] 郭建紅,魯傳敬,陳瑛,等.基于高階和高分辨率格式的自然空泡流數(shù)值模擬 [J].水動(dòng)力學(xué)研究與進(jìn)展,2009,24(2):224-231.

[10] 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.

[11] 李國(guó)杰,黃萌,陳斌.基于PISO算法的非結(jié)構(gòu)化網(wǎng)格VOF算法[J].工程熱物理學(xué)報(bào),2013,34 (3):476-479.

[12] 王聰,何春濤,權(quán)曉波,等.空氣壓強(qiáng)對(duì)垂直入水空泡影響的數(shù)值研究 [J].哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2012,44(5):14-19.

[13] 張兆順.湍流[M].北京:國(guó)防工業(yè)出版社,2002: 199-242.

Study on Jets of Gas-water Multiphase Flow in Small Holes Drived by Pressure Difference within the Boundary Layer

BAO Wen-chun,LIU Yuan-qing,LI Yan,WU Long-long
(Beijing Institute of Astronautical Systems Engineering,Beijing 100076,China)

In this paper,numerical simulations and steady water tunnel experimentals are carried out to study the jets of gas-water multiphase flow in small holes within the boundary layer.A computational model is constructed on the gas jet multiphase flow.Based on the comparison of numerical and experimental results,the correctness and accuracy of numerical model are verified by using typical parameters of small holes and conditions of gas-water flow.Combining the model of force and movement for fluid particle with the theory of plate boundary layer,the characteristics of pressure in the gas-water region is analyzed and the effect of the water region conditions on the process of jets in small holes is studied.Meanwhile,the influence of paremeters of small holes on the multiphase flowfield characteristics and jet flow is also studied.The study results of this article can provide methods for the gasjet mass flow rate calculation and water-seal devices design.

Small holes;Jets;Multiphase flow;Boundary layer;Water tunnel experimentals

O375.5

A

2096-4080(2017)04-0029-09

2017-07-28;

2017-10-11

裝備發(fā)展部共用技術(shù)基金(41406040402)

鮑文春(1988-),女,碩士,工程師,研究方向?yàn)榱黧w力學(xué)。E-mail:baowenchun@126.com