（上海理工大學(xué)能源與動(dòng)力工程學(xué)院,上海 200093）

（上海理工大學(xué)能源與動(dòng)力工程學(xué)院,上海 200093）

為了加深對(duì)水射流沖擊平板問(wèn)題的理解，將該問(wèn)題簡(jiǎn)化為二維模型,并對(duì)此二維模型進(jìn)行了理論分析和數(shù)值模擬。理論分析采用理想流體模型，并在忽略重力條件下運(yùn)用流體連續(xù)性方程、動(dòng)量方程和能量方程對(duì)該問(wèn)題進(jìn)行了研究。此外，采用fluent軟件中的VOF模型對(duì)該問(wèn)題進(jìn)行數(shù)值模擬。最后將兩種方法獲得的平板受力和流量分配情況進(jìn)行對(duì)比，兩者吻合良好。

射流沖擊；理論分析；數(shù)值模擬；VOF模型

0 引言

射流，是指流體從管口、孔口、狹縫中射出，并同周?chē)黧w摻混的一股流體流動(dòng)，是一種廣泛存在于自然界和工業(yè)生產(chǎn)中的現(xiàn)象，如噴泉、蒸汽泵、噴氣式飛機(jī)等都是射流現(xiàn)象的典型代表。

射流運(yùn)動(dòng)具有一定的能量，在工業(yè)上常利用射流沖擊材料實(shí)現(xiàn)清洗、剝層、切割等工序，而此時(shí)材料受射流沖擊力的大小以及射流流量分配為重要的參數(shù)。本文通過(guò)理論分析與數(shù)值模擬兩種方法，分別對(duì)平板所受水流沖擊力以及平板上的流量分配進(jìn)行了研究。

1 物理描述

本文考慮二維情況下的水射流沖擊平板問(wèn)題，并對(duì)其進(jìn)行了分析。二維模型如圖1所示,一股速度為 ，體積流量為 的水射流從管道入射，沖擊在一傾斜角為 的平板上，其后在平板上分為兩股流體，體積流量分別為 和 ，根據(jù)理論分析可得到平板所受沖擊力以及平板上流量的分配情況。

圖1 水射流沖擊平板問(wèn)題示意圖

2 理論計(jì)算

真實(shí)流體所具有的粘性會(huì)對(duì)理論分析帶來(lái)極大的困難，簡(jiǎn)便起見(jiàn)，這里考慮入射流體為處于環(huán)境壓力內(nèi)的定常不可壓縮理想流體，并忽略流體所受的重力，此時(shí)建立如圖1所示的坐標(biāo)系，設(shè)平板對(duì)水流的沖擊力方向如圖1所示。

由動(dòng)量方程得：

在定常流動(dòng)的情況下，上式右邊第一項(xiàng)為零。因此其分量式為：

由圖1所示可知其具體表達(dá)式分別為：

流量分配情況分析如下：

由連續(xù)性方程可得：Q0=Q1+Q2

射流處于環(huán)境壓力，即控制體各有效截面壓力均相等為環(huán)境壓力，

由能量方程得：

即：v0=v1=v2

帶入x方向動(dòng)量方程；

化簡(jiǎn)得：Q1-Q2=Q0cosθ

最終可得：

3 數(shù)值計(jì)算

3.1 數(shù)值計(jì)算方法

考慮空氣中水射流沖擊平板問(wèn)題，水流的流動(dòng)必然引起周?chē)諝獾牧鲃?dòng)，故為氣液兩相流問(wèn)題，計(jì)算中采用“流體體積”（Volume of Fluid-VOF）方法求解。

VOF方法作為一種解決兩相流問(wèn)題的方法，是通過(guò)體積函數(shù)來(lái)追蹤相間界面，具有物理概念清晰、容易實(shí)現(xiàn)、計(jì)算量小以及精度較高的優(yōu)點(diǎn)。

追蹤各相流體的體積分?jǐn)?shù)可以通過(guò)連續(xù)性方程來(lái)求解，例如第 相流體的體積分?jǐn)?shù)可以通過(guò)如下的連續(xù)性方程求解：

式中：ρk是第k相流體的密度；mzk是第k相流體流入第Z相流體中的質(zhì)量；mkz是第Z相流體流入第k相流體中的質(zhì)量；Sk為源項(xiàng)，如果多相流體間沒(méi)有特殊的質(zhì)量輸運(yùn)，該項(xiàng)一般取0。

3.2 計(jì)算區(qū)域及邊界條件

建立如圖2所示流動(dòng)區(qū)域進(jìn)行模擬，其中BC，AD為水流入射管道邊界，設(shè)為固壁，經(jīng) AB 邊界流入的水流沖擊 EFGH矩形平板，流動(dòng)區(qū)域邊界BJ、AK、IL、IJ、KL設(shè)為壓力出口。初始時(shí)流體區(qū)域均為靜止的空氣。

圖2 流動(dòng)區(qū)域圖

3.3 初始條件

AB為速度進(jìn)口，入射水流速 度v=1m/s， 密 度ρ=1000kg/m3，AB段長(zhǎng)度為10mm，入射平板傾角θ為45°，即有效射流面積A=0.01×sina45°=0.0071m2。平板厚10mm，長(zhǎng)300mm。由以上數(shù)據(jù)計(jì)算得進(jìn)口體積流量Q0為0.01414m3。

