任澤宇,孫龍泉,姚熊亮,趙紀(jì)鵬

(哈爾濱工程大學(xué)船舶工程學(xué)院, 哈爾濱 150000)

0 引言

水下航行體高速運動過程中，覆蓋于航行體表面的自然空泡經(jīng)歷初生、發(fā)展以及出水潰滅，改變了航行體表面的載荷分布，同時沖擊、腐蝕航行體表面，威脅航行體的結(jié)構(gòu)及彈道穩(wěn)定性[1]，改善自然空泡的負(fù)面影響成為兵器制造的核心問題之一[2]。通過注入不可凝氣體形成覆蓋于航行體表面的通氣空泡，即主動通氣技術(shù)，可增加氣泡穩(wěn)定性，降低航行體阻力，并調(diào)節(jié)航行體的姿態(tài)參數(shù)[3-4]。

通氣空泡流在本質(zhì)上是一種非線性非穩(wěn)態(tài)湍動的復(fù)雜多相流動問題[5]，國內(nèi)外學(xué)者在此方面進行了諸多研究。在通氣空泡發(fā)展過程中，楊武剛等[6]和張紀(jì)華等[7]發(fā)現(xiàn)通氣量是影響空泡幾何尺寸和表面穩(wěn)定性的重要因素，而空泡界面的波動和空泡的震蕩是影響通氣空泡表面不穩(wěn)定的主要因素[8]。龔瑞巖[9]通過自主設(shè)計的垂直約束式水下發(fā)射平臺開展水下航行體單排氣孔通氣試驗，對通氣空泡的流動形態(tài)變化規(guī)律進行了總結(jié)。Li等[5]采用邊界元法對三維通氣空泡進行了研究，結(jié)果表明入流速度是影響氣泡發(fā)展過程的主要因素，通氣量和孔口直徑對氣泡的發(fā)展過程有一定的影響。Zhang等[10]和劉濤濤等[4]通過湍流模型對通氣空泡流場特性進行研究，結(jié)果表明，在回射流與主流的相互作用下，反向旋轉(zhuǎn)的渦對導(dǎo)致了空泡的脫落。

為了形成覆蓋于航行體表面的通氣空泡，當(dāng)前工程上的主動通氣主要采用多孔排氣技術(shù)，但是通氣空泡在航行體表面不易融合，空泡形態(tài)不穩(wěn)定，龔瑞巖[9]在研究中發(fā)現(xiàn)，通氣孔形式下的通氣空泡受到通氣率、氣孔數(shù)量以及空泡外界流場等因素的影響，空泡的融合存在困難，未融合的通氣空泡表面凹凸不平存在大量褶皺，空泡尾部氣水混合嚴(yán)重，存在大量的空泡脫落及斷裂現(xiàn)象，不利于航行體的穩(wěn)定性，氣泡融合成為主動通氣技術(shù)的難點之一[11]。為解決空泡融合問題，改變航行體結(jié)構(gòu)的被動控制方法被大量應(yīng)用，胡少峰[12]和孫鐵志等[13]針對氣孔數(shù)量對通氣空泡融合的影響進行研究，發(fā)現(xiàn)氣孔數(shù)量的增加有利于空泡的融合，氣孔數(shù)量的增加實際上減小了氣孔間距，隨著氣孔間距的減小，氣泡軸向發(fā)展速度減緩，并開始向航行體周向發(fā)展并促進空泡的融合[14]。馬貴輝等[15]和時俊雯[16]通過CFD對雙排氣孔排氣的流動特性進行研究，發(fā)現(xiàn)由于前排氣孔的影響，氣體經(jīng)后排氣孔流出能較快地形成近似等壓區(qū)，有利于空泡的周向發(fā)展以及融合。但龔瑞巖[9]通過雙排氣孔排氣試驗發(fā)現(xiàn)，雖然雙排氣孔通氣促進了空泡的融合，但是通氣空泡表面不穩(wěn)定，尾部泄氣現(xiàn)象嚴(yán)重，不利于空泡的穩(wěn)定性。崔震宇[17]在航行體表面設(shè)置環(huán)形凹槽并在槽內(nèi)開孔通氣，研究表明該通氣方式能夠有效降低航行體表面波動，并促進空泡的融合，但是環(huán)形槽內(nèi)開孔又會造成航行體結(jié)構(gòu)強度的損失。

