[發(fā)明專利]一種基于數(shù)值模擬的孔板與洞塞劃分方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201711061447.X
• 申請(qǐng)日: 2017-11-01
• 公開(kāi)（公告）號(hào): CN107895070A
• 公開(kāi)（公告）日: 2018-04-10
• 發(fā)明（設(shè)計(jì)）人: 艾萬(wàn)政;王偉軍
• 申請(qǐng)（專利權(quán)）人: 浙江海洋大學(xué)
• 主分類號(hào): G06F17/50
• 分類號(hào): G06F17/50
• 代理公司: 杭州杭誠(chéng)專利事務(wù)所有限公司33109
• 代理人: 尉偉敏

1.一種基于數(shù)值模擬的孔板與洞塞劃分方法，其特征是，一種基于數(shù)值模擬的孔板與洞塞劃分方法具體采用如下步驟：

步驟一、選定空腔長(zhǎng)度L為考察量，空腔長(zhǎng)度L小于消能孔厚度T時(shí)為洞塞，空腔長(zhǎng)度L大于消能孔厚度T時(shí)為孔板；

步驟二、確定影響空腔長(zhǎng)度L的因素，包括泄洪洞內(nèi)水流的平均流速u，水流密度ρ，水流粘度μ，泄洪洞直徑D，消能孔直徑d；

步驟三、基于RNG k-ε模型建立水流通過(guò)消能孔的數(shù)值模擬模型；

步驟四、對(duì)數(shù)值模擬結(jié)果進(jìn)行分析擬合，確定空腔長(zhǎng)度L與相關(guān)因素的經(jīng)驗(yàn)表達(dá)式，根據(jù)獲得的空腔長(zhǎng)度L與消能孔厚度T進(jìn)行比較，劃分孔板和洞塞。

2.根據(jù)權(quán)利要求1所述的一種基于數(shù)值模擬的孔板與洞塞劃分方法，其特征是，步驟二中確定影響空腔長(zhǎng)度L的因素具體為通過(guò)量綱分析進(jìn)一步得：L/D＝f(β，Re)，其中Re＝uDρ/μ為雷諾數(shù)，β＝d/D為孔徑比。

3.根據(jù)權(quán)利要求1所述的一種基于數(shù)值模擬的孔板與洞塞劃分方法，其特征是，步驟三中基于RNG k-ε模型建立水流通過(guò)消能孔的數(shù)值模擬模型具體為：

(1)水流通過(guò)消能孔的數(shù)值模擬模型包括以下控制方程：

連續(xù)方程：

動(dòng)量守恒方程：

k-方程：

ε-方程：

其中：x_i(＝x，y)為軸向及徑向方向坐標(biāo)，u_i(＝u_x，u_y)為軸向或徑向方向的水流速度，ρ是流體的密度，p為壓強(qiáng)，v為水流動(dòng)力粘度，v_t是渦粘度，v_t＝C_μ(k²/ε)，k表示湍流動(dòng)能，ε為湍流動(dòng)能耗散率，C_μ＝0.085，其他參數(shù)取值分別表示為η＝Sk/ε，C₁＝1.42，η_o＝4.377，λ＝0.012，C₂＝1.68，α_k＝α_ε＝1.39；

(2)水流通過(guò)消能孔的數(shù)值模擬模型的邊界條件包括：入流邊界、出流邊界、對(duì)稱軸邊界以及壁面邊界，各邊界條件的處理方法是：入流邊界條件有入流平均流速、湍流動(dòng)能分布、湍流動(dòng)能耗散率分布，它們的數(shù)學(xué)表達(dá)式分別為u_in＝u₀，k＝0.0144u₀²，ε＝k^1.5/(0.5R)，其中u₀為入口平均流速，R為泄洪洞半徑；出流邊界處理方法為假定出流得到充分發(fā)展；對(duì)稱軸邊界的處理方法是假定徑向速度為0，而且各變量沿徑向的梯度也被認(rèn)為是0；壁面邊界的處理方法為邊界層流中采用了無(wú)滑移假定，壁面邊界速度與邊界節(jié)點(diǎn)的速度分量相等。

4.根據(jù)權(quán)利要求2所述的一種基于數(shù)值模擬的孔板與洞塞劃分方法，其特征是，步驟四中空腔長(zhǎng)度L與相關(guān)因素的經(jīng)驗(yàn)表達(dá)式具體為L(zhǎng)/D＝0.4427(d/D)^-0.3538，其中d/D＝0.4-0.8，Re＞10⁵。