3.4 網(wǎng)格劃分

采用不同間距對(duì)流動(dòng)區(qū)域網(wǎng)格進(jìn)行了四次劃分，并以平板受力大小為依據(jù)進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證，驗(yàn)證情況如表1：

表1 網(wǎng)格獨(dú)立性驗(yàn)證表

驗(yàn)證結(jié)果如上表所示，故選用網(wǎng)格數(shù)為31917個(gè)。對(duì)流體區(qū)域網(wǎng)格劃分具體情況如圖3所示。

圖3 網(wǎng)格劃分圖

3.5 計(jì)算結(jié)果

用CFD軟件對(duì)模型進(jìn)行仿真計(jì)算。圖4為射流流場(chǎng)中液相體積分?jǐn)?shù)分布，圖5為流場(chǎng)密度分布圖，圖6為流場(chǎng)中的流線(xiàn)示意圖，圖7為流體速度矢量圖。

圖4中紅色部分為入射的水流，由圖4可見(jiàn)，水流經(jīng)管道流出后，入射到平板上，在平板處分為體積流量不同的兩股流體。圖5的密度分布與圖4的水體積分布圖結(jié)果圖形相似，可以較好吻合。由圖6和圖7可以看出，水射流的流動(dòng)可以引起周?chē)諝獾牧鲃?dòng)，在平板射流分為兩股處空氣的流場(chǎng)會(huì)形成漩渦。

圖4 水體積分?jǐn)?shù)圖

圖5 密度分布圖

圖6 流線(xiàn)示意圖

圖7 速度矢量圖

以垂直于平板方向?yàn)閤軸方向，CFD仿真分析結(jié)果為作用在平板上的力為：

流量在平板處分配情況為：

4 結(jié)果對(duì)比及分析

數(shù)值模擬計(jì)算結(jié)果為：

理論計(jì)算結(jié)果為：

將v0=1m/s，Qv=0.01414m3/s，θ=45°，ρ=1000kg/m3代入理論計(jì)算得出公式：

表2 計(jì)算結(jié)果對(duì)比表

由上表可知理論計(jì)算結(jié)果與數(shù)值計(jì)算結(jié)果相差很小。

分析：

（1）在此水射流沖擊平板的二維模型中，數(shù)值模擬結(jié)果與理論計(jì)算結(jié)果可以較好吻合。

（2）由計(jì)算可知數(shù)值模擬值與理論值相接近，可知VOF方法對(duì)水射流沖擊平板的二維模型有較好的模擬。

（3）在理論計(jì)算的過(guò)程中，由于不考慮摩擦（即視為理想流體）和忽略了重力，故肯定與實(shí)際情況存在一定偏差。

（4）數(shù)值模擬結(jié)果可能存在誤差，分析造成模擬值有可能存在誤差的原因有：流體邊界值難以預(yù)測(cè)，邊界條件不一定準(zhǔn)確，網(wǎng)格劃分過(guò)程中存在截?cái)嗾`差，選用的模型中引入一些假設(shè)等。

（5）理論分析與數(shù)值模擬相比，更加簡(jiǎn)便快捷，但是計(jì)算結(jié)果較為簡(jiǎn)單，數(shù)值模擬可以獲得更多的流場(chǎng)信息。

5 總結(jié)

在學(xué)習(xí)過(guò)程中要善于思考，把學(xué)到的知識(shí)綜合運(yùn)用起來(lái)。本文將流體力學(xué)知識(shí)和計(jì)算流體知識(shí)相結(jié)合，互相驗(yàn)證，有利于對(duì)兩種知識(shí)更進(jìn)一步的理解和掌握。

各種數(shù)據(jù)顯示數(shù)值模擬的結(jié)果與實(shí)際情況存在一定的誤差，但對(duì)于本文所研究的二維問(wèn)題而言，數(shù)值模擬與理論分析誤差比較小，都可以較好的符合實(shí)際情況，大大了方便研究分析，雖然數(shù)值模擬的方法沒(méi)有理論分析快捷方便，但是可以顯示流線(xiàn)分布情況，速度矢量分布情況等大量的流場(chǎng)信息，便于更加全面、直觀的掌握流動(dòng)細(xì)節(jié)，是工程應(yīng)用中一種很方便的工具。

參考資料：

[1]倪玲英.工程流體力學(xué)[M].北京：中國(guó)石油大學(xué)出版社,2012.

[2]陳慶光，徐忠，張永建.軸對(duì)稱(chēng)湍流沖擊射流場(chǎng)的數(shù)值預(yù)測(cè)[J].動(dòng)力工程,2002,22（06）：2015-2019.

[3]趙永志，顧兆林，郁永章等.盆池渦渦動(dòng)過(guò)程數(shù)值研究[J].水利學(xué)報(bào),2002（12）：1-6.

[4]王福軍.計(jì)算流體動(dòng)力學(xué)分析——CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社,2004.

水射流沖擊平板問(wèn)題分析

白于于,楊 帆

白于于（1993—），女，山西呂梁人，本科，從事cfd數(shù)值模擬研究。