孫龍泉等[20]通過在航行體通氣孔后開槽的方式控制通氣空泡的融合,并對開槽控制空泡融合的機理進行了細(xì)致的分析，但對如何選擇合理的凹槽參數(shù)并未給出合理方案。本文通過有限體積法進行數(shù)值計算，采用VOF模型和動網(wǎng)格技術(shù)，分析了凹槽參數(shù)(槽寬和孔槽間距)對空泡融合控制效果的影響，給出最佳凹槽參數(shù)。

1 數(shù)值計算方法及驗證

1.1 幾何模型及邊界條件設(shè)置

數(shù)值模型示意如圖1所示。行進方向為航行體長度L方向(方向垂直于水平面)，零點設(shè)置在航行體尾部中點處，x為航行體截面到尾部的垂直距離。半球頭型的水下航行體直徑D=40 mm，總長度L=4.5D，凹槽寬度為d1，氣孔和凹槽的距離為d2，設(shè)置30個通氣孔均布于航行體頭部氣孔直徑d=0.01D。

圖1 幾何模型示意圖Fig.1 Geometric model

圖2為三維模型及計算域示意圖。計算域分為靜域、動域，靜域不參與航行體的運動，半徑為4.0D；通過滑移界面與靜域相連的為動域，半徑為1.5D，用相同的UDF定義動域和航行體的運動。計算域的上下表面設(shè)置為恒定壓力出口，側(cè)表面設(shè)置為壓力隨水深變化的壓力出口，由于本文計算不考慮航行體出水，所以出口處的回流均設(shè)置為水相，航行體表面設(shè)為無滑移、絕熱，氣孔邊界設(shè)置為恒定的質(zhì)量流量入口。由圖2可見，通過布置節(jié)點，實現(xiàn)孔內(nèi)網(wǎng)格到孔外網(wǎng)格由密到疏的過渡。為保證捕捉精細(xì)流場，設(shè)置半徑為1.0D的包裹域，給定10層邊界層，第一層網(wǎng)格高度設(shè)置為0.002D，邊界層內(nèi)網(wǎng)格增長率為1.2。經(jīng)過網(wǎng)格無關(guān)性驗證，同時考慮計算的準(zhǔn)確性和有效性，最終整體網(wǎng)格數(shù)量確定為260萬。

圖2 三維模型、計算域示意圖Fig.2 Three dimensional model and computational domain

將本文部分參數(shù)進行無量綱化處理，得到無量綱時間T，無量綱壓力P，通氣率Q，傅汝德數(shù)Fr，空泡長度Lp

(1)

式中，V為航行體的運動速度，p為航行體壁面壓力，Qbody為氣孔通氣量，lp為空泡長度。

1.2 數(shù)值計算方法及有效性驗證

計算采用有限體積法對控制微分方程進行求解，采用隱式VOF模型和Level Set方法進行相間界面的多相流計算，氣體設(shè)置為主相，水設(shè)置為副相(相函數(shù)α表征主相所占體積分?jǐn)?shù)，α=1時表示該區(qū)域只存在氣相，α=0時表示該區(qū)域只存在水相)。湍流模型采用RNGk-ε模型，對于雷諾數(shù)低、湍流發(fā)展不充分的近壁面區(qū)域，通過標(biāo)準(zhǔn)壁面函數(shù)進行近壁面處理；壓力速度耦合算法采用PISO算法。

