關(guān)于壁面函數(shù),無滑移邊界壁面,y+和第一層網(wǎng)格尺寸(轉(zhuǎn)部分)

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

壁面函數(shù)

1、 無論是標(biāo)準(zhǔn)k—ε模型、RNGk—ε模型,還是Realizable k—ε模型,都是針對充分發(fā)展的湍流才有效的,也就是說,這些模型均是高Re數(shù)的湍流模型。它們只能用于求解處于湍流核心區(qū)的流動。而壁面函數(shù)是對近壁區(qū)的半經(jīng)驗描述,是對某些湍流模型的補(bǔ)充(近壁區(qū)對整體流動影響較大和低雷諾數(shù)Re的情況),通過壁面函數(shù)法和低Re數(shù)k—ε模型與標(biāo)準(zhǔn)k—ε模型和RNGk—ε模型配合,成功解決整個整個管道的流動計算問題。

在壁面區(qū),流動情況變化很大。
解決這個問題目前有兩個途徑:

一、是不對粘性影響比較明顯的區(qū)域(粘性底層和過渡層)進(jìn)行求解,而是用一組半經(jīng)驗的公式(即壁面函數(shù))將壁面上的物理量與湍流核心區(qū)內(nèi)的相應(yīng)物理量聯(lián)系起來。這就是壁面函數(shù)法。在劃分網(wǎng)格的時候,不需要在壁面區(qū)加密,(只需要把第一個節(jié)點)布置在對數(shù)律成立的區(qū)域內(nèi),即配置在湍流充分發(fā)展區(qū)域。
如果要用到壁面函數(shù)的話,在define---modle--viscous面板里有near wall treatment一項。可以選擇標(biāo)準(zhǔn)壁面函數(shù)、不平衡壁面函數(shù)等。

二、是采用低Re數(shù)的k—ε模型來求解粘性底層和過渡層,此時需要在壁面區(qū)劃分比較細(xì)密的網(wǎng)格,越靠近壁面,網(wǎng)格越細(xì)。當(dāng)局部湍流的Re數(shù)小于150時,就應(yīng)該使用低Re數(shù)的k—ε模型。
總結(jié):相對于低Re數(shù)的k—ε模型,壁面函數(shù)法計算效率高,工程實用性強(qiáng)。但當(dāng)流動分離過大或近壁面流動處于高壓之下時,不是很理想。在劃分網(wǎng)格的時候,需要在壁面的位置設(shè)置邊界層網(wǎng)格,原因也是如此。

為什么要使用壁面函數(shù)呢?
首先,在CFD中應(yīng)用湍流模型并不一定需要使用壁面函數(shù),在粘性支層中可以對N-S方程直接求解。在粘性支層中,速度梯度很大,vorticity不為零,所以要直接求解,就必須在粘性支層中布置較多節(jié)點,一般要10層以上,這就是一般的低Re數(shù)湍流模型。當(dāng)然這樣將占用較多的計算資源。
而在邊界層中,是存在解析解的,如果在粘性支層內(nèi)不求解三維N-S方程,而用一維數(shù)學(xué)模型代替,將大大降低計算資源的使用,這就是壁面函數(shù)。一般高Re數(shù)湍流模型都使用壁面函數(shù)。第一層網(wǎng)格節(jié)點布置在粘性支層之外。那么你如何判斷你的邊界層網(wǎng)格節(jié)點布置是否合適呢?這就要檢查你的y+,y+就是第一層網(wǎng)格質(zhì)心到壁面的無量綱距離,與速度、粘度、剪應(yīng)力等等都有關(guān)系。對于y+的值,各個學(xué)者推薦的范圍是不一樣的,但一般在30-60之內(nèi)肯定是沒有問題的。也有推薦10-110甚至200的。y+的值合理,意味著你的第一層邊界網(wǎng)格布置比較合理,如果y+不合理,就要調(diào)整你的邊界層網(wǎng)格


無滑移邊界條件:壁面處法向和切向速度為零。

滑移邊界條件:壁面處法向速度為零和切向速度不為零



一個成功的湍流計算離不開好的網(wǎng)格。在許多的湍流中,空間的有效粘性系數(shù)不同,是平均動量和其它標(biāo)量輸運(yùn)的主要決定因素。因此,如果需要有足夠的精度,這就需要保證湍流量要比較精確求解。由于湍流與平均流動有較強(qiáng)的相互作用,因此求解湍流問題比求解層流時候更依賴網(wǎng)格。對于近壁網(wǎng)格而言,不同的近壁處理對網(wǎng)格要求也不同。下面對常見的幾種近壁處理的網(wǎng)格要求做個說明。采用壁面函數(shù)時候的近壁網(wǎng)格:第一網(wǎng)格到壁面距離要在對數(shù)區(qū)內(nèi)。對數(shù)區(qū)的y+ >30~60。FLUENT在y+ <12.225時候采用層流(線性)準(zhǔn)則,因此網(wǎng)格不必要太密,因為壁面函數(shù)在粘性底層更本不起作用。對數(shù)區(qū)與完全湍流的交界點隨壓力梯度和雷諾數(shù)變化。如果雷諾數(shù)增加,該點遠(yuǎn)離壁面。但在邊界層里,必須有幾個網(wǎng)格點。 壁面函數(shù)處理時網(wǎng)格劃分采用雙層模型時近壁網(wǎng)格要求當(dāng)采用雙層模型時,網(wǎng)格衡量參數(shù)是y+ ,并非y* 。最理想的網(wǎng)格劃分是需要第一網(wǎng)格在y+ =1位置。如果稍微大點,比如 =4~5,只要位于粘性底層內(nèi),都是可以接收的。理想的網(wǎng)格劃分需要在粘性影響的區(qū)域內(nèi)(Rey<200 )至少有十個網(wǎng)格,以便可以計算粘性區(qū)域內(nèi)的平均速度和湍流量。 采用雙層區(qū)模型時網(wǎng)格劃分采用Spalart-Allmaras 模型時的近壁網(wǎng)格要求該模型屬于低雷諾數(shù)模型。這就要求網(wǎng)格能滿足求解粘性影響區(qū)域內(nèi)的流動,引入了阻尼函數(shù),用以削弱粘性底層的湍流粘性影響。因此,理想的近壁網(wǎng)格要求和采用雙層模型時候的網(wǎng)格要求一致。采用大渦模擬的近壁網(wǎng)格要求對于大渦模擬,壁面條件采用了壁面法則,因此對近壁網(wǎng)格劃分沒有太多限制。但是,如果要得到比較好的結(jié)果,最好網(wǎng)格要細(xì),最近網(wǎng)格距離壁面在 y+=1的量級上。 for Hexa mesh, ==>Y+是第一層高度一半和 viscous length scale 的比值 for Tetra mesh==>Y+是第一層高度1/3和 viscous length scale 的比



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

相關(guān)標(biāo)簽搜索:關(guān)于壁面函數(shù),無滑移邊界壁面,y+和第一層網(wǎng)格尺寸(轉(zhuǎn)部分) Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

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

全國服務(wù)熱線

1358-032-9919

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