国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

空泡水筒流動(dòng)品質(zhì)數(shù)值模擬

2019-01-07 08:01:36恒,
關(guān)鍵詞:附面層軸流泵空泡

劉 恒, 伍 銳

(上海船舶運(yùn)輸科學(xué)研究所, 上海 200135)

0 引 言

空泡水筒是研究船舶和高速巡航體的水動(dòng)力學(xué)性能的重要試驗(yàn)設(shè)備之一,為推進(jìn)器水動(dòng)力性能測試[1-2]和超空泡減阻[3]等提供試驗(yàn)結(jié)果??张菟矐?yīng)能提供較佳的流場品質(zhì),主要指空泡水筒工作段流場參數(shù)在時(shí)間上和空間上的均勻性、紊流強(qiáng)度及附面層厚度均能滿足試驗(yàn)要求。以螺旋槳空泡性能試驗(yàn)為例,根據(jù)相似理論,物理模型試驗(yàn)不僅要求推進(jìn)器幾何相似,而且要求其必須滿足雷諾數(shù)、轉(zhuǎn)速空化數(shù)和推力系數(shù)相似,而在實(shí)際試驗(yàn)中這是不可能完全實(shí)現(xiàn)的。近年來,隨著計(jì)算機(jī)技術(shù)和計(jì)算流體動(dòng)力學(xué)(Computational Fluid Dynamics, CFD)技術(shù)的快速發(fā)展,采用CFD能較為真實(shí)地模擬試驗(yàn)段流場的特性[4]。因此,運(yùn)用CFD方法對(duì)整個(gè)空泡水筒進(jìn)行數(shù)值模擬,能給出整個(gè)空泡水筒的流場分布特征。

本文通過求解非定常不可壓縮RANS方程,對(duì)某空泡水筒流場進(jìn)行數(shù)值模擬,探究空泡水筒不同位置處的流場分布特征。對(duì)試驗(yàn)段流場的特性進(jìn)行分析,重點(diǎn)分析工作段不同斷面處的流速分布及流場參數(shù)在空間上和時(shí)間上的均勻性、湍流強(qiáng)度和附面層厚度等流場特征。采用皮托耙對(duì)工作段不同斷面進(jìn)行流場測量,以驗(yàn)證數(shù)值結(jié)果的準(zhǔn)確性。

1 幾何模型

水筒工作段長2.6 m,橫截面呈方形(帶圓角),尺寸為0.6 m×0.6 m(見圖1中的區(qū)域③)。收縮段長2.3 m,擴(kuò)散段長4.1 m(分別見圖1中的區(qū)域⑤和區(qū)域⑥)。蜂窩器長400 mm,橫截面為正方形,導(dǎo)片厚度為1 mm,位于收縮段前端(見圖1中的區(qū)域②),其三維模型見圖2。導(dǎo)流片位于空泡水筒4個(gè)拐角處(見圖1中的區(qū)域①)。軸流泵處于空泡水筒最底端(見圖1中的區(qū)域④),以免槳葉旋轉(zhuǎn)時(shí)產(chǎn)生空泡。圖1中,坐標(biāo)原點(diǎn)位于試驗(yàn)段進(jìn)口軸中心線處。

圖1 空泡水筒三維結(jié)構(gòu)示意

b) 軸流泵

圖2導(dǎo)流片、蜂窩器、軸系和軸流泵三維模型

2 網(wǎng)格劃分

本文采用Star ccm+軟件進(jìn)行網(wǎng)格劃分,該軟件具有較強(qiáng)的網(wǎng)格劃分能力,可為獲取高質(zhì)量網(wǎng)格提供保證。網(wǎng)格離散采用六面體核心網(wǎng)格(Trim)類型,Trim網(wǎng)格模型內(nèi)置有一個(gè)用于控制表面模板單元格尺寸的上限,該上限可避免幾何模型嚴(yán)重失真和生成質(zhì)量差的網(wǎng)格。鑒于本文計(jì)算域的尺寸,且在區(qū)域①~區(qū)域⑤進(jìn)行網(wǎng)格加密,整個(gè)計(jì)算域的總網(wǎng)格數(shù)約為920萬個(gè)(見圖3)。水筒壁、軸系、導(dǎo)流片和軸流泵均采用邊界層形式的計(jì)算網(wǎng)格,以提高邊界層流動(dòng)的模擬精度。壁面為普通壁面(wall),采用標(biāo)準(zhǔn)壁面函數(shù);在本文的計(jì)算工況下,y+均在30~300范圍內(nèi),滿足壁面函數(shù)的要求。由于蜂窩器導(dǎo)流片厚度較小,其壁面邊界設(shè)置為滑移壁面邊界。鑒于軸流泵轉(zhuǎn)速直接控制水筒工作段流速,單獨(dú)為軸流泵劃分一個(gè)旋轉(zhuǎn)的區(qū)域。采用滑移網(wǎng)格法進(jìn)行求解,轉(zhuǎn)速設(shè)定為106 r/min。為保證軸流泵內(nèi)外區(qū)域速度矢量的連續(xù)性,在2個(gè)子區(qū)域的交界面處進(jìn)行連續(xù)性設(shè)置。