為驗證算法的有效性，本文針對無槽多孔通氣進行仿真，并在哈爾濱工程大學(xué)減壓罐中進行排氣試驗，試驗裝置總體示意圖如圖3(a)所示。試驗系統(tǒng)包括3個主要部分：由空壓機和儲氣罐以及氣體流量控制柜組成的通氣控制系統(tǒng)，由高速攝像機組成的數(shù)據(jù)采集系統(tǒng)，由傳動控制柜和減壓罐組成的運動系統(tǒng)。試驗?zāi)Ｐ腿鐖D3(b)所示，模型由航行體本身、氣室、通氣管等部分組成，航行體直徑為D，長度為4.75D，氣孔直徑為0.02D。

(a)試驗操作系統(tǒng)

(b)實驗?zāi)Ｐ?/p>

在Fr=4.79，Q=0.025的工況進行試驗以及數(shù)值仿真，定義空泡長度Lp為從通氣孔到空泡尾部的距離，圖4給出了空泡無量綱長度的時歷曲線。通過對比發(fā)現(xiàn)，數(shù)值模擬結(jié)果和試驗結(jié)果吻合程度良好，可以認(rèn)為本文所采用的數(shù)值方法有效。

圖4 試驗和數(shù)值空泡無量綱長度對比Fig.4 Comparison of experimental and numerical cavity dimensionless length

2 凹槽對通氣空泡融合的控制作用

對孔后開槽的空泡流動進行研究，取槽寬d1=0.125D，孔槽間距d2=0.075D的典型工況進行分析(本文所有計算工況的傅汝德數(shù)Fr=4.75，通氣率Q=0.025)。圖5為空泡充分發(fā)展時的三維泡型；圖6為當(dāng)X/L=0.80，X/L=0.68和X/L=0.57時，這3個典型截面位置的氣水相體積分布;圖7為3個典型截面位置的渦量云圖。

渦量云圖中渦量值的大小反映該位置漩渦作用的強度和剪切作用的大小[17]。圖6和圖7分別為航行體截面的氣水相云圖和渦量云圖，可以看出，渦量云圖中數(shù)值較大的紅色區(qū)域大多位于氣水相云圖的純氣相區(qū)與純水相區(qū)域的交界位置，這說明在這一位置存在強度較大的漩渦作用或受到較大的流體剪切作用。在截面X/L=0.80位置，未開槽模型的渦量云圖呈現(xiàn)花瓣狀，渦量值較大的紅色區(qū)域位于辮狀空泡內(nèi)，因為在此區(qū)域剛從排氣孔排出的氣體沿軸線方向相對航行體速度較小，因此受到較強的水流黏性拖拽作用所致。對于開槽模型的渦量云圖渦量值較大的紅色區(qū)域周向聯(lián)通，因為此處通氣空泡周向融合充分，融合的空泡外表面受到流體黏性拖拽作用從而產(chǎn)生較大的剪切力所致。這說明通氣空泡由于同時受到浮力和空泡外流體黏性拖拽力的作用，在空泡外表面附近往往會產(chǎn)生較大的流體剪切作用，即渦量云圖中渦量值較大的位置實際反映的就是通氣空泡外表面或空泡內(nèi)純氣相區(qū)與泡內(nèi)水團的交界面。因此，本文主要以渦量云圖作為重要的通氣空泡融合判據(jù)，并結(jié)合氣水相云圖和流線圖對空泡融合進行分析。

如圖7所示，開槽模型渦量值較大的紅色區(qū)域聯(lián)通程度較好，尤其在截面位置X/L=0.68和0.58，相比紅色區(qū)域聯(lián)通較好的開槽模型，開槽模型的渦量分布更為分散，說明開槽模型空泡融合程度較好。

圖7 開槽模型和未開槽模型典型位置渦量圖Fig.7 Vorticity magnitude at typical positions of slotted model and unslotted model

