ICEM中不同類型網(wǎng)格交界面節(jié)點(diǎn)匹配的處理技巧

2016-09-01 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

眾所周知(很鄙視這個(gè)詞開頭的文章,一看就覺(jué)得作者沒(méi)水平),CFX和Fluent是目前比較流行的兩款CFD流體分析工具,而作為前處理的網(wǎng)格工具ICEM可以為它們提供優(yōu)秀的結(jié)構(gòu)化和非結(jié)構(gòu)化網(wǎng)格。通常在旋轉(zhuǎn)流體機(jī)械(如泵、攪拌器、風(fēng)力機(jī)等)的CFD分析中存在動(dòng)靜區(qū)域網(wǎng)格的交接面,有兩種方式定義交界面邊界條件,即interior和interface方式。下面首先對(duì)其進(jìn)行簡(jiǎn)單說(shuō)明(別嫌筆者像唐僧一樣啰嗦,很多網(wǎng)友都沒(méi)搞清楚哦)。

通俗地說(shuō)interior方式即用一個(gè)“假想的內(nèi)部墻”將計(jì)算域分隔,這個(gè)墻的臉皮薄到厚度為0,能使任何介質(zhì)穿過(guò),且墻面上只有一種網(wǎng)格形式,即所連接的計(jì)算域的網(wǎng)格節(jié)點(diǎn)是共用的。而interface方式相當(dāng)于“切蘋果”,因此分割后的兩個(gè)計(jì)算域均多了一個(gè)面邊界,即成對(duì)出現(xiàn)的interface。Interface方式的好處在于可以對(duì)連接的兩個(gè)計(jì)算域劃分不同形式的網(wǎng)格,比如一個(gè)劃分四面體網(wǎng)格,另一個(gè)劃分六面體網(wǎng)格,那么這兩個(gè)計(jì)算域的interface上分別為三角形網(wǎng)格和四邊形網(wǎng)格。這樣,不需要兩者在交界面上的網(wǎng)格節(jié)點(diǎn)重合,求解器依然可以進(jìn)行插值計(jì)算。另一方面,一對(duì)interface不要求完全相同,即一個(gè)大面和一個(gè)小面連接同樣可以計(jì)算。但是需要注意的是,為了提高計(jì)算精度,一對(duì)interface的面網(wǎng)格尺寸要盡量接近。interior方式的好處在于交界面(只有一個(gè)面)網(wǎng)格節(jié)點(diǎn)是共用的,因此避免了interface方式帶來(lái)的插值誤差,但在許多方面也存在很多缺點(diǎn)。比如用Fluent對(duì)旋轉(zhuǎn)機(jī)械的內(nèi)部流場(chǎng)進(jìn)行非穩(wěn)態(tài)計(jì)算時(shí),若采用interior方式連接交界面,則不能使用滑移網(wǎng)格模型導(dǎo)致求解失敗。據(jù)筆者了解,至少ANSYS Fluent14.5及以下版本均存在這個(gè)問(wèn)題,筆者嘗試過(guò)在fluent中將一個(gè)interior面分割成一對(duì)interface但失敗了(如果可行,求大神賜教!)。然而,若用CFX求解就不會(huì)存在這個(gè)問(wèn)題,因?yàn)镃FX認(rèn)為一個(gè)內(nèi)部面有兩個(gè)“side”,也就是說(shuō)CFX中可以將interior“分割成一對(duì)interface”。

由上述可知,不論網(wǎng)格節(jié)點(diǎn)對(duì)應(yīng)與否,似乎CFX求解更有優(yōu)勢(shì)一些。但是筆者的強(qiáng)迫癥已深入骨髓不可救藥,非要拿fluent軟件采用滑移網(wǎng)格模型來(lái)計(jì)算用不同網(wǎng)格形式連接但交界面上網(wǎng)格節(jié)點(diǎn)匹配的流場(chǎng)(對(duì)不起,哥哥不是神經(jīng)病)。那么問(wèn)題來(lái)了,挖掘機(jī)技術(shù)哪家強(qiáng)?

預(yù)知挖掘機(jī)到底哪家強(qiáng),不用下回分講,這就劇透給你,其實(shí)還是ICEM的網(wǎng)格技術(shù)強(qiáng)。沒(méi)事多睡覺(jué),無(wú)聊少扯淡。拋青磚,引紅磚,下面開始說(shuō)正點(diǎn)。天下網(wǎng)友靠不住,有病還得自己治。