3 三維流場數(shù)學(xué)模型

3.1 控制方程

圖3 空泡水筒網(wǎng)格示意(y=0截面)

空泡水筒運(yùn)行時(shí)水流為三維流動(dòng),假定流體不可壓縮,其控制方程包括連續(xù)性方程和動(dòng)量方程,有

(1)

雷諾平均N-S方程(RANS)為

(2)

式(1)和式(2)中:ui和uj為速度分量的時(shí)均值;P為壓力時(shí)均值;ρ為水的密度;μ為分子黏性系數(shù);gi為重力加速度分量。

3.2 湍流模型

在流體力學(xué)中,通常采用由歐拉方程延伸出的k-ε雙方程對(duì)湍流流場進(jìn)行描述。由于本文研究的空泡水筒流場中既有管道內(nèi)的流動(dòng),又有分離流動(dòng),因此在建立仿真方案時(shí)選取Realizablek-epsilon湍流模型[6]求解湍流黏性系數(shù)。該模型湍流黏度計(jì)算式中含有與旋轉(zhuǎn)和曲率相關(guān)的內(nèi)容,因此適于計(jì)算旋轉(zhuǎn)均勻的剪切流、管道內(nèi)流動(dòng)、邊界層流動(dòng)和帶有分離的流動(dòng)。

Realizablek-epsilon湍流模型方程如下:

1) 湍動(dòng)能k滿足

(3)

2) 湍流耗散率ε滿足

(4)

式(3)和式(4)中:xi和xj為坐標(biāo)分量;σk和σε分別為湍動(dòng)能k和耗散率ε對(duì)應(yīng)的Prantal數(shù);Gk為平均速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng);C1和C2為經(jīng)驗(yàn)常數(shù);μt為湍流黏性系數(shù)。湍流黏性系數(shù)μt的計(jì)算式為

(5)

(6)

4 皮托耙流速測量

采用梳狀皮托耙測量工作段不同斷面的流速場,以驗(yàn)證數(shù)值結(jié)果的準(zhǔn)確性。皮托耙共包含7個(gè)皮托管,各皮托管與軸中心線的距離分別為r1=49 mm、r2=63 mm、r3=78 mm、r4=95 mm、r5=114 mm、r6=135 mm、r7=152 mm,布置方式見圖4。皮托管是用來測量水流總壓與靜壓之差以確定水流速度的一種裝置,其結(jié)構(gòu)簡單且具有較高的測量精度,但會(huì)對(duì)流場產(chǎn)生擾動(dòng),無法測量水流紊動(dòng)強(qiáng)度。試驗(yàn)時(shí)設(shè)定軸流泵轉(zhuǎn)速為106 r/min,此時(shí)工作段進(jìn)口處(x=0)斷面的平均流速為4 m/s,皮托耙以10°為間隔,旋轉(zhuǎn)360°,共測量252個(gè)點(diǎn)在8 s內(nèi)的平均速度,進(jìn)而求解斷面x=0.5 m和x=1.0 m處的速度均勻性。

圖4 空泡水筒中測量伴流用的梳狀皮托耙

5 數(shù)值計(jì)算結(jié)果

當(dāng)軸流泵以106 r/min的轉(zhuǎn)速旋轉(zhuǎn)時(shí),空泡水筒試驗(yàn)段進(jìn)口處的平均流速為4 m/s,中心截面(y=0)處的速度場分布見圖5。由圖5可知,擴(kuò)散段下游和導(dǎo)流片①D下游均產(chǎn)生了流動(dòng)分離。水流經(jīng)過導(dǎo)流片①C之后,其速度幅值沿z軸的分布較不均勻,此時(shí)軸流泵工作于非均勻流場中,但軸流泵處于水筒底部,壓力較大,葉片未發(fā)生空化。軸流泵下游存在旋轉(zhuǎn)流,速度分布較不均勻,導(dǎo)流片①B附近產(chǎn)生劇烈的流動(dòng)分離,湍動(dòng)能較大。

圖5 空泡水筒速度分布(y=0截面)

為使流經(jīng)導(dǎo)流片①A的水流能均勻地沿水平方向流動(dòng),在①A與收縮段之間設(shè)置整流段,并在整流段安置蜂窩器,可對(duì)水流起導(dǎo)向作用,并能降低水流的橫向湍流強(qiáng)度。圖6給出蜂窩器附近不同斷面的湍動(dòng)能分布。從圖6中可看出,現(xiàn)有的整流蜂窩器能起到有效的整流作用,水流經(jīng)過蜂窩器之后湍動(dòng)能得到有效降低。

a) 蜂窩器上游0.2 m斷面

b) 蜂窩器下游0.2 m斷面

5.1 水筒工作段流場分布