如圖6所示，相比未開槽模型，開槽模型空泡的氣相整體分布均勻，空泡外緣形狀規(guī)律，空泡的融合效果較好，開槽改變了空泡的流場特性，泡內(nèi)流線對比如圖5(b)所示。通氣空泡流經(jīng)凹槽時，由于壁面的突然擴張，空泡與壁面發(fā)生分離并形成自由剪切層[19-20]，與凹槽發(fā)生碰撞后形成旋渦，流經(jīng)凹槽的氣體受到旋渦的卷吸作用進入槽中。來流的沖刷對槽內(nèi)氣體流出凹槽造成困難，導(dǎo)致槽內(nèi)氣體會向周向擴散，促進上游通氣空泡的融合，如圖5(a)所示。與無槽模型的空泡相比，有槽模型的空泡在凹槽附近已經(jīng)發(fā)生融合并向空泡下游擴張，空泡整體的融合程度較好；凹槽的存在會增加近壁面處的速度梯度，增大流體的動能，阻礙空泡尾部回射流的發(fā)展，降低回射流的速度[20]，如圖5(a)所示。與無槽模型的相比較，由于開槽模型的回射流強度低，泡尾無大量凸起且邊界形狀較為規(guī)則。

(a)空泡形態(tài)

(b)對稱面流線圖5 開槽模型和未開槽模型通氣空泡對比Fig.5 Comparison of bubble between slotted model and unslotted model

圖6 開槽模型和未開槽模型典型位置氣水相分布圖Fig.6 Gas-water phase distribution at typical positions of slotted model and unslotted model

空泡的融合程度和回射流的強度影響航行體表面的壓力波動，進而影響航行體的結(jié)構(gòu)穩(wěn)定性[18]，為簡化航行體表面的壓力波動描述，引入航行體周向壓力標(biāo)準(zhǔn)差S

(2)

分別計算截面位置X/L=0.81，0.70，0.58時，3個典型位置的兩種模型的壓力標(biāo)準(zhǔn)差S，所得結(jié)果如圖8所示。可以發(fā)現(xiàn)，隨著截面位置向空泡下游推進，未開槽模型的壓力標(biāo)準(zhǔn)差S逐漸增大，航行體表面的整體壓力波動較大，因為未開槽模型的空泡融合效果較差。如圖6所示，隨著截面位置向泡尾推進，泡內(nèi)水相占比逐漸增加，空泡融合程度較差；而開槽模型的穩(wěn)壓效果較好，隨著截面位置向下游推進，壓力標(biāo)準(zhǔn)差S逐漸降低，在截面位置X/L=0.68，空泡已經(jīng)發(fā)生了一定程度上的融合，氣相分布較為均勻，在截面位置X/L=0.58，雖然受到了尾部回射流的影響，泡內(nèi)存在水相，但空泡總體融合程度高。

圖8 開槽和未開槽模型典型截面位置壓力標(biāo)準(zhǔn)差S對比Fig.8 Comparison of standard deviation of pressure at typical positions of slotted and unslotted models

開槽模型3個截面位置的壓力標(biāo)準(zhǔn)差S均小于未開槽模型的壓力標(biāo)準(zhǔn)差S，說明凹槽可以降低航行體表面的壓力波動。

3 不同開槽方案下的通氣空泡融合控制分析

通過氣孔后開槽的被動控制方式能促進通氣空泡的融合，但凹槽尺寸和開槽位置是影響控制效果的關(guān)鍵因素，選擇合理的凹槽尺寸以及開槽位置不僅有利于空泡的融合，而且能夠降低開槽對航行體結(jié)構(gòu)強度的負(fù)影響。本節(jié)主要針對凹槽寬度d1和孔槽間距d2對通氣空泡融合的影響進行分析，并確定最佳的開槽方案。

3.1 凹槽寬度對通氣空泡融合的影響

槽寬是影響空泡融合控制效果的重要因素，但在航行體表面開槽勢必會對航行體的結(jié)構(gòu)產(chǎn)生負(fù)面影響，所以在保證空泡融合效果良好的前提下，盡量減小槽寬。為探究槽寬對通氣空泡融合的影響，分別取0.075D,0.100D,0.125D和0.150D這4種槽寬進行計算，孔槽間距d2=0.075D。

