基于ICEM的CFX網(wǎng)格重劃分

2016-10-08 by:CAE仿真在線來源:互聯(lián)網(wǎng)

在很多流-構(gòu)交互(FSI)的CFD問題中,大變形非常普遍。通常CFD求解器允許網(wǎng)格變形以適應(yīng)這種變形。然而,如果這種變形太大,不破壞這些網(wǎng)格或不使太多單元發(fā)生扭曲幾乎是不可能的?？朔@個問題的方法就是當(dāng)網(wǎng)格質(zhì)量變得很差時,進行網(wǎng)格重劃分。
創(chuàng)建輸入網(wǎng)格
這部分工作需要在ICEMCFD中完成,因為我們將要用它進行網(wǎng)格重劃分。在這個例子中,我要讓一個矩形盒子在方形通道中向下移動
圖1顯示的幾何最初在DesignModeler中創(chuàng)建。首先我們需要創(chuàng)建幾個CFD仿真中要用到幾個邊界。位于中心的盒子是向下移動的目標(biāo)。給它一個名字:BOX。就像所有ICEM CFD中的網(wǎng)格,你必須制定必要的仿真標(biāo)準(zhǔn)。

其次,在ICEM CFD中的Replay Manager執(zhí)行必要的網(wǎng)格劃分,將這些命令保存為一個重放(回放)文件(.rpl),這個重放文件以后要在CFX中用到。然后再將創(chuàng)建好的網(wǎng)格以CFX格式導(dǎo)出,并保存這個項目,記住ICEM CFD文件(.tin和.rpl)所在的位置。

設(shè)置CFX CASE
一般設(shè)置
--在這個例子中BOX隨時間向下移動。
--用希望的物理屬性定義域中的流體流動,并允許網(wǎng)格變形(mesh deformation),將網(wǎng)格剛度設(shè)置為increase in small volumes。這樣將保持小的棱鏡單元盡可能地靠近實體壁面,主要變形會被遠(yuǎn)離Body的部分吸收。

邊界細(xì)節(jié)
--頂面、底面使用開放條件(Opening),使流體可以按照需要進入或流出計算域。
--這種情況下,基本部分允許網(wǎng)格變形并偵測網(wǎng)格質(zhì)量。當(dāng)網(wǎng)格質(zhì)量變得太低時,實現(xiàn)網(wǎng)格重劃分的調(diào)用。這里,壁面運動用一個表達式來描述。

--創(chuàng)建一個表達式描述實體的運動。這里使用一個隨時間增加的新型運動。表達式的名稱命名為"motion"。之所以要使用Time This Run是因為網(wǎng)格重劃分后求解器重新啟動會獲得大跳躍。另外,可對這個變化的網(wǎng)格重新初始化時間。

創(chuàng)建移動壁之外的其它邊界,并將網(wǎng)格運動設(shè)為unspecified。這將允許這些邊界面上的節(jié)點在需要的時候可以移動。

BOX的邊界與其它邊界有所不同。在這里運動表達式應(yīng)該設(shè)定為z-axis位移。因為BOX將向下移動,因此應(yīng)該在表達式前加上負(fù)號"-"。

中斷控制設(shè)置
接下來設(shè)置一個停止控制,這將偵測網(wǎng)格質(zhì)量并停止求解器的計算。
--首先創(chuàng)建一個邏輯表達式來偵測網(wǎng)格的正交性。表達式命名為:remeshingcond。如果正交性變得小于5度,則表達式返回值為1,否則返回值為0。

--在求解器控制標(biāo)簽下的中斷控制中,添加這個表達式。

--指定一個配置允許仿真采用這個指定的事件。例如,一個配置可以允許一系列的仿真或網(wǎng)格重劃。

Re-mesh控制設(shè)置
在定義re-meshing之前需要定一兩個輸出偵測。因為Re-meshing特征要使用輸出偵測作為ICEM CFD的輸入。
1)第一個輸出控制為GeometryScale factor。這將衡量CFX輸出,使其與ICEM CFD中的模型相匹配。這里geometry scale設(shè)為1mm,因為ICEM模型是以mm為單位建立的。如果這個模型是按照相同的尺寸建立的,這個控制就不需要了。
2)第二個輸出偵測為WALLDisplacement z。這個輸出偵測跟蹤移動邊界總平均質(zhì)心位移。在12版中增加了這個變量,從仿真開始便開始跟蹤曲面的質(zhì)心。因此,這個偵測從仿真的第一個時間步開始就會給我們提供BOX的總運動。

利用定義的這兩個輸出控制現(xiàn)在就能指定網(wǎng)格重劃分(re-meshing)了:
3)利用ICEM CFDReplay選項創(chuàng)建一個新的re-mesh。
4)這個re-mesh應(yīng)該由我們前面創(chuàng)建的中斷條件激活。
5)在這里應(yīng)該被重劃網(wǎng)格的域是:BODY。這個Body對應(yīng)ICEM CFD中的體(volume)。
6)幾何文件應(yīng)該是在網(wǎng)格化分時創(chuàng)建的.tin文件。
7)回放文件是劃分網(wǎng)格時存儲分網(wǎng)步驟的那個.rpl文件。
接下來要創(chuàng)建一個ICEMCFD Part Map。這將使CFX中的邊界與ICEM CFD中的邊界配對。這種情況下我們只需要創(chuàng)建一個part映射,因為我們只有一個移動邊界。

最后需要創(chuàng)建兩個標(biāo)量參數(shù)。這些標(biāo)量參數(shù)的名稱是敏感的,因為它們要由CFX安裝的腳本識別。
第一個命名為DISPLACEMENTZ BOX,并應(yīng)使用上面定義的WALLDisplacement z輸出控制。
第二個標(biāo)量參數(shù)應(yīng)該命名為ICEMCFD Geometry Scale,并設(shè)置相應(yīng)的偵測點。

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