水流的黏性工作段筒壁附近會(huì)形成具有較大速度梯度的邊界層,為能較為準(zhǔn)確地捕捉該邊界層的厚度變化趨勢,對(duì)工作段的邊界層網(wǎng)格進(jìn)行加密,第一層網(wǎng)格的厚度為0.5 mm,共設(shè)置15層。圖7給出各網(wǎng)格節(jié)點(diǎn)位置處速度的無量綱分布,工作段進(jìn)口處軸向平均速度u*為4 m/s,工作段斷面寬度D為600 mm。由圖7可知,第一個(gè)網(wǎng)格可捕捉到平均速度的55%左右,其y+處于 30~55 范圍內(nèi),可準(zhǔn)確模擬筒壁附近的附面層厚度。工作段進(jìn)口至出口的附面層厚度發(fā)展分布見圖8中的陰影部分,附面層厚度沿水流方向逐漸增大,出口處附面層的厚度約為50 mm,在進(jìn)行船舶推進(jìn)器水動(dòng)力性能試驗(yàn)時(shí),物理模型位于0.5~1.0 m,該范圍內(nèi)附面層的厚度約為25 mm。圖9為不同位置處的平均瞬時(shí)速度。從圖9中可看出,速度有微小波動(dòng),波動(dòng)頻率與軸流泵軸頻相同,均為1.77 Hz。軸流泵推力沿時(shí)間歷程變化見圖10。

圖7 試驗(yàn)段速度分布

圖8 不同斷面處速度分布(y=0截面)

a) x=0.5 m處斷面平均瞬時(shí)流速

b) 坐標(biāo)點(diǎn)(0.5 m, 0.0 m, 0.3 m)處瞬時(shí)流速

圖10 軸流泵推力沿時(shí)間歷程變化

5.2 工作段軸向速度均勻度

速度均勻度定義為

(7)

表1 數(shù)值模擬與試驗(yàn)速度均勻度結(jié)果對(duì)比

5.3 工作段湍流強(qiáng)度分布

在數(shù)值計(jì)算過程中,共監(jiān)測11個(gè)軸中心位置處的湍動(dòng)能,試驗(yàn)段軸中心線位置的湍流強(qiáng)度分布見圖11。湍流強(qiáng)度定義為

(8)

式(8)中:k為湍動(dòng)能。由圖11可知,整個(gè)工作段內(nèi)軸向湍流強(qiáng)度?1%,湍流強(qiáng)度沿水流方向逐漸減弱并趨于定值0.26%。

圖11 試驗(yàn)段湍流強(qiáng)度分布

6 結(jié) 語

通過求解非定常不可壓縮RANS方程,對(duì)上海船舶運(yùn)輸科學(xué)研究所的空泡水筒流場進(jìn)行數(shù)值模擬,探究空泡水筒不同位置處的流場特性,重點(diǎn)分析工作段內(nèi)的流場分布,得到以下結(jié)論:

1) 擴(kuò)散段下游、導(dǎo)流片①D下游和導(dǎo)流片①B附近產(chǎn)生流動(dòng)分離,整流蜂窩器能有效發(fā)揮整流作用,湍流強(qiáng)度減弱;

2) 工作段附面層厚度沿水流方向逐漸增大,當(dāng)截面位置處于0.5~1.0 m時(shí),附面層厚度約為25 mm,在出口處達(dá)到最大值50 mm;

3) 工作段軸向速度均勻性較佳,數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果的吻合度良好,驗(yàn)證了數(shù)值模擬結(jié)果的準(zhǔn)確性;

4) 工作段軸向湍流強(qiáng)度?1%,且沿水流方向逐漸減弱并趨于定值0.26%。

猜你喜歡
附面層軸流泵空泡
潛水軸流泵運(yùn)行故障分析與排除研究
基于網(wǎng)格框架的非結(jié)構(gòu)附面層網(wǎng)格生成技術(shù)
潛水軸流泵電機(jī)運(yùn)行工況的特點(diǎn)及可靠性探討
水下航行體雙空泡相互作用數(shù)值模擬研究
基于數(shù)值模擬的流場附面層邊緣識(shí)別方法
基于LPV的超空泡航行體H∞抗飽和控制
濃縮軸流泵干氣密封改造
基于CFD的對(duì)轉(zhuǎn)槳無空泡噪聲的仿真預(yù)報(bào)
船海工程(2015年4期)2016-01-05 15:53:28
超聲壓氣機(jī)葉柵流場的數(shù)值模擬與試驗(yàn)驗(yàn)證
SPH在水下高速物體空泡發(fā)展模擬中的應(yīng)用
404 Not Found

404 Not Found


nginx
桃园市| 泸水县| 新巴尔虎右旗| 仁布县| 洪泽县| 江安县| 濮阳县| 项城市| 潍坊市| 宝山区| 尼玛县| 信阳市| 东城区| 大余县| 台北县| 南阳市| 壤塘县| 亚东县| 阿克苏市| 博白县| 永胜县| 遂宁市| 石嘴山市| 杭锦后旗| 额济纳旗| 芮城县| 彰化县| 安顺市| 延津县| 临汾市| 钦州市| 枞阳县| 玉门市| 南丹县| 淮北市| 四子王旗| 年辖:市辖区| 利辛县| 随州市| 洞口县| 中江县|