基于ICEM的CFX網格重劃分

2016-10-08 by:CAE仿真在線來源:互聯網

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

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

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

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

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

創建移動壁之外的其它邊界,并將網格運動設為unspecified。這將允許這些邊界面上的節點在需要的時候可以移動。

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

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

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

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

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

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

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

開放分享：優質有限元技術文章,助你自學成才