簡(jiǎn)化一下問(wèn)題吧,如圖1所示的模型,計(jì)算域分為左右兩部分,流體左進(jìn)右出?，F(xiàn)在要對(duì)左邊小域劃分四面體非結(jié)構(gòu)網(wǎng)格,右邊大域劃分六面體結(jié)構(gòu)化網(wǎng)格,中間的交界面coupling上要求節(jié)點(diǎn)重合,但不是interior方式連接,而是一對(duì)interface方式連接。

圖1 計(jì)算域示意圖

直奔主題太過(guò)倉(cāng)促,畢竟筆者研究這個(gè)也花了一些心思,還是娓娓道來(lái)比較好,多讓我感受一下小小成就感吧。

一、節(jié)點(diǎn)匹配的同種類型網(wǎng)格

創(chuàng)建兩個(gè)body,對(duì)兩段計(jì)算域均同時(shí)劃分四面體網(wǎng)格(或者六面體網(wǎng)格),則coupling面上都只存在一種網(wǎng)格,即節(jié)點(diǎn)匹配。輸出msh并導(dǎo)入fluent后,邊界條件里面的交界面除了coupling面以外,還會(huì)自動(dòng)生成一個(gè)coupling-shadow面,兩者均默認(rèn)為wall類型,但均不能被設(shè)置為interface類型,將coupling邊界設(shè)置為interior后,coupling-shadow自動(dòng)消失。這種方法比較常用,因此不予贅述。

二、節(jié)點(diǎn)不匹配的不同網(wǎng)格類型

對(duì)左邊計(jì)算域劃分四面體網(wǎng)格,右邊域劃分六面體結(jié)構(gòu)網(wǎng)格(如圖2),可以發(fā)現(xiàn)coupling面上的左右節(jié)點(diǎn)是不重合的,即coupling面上存在兩種網(wǎng)格,一種三角形網(wǎng)格,一種四邊形網(wǎng)格。輸出msh并導(dǎo)入fluent后,邊界條件里面的交界面除了coupling面以外,還會(huì)自動(dòng)生成一個(gè)coupling:002面,均默認(rèn)為wall,一個(gè)是三角形網(wǎng)格,一個(gè)是四邊形網(wǎng)格?？梢詫烧呔O(shè)置為interface類型并在“Mesh Interfaces”選項(xiàng)中建立交界面數(shù)據(jù)傳遞。

主要步驟:

(1)對(duì)計(jì)算域各個(gè)邊界建立part后,僅對(duì)左邊域建立body1,這樣生成四面體網(wǎng)格的時(shí)候軟件只對(duì)建立了body的部分劃分網(wǎng)格,檢查質(zhì)量并優(yōu)化;

(2)創(chuàng)建初始block,命名為body2,選擇右邊域的柱面,僅對(duì)右邊域創(chuàng)建初始block;

(3)創(chuàng)建O-block,建立映射,生成右邊域的六面體結(jié)構(gòu)網(wǎng)格,檢查質(zhì)量;

(4)轉(zhuǎn)換為非結(jié)構(gòu)網(wǎng)格,并與左邊四面體網(wǎng)格合并。

a.全局網(wǎng)格 b. coupling處局部顯示 c.交界面的兩種網(wǎng)格

圖2計(jì)算域網(wǎng)格

三、節(jié)點(diǎn)匹配的不同網(wǎng)格類型

對(duì)左邊計(jì)算域劃分四面體網(wǎng)格,右邊域劃分六面體結(jié)構(gòu)網(wǎng)格,但通過(guò)節(jié)點(diǎn)匹配使交界面上只有一種網(wǎng)格,這樣在左邊域中會(huì)形成五面體網(wǎng)格。網(wǎng)格生成步驟為上述4步之后再進(jìn)行如下操作。

a)合并節(jié)點(diǎn),執(zhí)行Edit Mesh—Merge Nodes—Merge Meshes命令,勾選交界面(coupling)的part,點(diǎn)擊Apply。

b)由于合并節(jié)點(diǎn)后,coupling面上僅剩下四邊形網(wǎng)格,左邊域中生成五面體網(wǎng)格以連接兩段計(jì)算域。這將導(dǎo)致網(wǎng)格質(zhì)量下降,因此對(duì)其進(jìn)行光順。最后,節(jié)點(diǎn)合并后的網(wǎng)格如圖3所示。

圖3計(jì)算域內(nèi)部網(wǎng)格