不同槽寬模型的渦量云圖如圖9所示。按前文渦量云圖的分析方式，各槽寬模型渦量分布存在明顯差別，特別是在截面位置X/L=0.80處，截面位置位于各開槽模型凹槽上緣處?？梢钥吹?，此處槽寬0.100D,0.125D和0.150D模型渦量云圖紅色區(qū)域基本聯(lián)通，說明空泡融合狀態(tài)較好，特別是槽寬0.125D效果最佳。在X/L=0.68和0.58位置，槽寬0.125D模型渦量值較大的紅色區(qū)域聯(lián)通效果最好說明該模型空泡融合最佳，槽寬0.150D模型較其他模型也有較好的改善。所以，槽寬0.125D對空泡融合的控制效果最佳。

圖9 不同槽寬模型的渦量云圖Fig.9 Vorticity magnitude at typical positions under different groove widths

氣水相分布如圖10所示，在截面位置X/L=0.80處，該位置在氣孔下游附近，氣體從氣孔中排出形成的空泡呈辮狀，槽寬0.150D,0.125D和0.100D模型的相鄰辮狀空泡受到凹槽的影響發(fā)生連通，而槽寬0.075模型的相鄰辮狀空泡的連通現(xiàn)象不明顯，如圖11(a)所示。槽寬0.150D,0.125D和0.100D模型的空泡在凹槽后發(fā)生了融

圖10 不同槽寬模型的氣水相體積分布Fig.10 Gas-water phase distribution at typical positions under different groove widths

(a)空泡形態(tài)

(b)對稱面流線圖11 不同槽寬下通氣空泡對比(從上到下槽寬依次為0.150D，0.125D，0.100D，0.075D)Fig.11 Comparison of bub ble under different groove widths

合，而槽寬0.075D模型的空泡未發(fā)生融合，航行體壁面裸露在外，各個辮狀空泡相對較為獨立。在截面位置X/L=0.68和0.58處，凹槽對空泡融合的促進作用明顯，各個模型的空泡均已經(jīng)融合，但觀察圖10的氣水相分布云圖可以發(fā)現(xiàn)，相比于槽寬0.100D和0.075D模型，槽寬0.150D和0.125D模型的氣相分布更均勻，泡型更穩(wěn)定，如圖11(a)所示，槽寬0.100D和0.075D模型的空泡中部存在凸起，而寬0.150D和0.125D模型的空泡中部形狀更為規(guī)則。

在截面位置X/L=0.58處，該位置距離泡尾最近，空泡受到尾部回射流的影響，融合的空泡中混入水相，如圖10所示。但比于槽寬0.100D和0.075D模型，槽寬0.150D和0.125D模型的水相占比更小，泡型更加平穩(wěn)。

閆宸鋒[19]在研究中發(fā)現(xiàn)，流體順?biāo)俣确较蛴龅綑M截面突然擴張時，會與壁面分離并在分離區(qū)產(chǎn)生旋渦，旋渦經(jīng)歷生成、發(fā)展、遷移和破碎4個階段。如圖11(b)所示，槽寬0.100D和0.075D模型由于槽寬太小，產(chǎn)生的旋渦在發(fā)展之前就會接觸到凹槽末端的壁面而發(fā)生破碎，而槽寬0.150D和0.125D模型在凹槽產(chǎn)生的旋渦有一定程度上的發(fā)展，旋渦狀流動特征更明顯，這樣凹槽的卷吸作用明顯，促進了通氣空泡的軸向分離，對通氣空泡的軸向發(fā)展產(chǎn)生阻礙，槽外流經(jīng)的氣體更容易匯入槽內(nèi)并向周向發(fā)展，相鄰辮狀空泡更易周向發(fā)展融合，空泡的融合效果更好。

