[案例分析]基于商用軟件FLUENT的LES(大渦模擬)計算教學

寫在前面

LES的計算中,實際上對網(wǎng)格是有要求滴,這方面內(nèi)容可以從相關文獻中找到,本文只是針對LES的計算設置進行一個簡單的2D圓柱擾流講解,不涉及網(wǎng)格要求方面,童鞋們要注意這一點哇!

文主花了兩天時間學習FLUENT中的LES計算,所以,這實際上還是蠻簡單的,只要大體思路成型了,剩下的只是細節(jié)工作。

這個就是文主計算的結(jié)果(延伸段實際上應該更長一些)

寫稿初衷

本文的寫稿初衷是因為當初在各大網(wǎng)站上苦苦找尋LES設置算例,然而。。。木有找到,因此就想做一個基于商用軟件FLUENT的LES教程。

選擇FLUENT的原因是因為目前大多數(shù)童鞋都比較喜歡使用FLUENT來進行流動數(shù)值模擬。鑒于FLUENT的受歡迎程度以及初學者們的需求,本文就基于FLUENT做一次LES計算的教程。

適宜人群:想學習LES計算、流動非定常計算、FLUENT的筒子們

文主使用的軟件:ICEM CFD15.0、FLUENT 15.0

算例:二維(2D)圓柱擾流計算

  • First Step:前處理

前處理用一句話來概括就是:準備計算網(wǎng)格!

網(wǎng)格這一塊不是LES計算的重點,因為任何計算都要畫網(wǎng)格,因此建議童鞋們可以參考其他教程單獨學習畫網(wǎng)格,本教程只作簡單的介紹。

由于圓柱擾流問題比較簡單,因此可以直接在ICEM上畫,思想是:由點構(gòu)成線,再由線構(gòu)成面。

步驟是:

  1. 先給出幾何點(比如圓柱的圓心、流域的邊界點);

  2. 通過點連線,最終得到面(如下圖所示)

到這里就可以畫網(wǎng)格了,可以使用非結(jié)構(gòu)網(wǎng)格劃分或者結(jié)構(gòu)化網(wǎng)格劃分,兩種網(wǎng)格FLUENT均能計算。

文主一般比較喜歡用結(jié)構(gòu)化網(wǎng)格,所以在這里展示一下如何畫結(jié)構(gòu)化網(wǎng)格。

總體思路是:先給出塊(Block),再分割塊(Block),然后給網(wǎng)格節(jié)點,最后生成網(wǎng)格并轉(zhuǎn)成非結(jié)構(gòu)化網(wǎng)格(ICEM只能導出非結(jié)構(gòu)化網(wǎng)格)。

步驟是:

  1. 給出塊(Block),用一個Block把所有的面(Surface)包住;

  2. 分割塊(Block),這里解釋一下為什么要分割:圓柱擾流問題中,圓柱(固體)部分是不計算的,因此流體域中肯定會缺一個圓,這個缺口附近的網(wǎng)格很難把握,為了保證網(wǎng)格質(zhì)量,因此建議用一個正方形的Block包住這個圓,所以把長方形的Block分成一個正方形和一個小長方形;

  3. OGrid Block劃分圓柱附近的網(wǎng)格,并刪除掉中心區(qū)域(中心區(qū)域為圓柱);

  4. Block的Edge(也就是塊的邊)與幾何的邊(Curve,中文意思曲線)關聯(lián)起來

  5. 設置Edge的網(wǎng)格節(jié)點數(shù)以及分布;

  6. 最后一步很關鍵,將結(jié)構(gòu)化網(wǎng)格轉(zhuǎn)化為非結(jié)構(gòu)化網(wǎng)格(文主曾經(jīng)被困在這里很久,后來才知道ICEM只能導出非結(jié)構(gòu)化網(wǎng)格到FLUENT里面)

  7. 最終生成網(wǎng)格如下圖:

[案例分析]基于商用軟件FLUENT的LES(大渦模擬)計算教學的圖1


接下來就是導入FLUENT進行計算。

  • Second Step:FLUENT計算

總的思路是:

  1. 先導入網(wǎng)格,然后設置邊界條件;

  2. 進行定常計算(Steady),提供一個比較好的初始流場;

  3. 在選擇非定常計算(Transient),選擇LES-2D模式(這里后面會講);

  4. 然后將非定常的結(jié)果輸出就好的(這里講的時候要細看哇)

下面是具體操作。

打開FLUENT,因為是計算2D圓柱擾流,因此這里選擇2D格式。

[案例分析]基于商用軟件FLUENT的LES(大渦模擬)計算教學的圖2

導入網(wǎng)格并顯示

[案例分析]基于商用軟件FLUENT的LES(大渦模擬)計算教學的圖3

設置邊界條件(比如進口、出口、壁面等)

[案例分析]基于商用軟件FLUENT的LES(大渦模擬)計算教學的圖4

下面就可以開始進行定常計算。一般按照默認設置就OK,我這里給的定常迭代次數(shù)是100步,計算完成后,我們就有一個好的初始流場,然后關鍵的兩步來了

首先,把General中的“Steady”換成“Transient”,然后在FLUENT里面輸入

“(rpsetvar 'les-2d? #t)”

這里的問題是:注意開頭木有空格!開頭木有空格!開頭木有空格!否則FLUENT不能識別,接下來按回車,再點擊Models中的Viscious就可以看到LES,好開心有木有。

[案例分析]基于商用軟件FLUENT的LES(大渦模擬)計算教學的圖5

這是第一個需要注意的問題,因為LES默認是3D的。

第二個需要注意的問題是:為了能夠看到流場隨時間的變化,這里一定要把非定常計算過程中流場給保存下來,這是在Caculation activities中設置,

[案例分析]基于商用軟件FLUENT的LES(大渦模擬)計算教學的圖6

點擊Creat中的Solution Data Export,然后按Flow Time(物理時間)進行文件保存,這樣每一個FlowTime的流場都能保存下來啦!

最后把每個時刻的流場進行處理沒就能得到每個時刻的非定常流場啦!

  • 后處理可以使用Tecplot得到動畫,這方面不作展開啦!

有問題咨詢或者有商業(yè)合作可以聯(lián)系文主,聯(lián)系郵箱為1403476427@mail.nwpu.edu.cn。本文轉(zhuǎn)自 知乎專欄《技術鄰CAE學院》,原帖地址:https://zhuanlan.zhihu.com/p/35832439,感謝原作者。

登錄后免費查看全文
立即登錄
App下載
技術鄰APP
工程師必備
  • 項目客服
  • 培訓客服
  • 平臺客服

TOP

17
2
2