讓網(wǎng)格動(dòng)起來(lái)（9）—彈簧光順方法實(shí)例

上次講到,光順方法計(jì)算網(wǎng)格節(jié)點(diǎn)位移并不是依賴于網(wǎng)格質(zhì)量。因此,僅使用光順方法,只能處理一些小變形的情況 。倘若變形量過(guò)大,則會(huì)造成網(wǎng)格質(zhì)量降低,輕則影響收斂,重則出現(xiàn)負(fù)網(wǎng)格,導(dǎo)致計(jì)算出錯(cuò)。

1、問題描述

邊長(zhǎng)120mm的正方形空腔,內(nèi)有一個(gè)邊長(zhǎng)30mm的正方形固體。研究固體以角速度1rad/s旋轉(zhuǎn)時(shí)對(duì)內(nèi)部流場(chǎng)的擾動(dòng)。

1、網(wǎng)格方案

由于四邊形網(wǎng)格在剛體變形時(shí)會(huì)發(fā)生較大的扭曲,不太適合于使用彈簧光順方法,因此本例中使用三角形網(wǎng)格。使用全局尺寸2mm,生成網(wǎng)格如下圖所示。

2、指定運(yùn)動(dòng)

可以使用profile文件或DEFINE_CG_MOTION宏來(lái)定義運(yùn)動(dòng)。我們這里利用udf宏來(lái)定義運(yùn)動(dòng)。

UDF宏文件如下:

#include “udf.h”

DEFINE_CG_MOTION(RotationVel , dt , vel , omega , time , dtime)

{

omega[2] = 1;

}

omega[2]表示z方向的角速度 。本例中以恒定的角速度旋轉(zhuǎn)。

將內(nèi)部四條邊定義為rigid body。并把賦予運(yùn)動(dòng)速度。

3、動(dòng)網(wǎng)格參數(shù)

勾選Dynamic Mesh選項(xiàng),激活動(dòng)網(wǎng)格。本例只考慮網(wǎng)格光順,因此只勾選Smoothing。(說(shuō)明:通常動(dòng)網(wǎng)格應(yīng)用中不會(huì)只使用一種網(wǎng)格更新方法,此處只是考慮光順,所以才會(huì)這樣。正常情況下本例應(yīng)當(dāng)使用到網(wǎng)格重構(gòu))。

點(diǎn)擊setting…按鍵進(jìn)行smoothing設(shè)置面板 。如右下圖所示。

采用彈簧光順需要設(shè)置的參數(shù)一共有四個(gè):彈簧常數(shù)因子、邊界節(jié)點(diǎn)松弛、收斂精度及迭代次數(shù)。

通常需要設(shè)置的參數(shù)為彈簧常數(shù)因子及邊界節(jié)點(diǎn)松弛,其它兩個(gè)參數(shù)采用默認(rèn)即可。

(1)彈簧常數(shù)因子:該參數(shù)取值范圍[0,1],取值越大表示阻尼越大,也表示邊界運(yùn)動(dòng)對(duì)內(nèi)部網(wǎng)格的影響越小。

(2)邊界節(jié)點(diǎn)松弛:取值范圍[0,1],1表示無(wú)松弛,0表示禁止節(jié)點(diǎn)位移。

這里先將兩參數(shù)取中間值0.5。

3、其它設(shè)置

創(chuàng)建運(yùn)動(dòng)區(qū)域。設(shè)置中間四條邊為rigid body,設(shè)置運(yùn)動(dòng)udf(注意UDF必須為編譯型)。如下圖所示。

4、預(yù)覽網(wǎng)格運(yùn)動(dòng)

(1)觀察0.5s后的網(wǎng)格情況

可以看到,四角位置網(wǎng)格已經(jīng)扭曲得很厲害了。其它參數(shù)有興趣的同學(xué)可以自己嘗試。

5、計(jì)算結(jié)果演示