不同槽寬下經(jīng)典位置的壓力標(biāo)準(zhǔn)差S如圖12所示。隨著截面位置向空泡尾部推進，壓力標(biāo)準(zhǔn)差S逐漸降低，在截面位置X/L=0.80處，雖然辮狀空泡受到凹槽的影響在航行體周向產(chǎn)生了一定程度上的連通，但是空泡融合程度低，受到流體擾動嚴(yán)重，導(dǎo)致壓力標(biāo)準(zhǔn)差S較大；在截面位置X/L=0.68處，空泡已經(jīng)發(fā)生了融合，但是融合程度較低，如圖10所示，該位置的空泡外緣仍然呈辮狀，導(dǎo)致該位置的壓力標(biāo)準(zhǔn)差S只相對于截面位置X/L=0.80處有所降低；在截面位置X/L=0.58處，該位置空泡的融合程度較好，雖然受到了回射流的影響，泡內(nèi)存在一定的水相，但是總體上穩(wěn)壓效果較好，壓力標(biāo)準(zhǔn)差S最低。

圖12 不同槽寬模型典型截面位置壓力標(biāo)準(zhǔn)差S對比Fig.12 Comparison of standard deviation of pressure at typical positions under different groove widths

各個截面位置的壓力標(biāo)準(zhǔn)差S受到空泡融合程度的影響在穩(wěn)壓效果上存在差異，相比于槽寬0.100D和0.075D模型，槽寬0.150D和0.125D模型對航行體表面穩(wěn)壓效果較好而且二者的穩(wěn)壓效果相近，同時通過對T=20時的空泡形態(tài)和氣水相分布的分析發(fā)現(xiàn)，二者對通氣空泡融合的控制效果也近似，但槽寬的增加會降低航行體的結(jié)構(gòu)穩(wěn)定性，取槽寬0.125D為最佳。

3.2 孔槽間距對通氣空泡融合的影響

孔槽間距也是影響空泡融合控制效果的關(guān)鍵因素，合理選擇開槽位置可以有效避開水下航行體的關(guān)鍵部位，從而降低開槽對航行體局部強度的負(fù)影響。為探究孔槽間距對空泡融合的影響，分別取0.025D，0.075D，0.150D，0.300D這4種不同的孔槽間距，槽寬d1=0.125D。

如圖13所示，無論在任何截面位置，各模型渦量分布存在明顯差別，特別是在截面位置X/L=0.68和0.57處?？梢钥吹?，只有孔槽間距0.075D模型渦量云圖紅色區(qū)域聯(lián)通，其他模型的渦量分布散亂，說明孔槽間距0.075D對空泡融合的促進效果較好。因此，孔槽間距0.075D對空泡融合的控制效果最佳。

圖13 不同孔槽間距模型的渦量云圖Fig.13 Vorticity magnitude at at typical positions under different hole-groove spacing

氣水相分布如圖14所示?？撞坶g距0.025D模型3個典型位置的氣水混合現(xiàn)象嚴(yán)重而且分布散亂，如圖15(a)所示，通氣空泡整體的融合程度較差，空泡尾部存在嚴(yán)重的空泡脫落現(xiàn)象，這是因為從氣孔排出的氣體需要一定發(fā)展距離，而孔槽間距0.025D模型的氣孔和凹槽的距離過近，如圖15(b)所示，流經(jīng)凹槽的氣體未在槽內(nèi)形成較強的漩渦將氣體卷入凹槽內(nèi)，相比于其他模型，孔槽間距0.075D模型的槽內(nèi)存在大量水相。

觀察圖14可知，相比于孔槽間距0.025D模型，孔槽間距0.300D，0.150D和0.075D模型對空泡融合的控制效果較好，不同截面位置的空泡都已經(jīng)發(fā)生了不同程度的融合，但孔槽間距0.300D模型空泡在周向分布不均勻，在圖15(a)中表現(xiàn)為該模型的空泡表面存在凸起，泡尾邊界參差不齊，這是因為孔槽間距過大，在回射流的影響下凹槽位置附近存在大小不一的漩渦，如圖15(b)所示。從而影響了槽內(nèi)漩渦對氣體的卷吸作用，造成空泡周向分布不均的現(xiàn)象。

