??Fluent版本：19.0（其他版本應(yīng)該也適用）

這里我們用一個(gè)簡單的算例（同心環(huán)中的自然對(duì)流）來說明

算例來自《ANSYS Fluid Dynamics Verification Manual》中的VMFL009: Natural Convection in a Concentric Annulus.?

外環(huán)溫度為327K，內(nèi)環(huán)溫度為373K

圓環(huán)內(nèi)流體的物性參數(shù)為：

?

通過計(jì)算瑞利數(shù)，本算例應(yīng)該采用層流模型進(jìn)行計(jì)算

網(wǎng)格模型如下：

?

首先我們不考慮重力

打開Fluent，讀入網(wǎng)格（建議在ICEM導(dǎo)出網(wǎng)格的時(shí)候就進(jìn)行縮放，這樣可以避免在Fluent當(dāng)中scale了）

?

計(jì)算很快完成

溫度云圖如下：

?

?　　我們可以看出自然對(duì)流在沒有施加重力的作用下，和導(dǎo)熱幾乎是一樣的

?接下來我們考慮重力對(duì)自然對(duì)流的影響，我們?cè)?/span>Y的負(fù)向施加9.81m/s²的重力加速度，只需在上面算例的基礎(chǔ)上勾選重力選項(xiàng)，然后輸入對(duì)應(yīng)的重力加速度

?那么如果我們重力是隨時(shí)間變化的怎么辦呢？我們可以采用下面三種方式來讓重力隨時(shí)間變化

方法一：采用UDF當(dāng)中的DEFINE_ADJUST宏來實(shí)現(xiàn)重力隨時(shí)間變化，我們?nèi)匀徊捎帽舅憷?#xff0c;我們?nèi)匀还催x重力選項(xiàng)（如果不勾選，則不能改變重力），但是我們不輸入重力加速度的值，而采用默認(rèn)值0，在UDF中指定重力，我們編寫如下的代碼：

編譯該UDF，掛載該UDF

然后

如果沒有上面的步驟，否則UDF指定的重力將不生效

通過對(duì)面我們可以看出GUI施加重力和UDF施加重力的效果是一樣的，計(jì)算結(jié)果是完全一致的，如果我們想要重力隨時(shí)間變化，我們可以在UDF中采用RP_Get_Real(“flow-time”)或者CURRENT_TIME來根據(jù)時(shí)間變化調(diào)整重力方向

方法二：通過UDF在Fluent當(dāng)中添加到動(dòng)量源項(xiàng)的方式

我們編寫如下的代碼：

編譯該UDF

掛載該UDF

然后計(jì)算

我們可以看到通過添加動(dòng)量源項(xiàng)的方式與GUI施加重力的效果也是一致的，施加動(dòng)量源項(xiàng)需要注意一下量綱，返回源項(xiàng)的值時(shí)，需要重力與源項(xiàng)之間相差一個(gè)密度，需要用重力乘密度以后再返回。如果想要重力隨時(shí)間變化，同上，可以在UDF中采用RP_Get_Real(“flow-time”)或者CURRENT_TIME來根據(jù)時(shí)間變化調(diào)整重力方向。

方法三：通過編寫scheme來實(shí)現(xiàn)改變重力

我們編寫如下的scheme代碼：

(define (ChangeGravity)

(begin

??? (rpsetvar 'gravity/y -9.81)

??? (%models-changed)

)

)

然后將文本文件重命名為以.scm為后綴名（我命名為demo.scm）

我們?nèi)匀还催x重力選項(xiàng)（如果不勾選，則不能改變重力），但是我們不輸入重力加速度的值，而采用默認(rèn)值0。

然后計(jì)算

從計(jì)算結(jié)果我們看到和scheme施加重力的方式和GUI施加重力的方式也是一致的，隨時(shí)間變化在scheme當(dāng)中為(rpgetvar 'flow-time)，然后根據(jù)獲取的時(shí)間來改變重力

轉(zhuǎn)載于:https://www.cnblogs.com/liusuanyatong/p/11259629.html

總結(jié)

以上是生活随笔為你收集整理的fluent在运行时改变重力方向方法总结的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。