這時(shí),如果導(dǎo)出msh并導(dǎo)入fluent,交界面將同第一種類型,不能被改成interface類型,與其不同的僅僅是coupling面上為四邊形網(wǎng)格而已。這種方法實(shí)用性很強(qiáng),對(duì)多區(qū)計(jì)算域中拓?fù)浣Y(jié)構(gòu)比較簡(jiǎn)單的部分劃分六面體結(jié)構(gòu)化網(wǎng)格(如旋轉(zhuǎn)機(jī)械流體計(jì)算域進(jìn)出口延伸管段和翼型繞流的外場(chǎng)等區(qū)域)可以在保證節(jié)點(diǎn)匹配的情況下盡可能減少網(wǎng)格數(shù)目。但是,只要交界面類型設(shè)置為interior(或者說(shuō)不能被改為一對(duì)interface)在fluent中就不能采用滑移網(wǎng)格模型而不能順利進(jìn)行瞬態(tài)流場(chǎng)計(jì)算。鑒于此,還需將交界面網(wǎng)格進(jìn)行分割,處理為interface方式。方法如下。

在執(zhí)行上述命令后,執(zhí)行Edit Mesh—Split Mesh—Split Internal Wall命令,勾選coupling的part,點(diǎn)擊Apply。這時(shí),ICEM的模型樹中,COUPLING下面多了一個(gè)part(叫COUPLING_BACK),通過(guò)觀察可以發(fā)現(xiàn)交界面被分割成了兩部分,導(dǎo)入fluent后,可以將二者設(shè)置為interface。至此,就能保證網(wǎng)格模型能夠適應(yīng)fluent進(jìn)行瞬態(tài)流場(chǎng)計(jì)算的需要了。

實(shí)際上,執(zhí)行這一步“Split”命令的效果等同于在Gambit網(wǎng)格劃分中不勾選“connect”而生成的網(wǎng)格模型。順便說(shuō)明一下,在前述第一種類型下執(zhí)行此命令就可以獲得一對(duì)interface的效果,博主故意放在最后只是想考驗(yàn)下大家的耐心而已,并虛心并熱情地?fù)肀СＳ么朔椒ǖ耐械膹?qiáng)烈鄙視。

很榮幸您已經(jīng)讀到這里了,為了對(duì)您表示感謝,順便再向您展示一點(diǎn)我最近剛學(xué)的小技巧,可能是多數(shù)人已經(jīng)用慣了的老知識(shí)吧。歡迎吐槽!

Fluent轉(zhuǎn)換多面體網(wǎng)格的方法(但這個(gè)方法僅對(duì)四面體網(wǎng)格有效哦!)

采用多面體網(wǎng)格最大的好處是能夠減少網(wǎng)格數(shù)目,但是至今不是很流行想必還有很多原因,希望慎用!多面體網(wǎng)格生成原理如下:我們知道對(duì)一個(gè)三角形進(jìn)行Y型劃分(類似通過(guò)內(nèi)心向各邊做垂線)可以分割為三個(gè)四邊形,那么對(duì)一系列相鄰的三角形均做Y型劃分,這樣原來(lái)三角形的端點(diǎn)(即四面體網(wǎng)格的節(jié)點(diǎn))就被進(jìn)行Y型劃分后的邊所組成的多邊形包圍,然后刪除原來(lái)的節(jié)點(diǎn)并由多邊形的端點(diǎn)構(gòu)成新的網(wǎng)格節(jié)點(diǎn),這樣原來(lái)的三角形網(wǎng)格就被多邊形網(wǎng)格所替代,同理,四面體網(wǎng)格轉(zhuǎn)換為多面體網(wǎng)格。沒(méi)事可以畫著玩玩,如果你媽不打你!

方法:對(duì)導(dǎo)入的msh文件,執(zhí)行Mesh—Polyhedra—Convert Domain命令,將四面體網(wǎng)格轉(zhuǎn)換為多面體網(wǎng)格,可以減少網(wǎng)格數(shù)至原來(lái)的1/5~1/3,對(duì)于計(jì)算機(jī)配置較低的童鞋很實(shí)用哦!轉(zhuǎn)換完之后,再執(zhí)行Mesh—Reorder—Domain命令對(duì)網(wǎng)格重排序,直到Bandwidth reduction≤1(帶寬壓縮,偶也不知道這是啥玩意),能夠使網(wǎng)格排序更能適應(yīng)流動(dòng)結(jié)構(gòu),提高收斂性和計(jì)算精度。最后,附上網(wǎng)格圖一張。

附圖多面圖網(wǎng)格

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