圖14 不同孔槽間距模型的氣水相體積分布Fig.14 Gas-water phase distribution at typical positions under different hole-groove spacing

(a)空泡形態(tài)

(b)對稱面流線圖15 不同孔槽間距下通氣空泡對比(從上到下孔槽間距依次為0.300D，0.150D，0.075D，0.025D)Fig.15 Comparison of bubble under different hole-groove spacing

不同孔槽間距下經(jīng)典位置的壓力標(biāo)準(zhǔn)差S如圖16所示，相比于槽寬對壓力標(biāo)準(zhǔn)差S的影響，孔槽間距改變了壓力標(biāo)準(zhǔn)差S沿空泡下游推進的規(guī)律?？撞坶g距0.025D和0.300D模型壓力標(biāo)準(zhǔn)差S變化規(guī)律和相同，在截面位置X/L=0.80處各個辮狀空泡相對獨立，氣水相混合嚴(yán)重，壓力標(biāo)準(zhǔn)差S較大；在截面位置X/L=0.68處，空泡發(fā)生融合，壓力標(biāo)準(zhǔn)差S較截面位置X/L=0.80處有所下降；在截面位置X/L=0.57處，該位置距離空泡尾部較近，受到回射流的影響程度大，空泡內(nèi)部存在大量水相，故壓力標(biāo)準(zhǔn)差S最大。但相比孔槽間距0.025D模型，孔槽間距0.300D模型的氣相分布更加均勻，故整體的壓力標(biāo)準(zhǔn)S差較小。

圖16 不同孔槽間距模型典型截面位置壓力標(biāo)準(zhǔn)差S對比 Fig.16 Comparison of standard deviation of pressure at typical positions under different hole-groove spacing

相比于空泡融合程度最好的孔槽間距0.075D模型孔，孔槽間距0.150D模型的凹槽參數(shù)使空泡完全融合區(qū)向空泡下游推移，截面位置X/L=0.80處，辮狀空泡在航行體周向發(fā)生了一定程度的連通，而在截面位置X/L=0.68處，該位置位于空泡融合的上邊界附近，氣水混合復(fù)雜，故該位置處壓力標(biāo)準(zhǔn)差S有所提高；在截面位置X/L=0.57處，氣相分布均勻，空泡的融合程度較好，故壓力標(biāo)準(zhǔn)差S降低。

通過對比發(fā)現(xiàn)孔槽間距0.075D模型的穩(wěn)壓效果和對空泡融合的促進作用最好，所以孔槽間距選取0.075D最佳。

4 結(jié)論

本文針對凹槽對通氣空泡融合的控制效果，開展了計算分析，得出以下結(jié)論：

1)凹槽改變了航行體表面的流場特性，在凹槽內(nèi)部產(chǎn)生旋渦，將流經(jīng)凹槽的氣體卷吸到槽內(nèi)，槽內(nèi)氣體向航行體周向擴散促進空泡的融合，能夠降低航行體表面的壓力波動。

2)槽寬d1取0.125D時，對空泡融合的促進以及穩(wěn)壓效果最佳。槽寬過小時，槽內(nèi)旋渦沒有足夠的空間發(fā)展就已經(jīng)發(fā)生破碎，對空泡融合的促進作用不明顯；槽寬過大時，相比于槽寬0.125D，對空泡的融合促進作用以及穩(wěn)壓效果提升不明顯，而且不利于航行體的結(jié)構(gòu)穩(wěn)定性。

3)孔槽間距d2取0.075D時，對空泡融合的促進以及穩(wěn)壓效果最佳?？撞坶g距過小時，無法誘導(dǎo)空泡在槽內(nèi)產(chǎn)生旋渦并促進空泡融合；孔槽間距過大時，在回射流的影響下凹槽位置附近存在大小不一的漩渦，影響了凹槽對氣體的卷吸作用，降低了凹槽對空泡融合的促進作用。