關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)

2017-02-07  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

1、為什么會存在第一層網(wǎng)格的問題

關(guān)于這個問題,實(shí)際上要從邊界層說起。實(shí)驗(yàn)表明,邊界層內(nèi)根據(jù)流動狀態(tài)的不同可以分為三層,自壁面向流動核心區(qū)分別為:粘性子層、過渡層和湍流核心層。邊界層很薄,一般都是毫米~微米級,因此,若采用劃分網(wǎng)格進(jìn)而利用數(shù)值方法求解的話,勢必會大大增加計(jì)算網(wǎng)格的數(shù)量,從而急劇增加計(jì)算工作量。又有實(shí)驗(yàn)發(fā)現(xiàn),在粘性子層和過渡層內(nèi),主要是粘性力在起主導(dǎo)作用,慣性力的作用幾乎可以忽略。在該區(qū)域內(nèi),粘性力與速度梯度成線性關(guān)系,因此在于核心層為高雷諾數(shù)湍流流動的情況下,過渡層與粘性子層內(nèi)的速度分布可以通過經(jīng)驗(yàn)公式直接計(jì)算得到,而無需劃分網(wǎng)格,換句話說,在這種情況下,可以將計(jì)算節(jié)點(diǎn)的第一層網(wǎng)格節(jié)點(diǎn)放置在湍流核心區(qū)內(nèi),而過渡層與粘性子層中則無需要任何網(wǎng)格。這部分區(qū)域中的物理量分布采用壁面函數(shù)(wall function)來計(jì)算完成。需要用到壁面函數(shù)的湍流模型包括:k-epsilon模型,雷諾應(yīng)力模型。

另外一種低雷諾數(shù)湍流模型的情況則與之不同,其不采用壁面函數(shù)來求解粘性子層與過渡層中的流動物理量分布,而是采用NS方程離散求解,與核心區(qū)域求解方式一樣,如K-W模型,SA模型等。

2、Y+的問題

Y+是什么玩意兒?Y+其實(shí)是一個無量綱量,其定義為:

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent流體分析圖片1

式中,u*為近壁面摩擦速度(friction velocity),為第一層網(wǎng)格節(jié)點(diǎn)與壁面的間距,υ為流體的運(yùn)動粘度。

其中壁面摩擦速度關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent流體分析圖片2
式中τw為壁面剪切應(yīng)力,其值為關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent分析圖片3,其中μ為動力粘度,為第一層網(wǎng)格間距。

因此可以估算第一層網(wǎng)格間距:

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent分析圖片4

3、更簡單的計(jì)算方式

采用上式進(jìn)行第一層網(wǎng)格間距計(jì)算比較麻煩,因?yàn)樾枰?jì)算u*,而u*的計(jì)算又涉及到壁面剪切應(yīng)力的計(jì)算,壁面剪切應(yīng)力的計(jì)算又涉及到速度梯度的計(jì)算。麻煩的事情在于壁面法向速度梯度在劃分網(wǎng)格的時(shí)候是未知的,只有在計(jì)算完畢后才能得到,這實(shí)在是打臉的行為。那用什么辦法去補(bǔ)救呢?工程應(yīng)用中,壁面剪切應(yīng)力通常采用估計(jì)值。

引入一個新的物理量:壁面摩擦系數(shù)(skin friction coefficient):

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent分析圖片5

通過計(jì)算Cf的值,從而可以計(jì)算出壁面剪切應(yīng)力τw

Cf的計(jì)算方式有很多種,如下表所示。

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent分析圖片6

上面計(jì)算公式中,有顯式也有隱式,目的都是計(jì)算Cf。顯式可以直接計(jì)算,隱式可以采用迭代法求解計(jì)算。

式中Re為雷諾數(shù),其表達(dá)式為:

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent培訓(xùn)的效果圖片7

式中ρ為流體密度,U為速度,L為邊界層參考尺寸,μ為流體的動力粘度。

計(jì)算得到Cf后,即可計(jì)算壁面剪切應(yīng)力τw:

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent培訓(xùn)的效果圖片8

從而可以計(jì)算:

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent培訓(xùn)的效果圖片9

進(jìn)而可以計(jì)算出第一層網(wǎng)格高度:


關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent培訓(xùn)的效果圖片10

4、Y+的取值

從上面計(jì)算第一層網(wǎng)格高度的公式可以看出,我們需要自己提供Y+值,那么這個值應(yīng)該給多少呢?

一般來說,對于高雷諾數(shù)模型(如k-Epsilon模型、雷諾應(yīng)力模型等),需要滿足

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent培訓(xùn)課程圖片11

一般以接近30為佳。

對于低雷諾數(shù)模型(如k-w模型,SA模型,LES等),需要滿足,以接近于1為佳。

所以在估算第一層網(wǎng)格時(shí),按選擇使用的湍流模型的不同,通常取301進(jìn)行估算。

5、更簡單的方法

萬能的互聯(lián)網(wǎng)給我們提供了很多的便利,實(shí)際上網(wǎng)絡(luò)上有很多現(xiàn)成的計(jì)算Y+的工具,這里推薦PointwiseCFD-online的計(jì)算器。

(1)pointwisey+計(jì)算工具網(wǎng)址為http://www.pointwise.com/yplus/,打開后如下圖所示,輸入速度、密度、粘度、特征尺寸以及Y+,網(wǎng)頁會計(jì)算出第一層網(wǎng)格高度與雷諾數(shù)。

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent培訓(xùn)課程圖片12

(2)CFD-Online的計(jì)算工具網(wǎng)址為http://www.cfd-online.com/Tools/yplus.php,打開后如下圖所示。

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent培訓(xùn)課程圖片13

網(wǎng)絡(luò)上常見的NASAY+計(jì)算器http://geolab.larc.nasa.gov/APPS/YPlus/,如果不是做空氣外流場計(jì)算的話,不推薦使用。從下圖可以看出,該計(jì)算器并沒有密度和粘度項(xiàng),其實(shí)其默認(rèn)了介質(zhì)為空氣,粘度為1.7894e-5,溫度273.15,ν =1.4,密度通過狀態(tài)方程計(jì)算得到。不信的話自己看源代碼。

關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì)fluent分析案例圖片14

6、遺留的問題

這里談的第一層網(wǎng)格間距估計(jì)是在輸入已知Y+的情況下獲得的,實(shí)際上在計(jì)算完后還需要檢查壁面的Y+分布,看是否滿足湍流模型的要求,如果不滿足的話,還需要重新劃分網(wǎng)格,重新計(jì)算,重新檢查。不斷進(jìn)行下去,直到滿足Y+要求(高雷諾數(shù)30~300,低雷諾數(shù)<1)。

7、題外話

Y+這個東西其實(shí)挺麻煩的,實(shí)際工程應(yīng)用中沒那么嚴(yán)格,除非研究型的計(jì)算。(待商榷)




開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:關(guān)于CFD中的Y+與第一層網(wǎng)格厚度估計(jì) Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報(bào)名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進(jìn)行在線報(bào)名,我們會及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com