流固熱耦合仿真的實(shí)例教程
FLUENT流-固-熱耦合分析
ANSYS FLUENT軟件自V2019版本起,新增了Structure結(jié)構(gòu)求解功能,能夠基于Fluent軟件進(jìn)行簡(jiǎn)單模型的結(jié)構(gòu)應(yīng)力、變形分析,具備線性及非線性結(jié)構(gòu)分析功能。本案例基于ANSYS FLUENT 2020R1進(jìn)行管道閥門(mén)流-固-熱三場(chǎng)耦合分析。
1 模型描述
如圖所示尺寸的三維管道模型,管道模型中存在4個(gè)簡(jiǎn)化的閥瓣模型,給定管道入口氣體流速為10m/s,閥板內(nèi)給定體積熱源為2000000w/m^3;
閥瓣模型材料參數(shù):
密度:2700kg/m^3;
比熱:871J/kg.K;
熱傳導(dǎo)系數(shù):202W/m^2.K;
楊氏模量:2.5E7Pa;
泊松比:0.37;
2 網(wǎng)格劃分
本案例網(wǎng)格基于ANSYS ICEM CFD進(jìn)行全六面體網(wǎng)格劃分,網(wǎng)格如下圖所示:
流體區(qū)域:480000六面體網(wǎng)格;
固體區(qū)域:3800六面體網(wǎng)格。
3 FLUENT求解設(shè)置
求解計(jì)算分兩步完成,首先不考慮結(jié)構(gòu)變形對(duì)流體-固體進(jìn)行穩(wěn)態(tài)共軛傳熱分析,然后基于上一步仿真計(jì)算結(jié)果考慮流固耦合作用實(shí)現(xiàn)瞬態(tài)流-固-熱耦合仿真分析。
3.1流固共軛傳熱仿真
? 啟動(dòng)FLUENT軟件,利用菜單File>>Read case….打開(kāi)文件對(duì)話框,讀入網(wǎng)格文件vavle_test.msh;新版本顯式界面如下:
? 新版本的FLUENT軟件默認(rèn)選擇k-w sst湍流模型,本案例不做修改;
? 激活能量方程
? 邊界條件設(shè)置
1)固體區(qū)域熱源:2000000W/m^3;選擇對(duì)應(yīng)的固體區(qū)域,勾選source terms加載能量源項(xiàng)。
展開(kāi) 軸承和軸在仿真中的意義也不明顯,因此我們都予以簡(jiǎn)化。
分析傳熱模型,齒輪摩擦生熱是熱源,這些熱量通過(guò)幾種方式傳播:
1.熱傳導(dǎo)——從齒緣往齒輪中心傳導(dǎo)
2.熱對(duì)流——齒輪和潤(rùn)滑油,潤(rùn)滑油和空氣,又稱(chēng)為共軛傳熱
3.熱輻射——溫度不高,輻射量小可忽略
因此,滑油和空氣是傳熱的介質(zhì),必須在模型中考慮進(jìn)去(事實(shí)上這部分傳熱達(dá)到91%)。滑油和空氣是兩相,因此要使用到fluent的多相流模型;要模擬甩油過(guò)程,要使用動(dòng)網(wǎng)格模型;要模擬傳熱過(guò)程,利用fluent內(nèi)建的傳熱模型。這三者是本案例的核心。
這里不得不提到兩位外國(guó)學(xué)者,Guillaume Houzeaux對(duì)齒輪泵進(jìn)行了仿真,并且關(guān)注局部網(wǎng)格,這可能是最早對(duì)齒輪+流體進(jìn)行仿真;而F.Lemfeld率先采用兩相流模型捕捉了齒輪箱內(nèi)的流體瞬態(tài)變化情況,但他在網(wǎng)格方面的處理比較簡(jiǎn)單,對(duì)齒輪齒形進(jìn)行了切除,同時(shí)使用一定的壁面粗糙度值模擬齒形的存在,使齒輪能夠甩油。
說(shuō)了這么多廢話,現(xiàn)在回到主題。
圖3 流固熱耦合仿真流程
本例需要用到的模塊包括fluent模塊,其中又集成了ansys自帶的幾何處理與網(wǎng)格劃分工具。后面與fluent共享結(jié)果的是穩(wěn)態(tài)熱分析模塊,以及靜力結(jié)構(gòu)模塊,用來(lái)分析熱應(yīng)力對(duì)結(jié)構(gòu)的影響,如用來(lái)分析熱變形,限于篇幅本例不涉及。本例實(shí)際流程可以簡(jiǎn)化如下,我個(gè)人喜歡拆分不同的模塊,這樣方便“故障隔離”:
圖4 流體仿真流程
一、模型簡(jiǎn)化與網(wǎng)格劃分
由于復(fù)雜的三維結(jié)構(gòu)會(huì)增加網(wǎng)格劃分的難度,會(huì)導(dǎo)致網(wǎng)格數(shù)目的無(wú)謂增加,加大計(jì)算量,因此對(duì)齒輪減速器三維模型進(jìn)行簡(jiǎn)化:殼體的凸臺(tái)、通孔、墊圈等予以去除;統(tǒng)一壁面厚度;滾動(dòng)軸承結(jié)構(gòu)在對(duì)應(yīng)位置采取同心圓環(huán)來(lái)表示,方便施加熱流。這里的模型簡(jiǎn)化工作是用SpaceClaim做的。
展開(kāi) 隨著計(jì)算流體力學(xué)的發(fā)展以及計(jì)算性能的提升,對(duì)航空發(fā)動(dòng)機(jī)整機(jī)仿真成為了可能,本教程對(duì)KJ66航空發(fā)動(dòng)機(jī)進(jìn)行整機(jī)仿真,整機(jī)仿真結(jié)合氣動(dòng)、傳熱、燃燒、多相流、固體應(yīng)力,將航空發(fā)動(dòng)機(jī)從冷態(tài)計(jì)算至熱態(tài),即仿真始于冷態(tài),終于熱態(tài)。
KJ66航空發(fā)動(dòng)機(jī)幾何模型如圖,對(duì)航空發(fā)動(dòng)機(jī)氣熱彈耦合仿真,計(jì)算采用穩(wěn)態(tài),氣動(dòng)的計(jì)算采用求解粘性N-S方程的方法,燃油的噴射計(jì)算采用拉格朗日多相流,燃燒的計(jì)算采用有限速率的渦耗散模型,流體與結(jié)構(gòu)的相互作用(FSI)采用雙向耦合的方式。
流體結(jié)構(gòu)相互作用 (FSI)是指一種耦合的表面問(wèn)題,其中流體模型的狀態(tài)取決于結(jié)構(gòu)模型的狀態(tài),反之亦然。這種相互關(guān)系可以是對(duì)稱(chēng)或非對(duì)稱(chēng)的。非對(duì)稱(chēng)問(wèn)題通常指單向耦合問(wèn)題,表示其中一個(gè)模型是獨(dú)立的,另一個(gè)模型則具有關(guān)聯(lián)性。
流體結(jié)構(gòu)相互作用(FSI)耦合交界面處的對(duì)應(yīng)流體和固體移動(dòng)時(shí)運(yùn)動(dòng)學(xué)特性(位置、速度和加速度)相同,受到的力也相同。
從流體傳遞到固體的信息是流體拉力,它由流體壓力和壁面剪切應(yīng)力組成的。此傳遞發(fā)生在耦合壁面邊界流體-結(jié)構(gòu)交界面)上。
從固體傳遞到流體的信息是固體的變形,尤其是流體-結(jié)構(gòu)交界面的變形。
一般情況下,F(xiàn)SI模擬在運(yùn)動(dòng)學(xué)和力方面保持一致,稱(chēng)為雙向耦合,在STAR-CCM+中,雙向耦合FSI問(wèn)題是指從流體到固體和從固體到流體的交換的綜合采用并行求解方法。
進(jìn)行航空發(fā)動(dòng)機(jī)整機(jī)氣熱彈耦合仿真的STAR-CCM+版本為STAR-CCM+ 2206.
將航空發(fā)動(dòng)機(jī)整機(jī)從冷態(tài)模型計(jì)算至熱態(tài)模型后發(fā)動(dòng)機(jī)伸長(zhǎng)約1mm。
詳細(xì)計(jì)算結(jié)果如下:
速度
溫度
溫度
位移
固體應(yīng)力
文章來(lái)源:STAR CCM仿真學(xué)堂
展開(kāi) 為了更好地了解風(fēng)機(jī)的結(jié)構(gòu)及特點(diǎn),提高風(fēng)機(jī)的總體設(shè)計(jì)水平與使用效能,可通過(guò)自建高性能并行集群仿真平臺(tái), 利用OpenFOAM開(kāi)源軟件進(jìn)行計(jì)算, 考慮流固耦合方式對(duì)風(fēng)機(jī)葉片上的氣動(dòng)載荷進(jìn)行分析。 下圖為數(shù)值模擬結(jié)果。
風(fēng)機(jī)在計(jì)算域中的示意圖
風(fēng)機(jī)在計(jì)算域中的示意圖
風(fēng)機(jī)在簡(jiǎn)化氣動(dòng)力下轉(zhuǎn)動(dòng)效果
流固耦合條件下模擬,可以考慮風(fēng)機(jī)塔架、機(jī)艙的振動(dòng)響應(yīng)。
在此種模擬方法下,可以輸出風(fēng)場(chǎng)縱剖面速度云圖,考慮風(fēng)機(jī)的尾流效應(yīng)。
單風(fēng)機(jī)尾渦效果展示
雙風(fēng)機(jī)尾渦效果展示
葉片是風(fēng)力發(fā)電機(jī)中最基礎(chǔ)和最關(guān)鍵的部件,其良好的設(shè)計(jì),可靠的質(zhì)量和優(yōu)越的性能是保證機(jī)組正常穩(wěn)定運(yùn)行的決定因素。考慮流固耦合方式對(duì)風(fēng)機(jī)葉片上的氣動(dòng)載荷進(jìn)行分析,可以為風(fēng)機(jī)的總體設(shè)計(jì)提供一個(gè)較為全面的建議及分析方法。
展開(kāi) TASK
熱流固耦合分析中包括沸騰模型修正、接觸熱阻 計(jì)算和位移的傅立葉分解這三項(xiàng)功能。由于模型中部分區(qū)域發(fā)生了沸騰,非沸騰態(tài)下的換熱公式不再適用于計(jì)算沸騰態(tài)下的換熱量,因此需要對(duì)模型的換熱系數(shù)進(jìn)行修正;接觸熱阻程序?qū)崿F(xiàn)的功能是根據(jù)接觸面之間的實(shí)際接觸面積、接觸表面的材料、接觸面間隙中介質(zhì)的導(dǎo)熱系數(shù)和接觸面的壓力計(jì)算接觸面的接觸熱阻;對(duì)于已知的位移結(jié)果,可以在二維坐標(biāo)系下可以將該平面內(nèi)的位移進(jìn)行傅立葉分解,展開(kāi)成多階傅立葉級(jí)數(shù)。
SOLUTION
主要技術(shù)挑戰(zhàn):
沸騰修正涉及結(jié)構(gòu)模型和流體模型之間網(wǎng)格的插值和數(shù)據(jù)傳遞;
接觸熱阻公式較復(fù)雜,涉及物理量較多;
位移傅立葉分解計(jì)算較復(fù)雜;
解決方案:
開(kāi)發(fā)沸騰修正模板,實(shí)現(xiàn)插值和模型修正功能;
開(kāi)發(fā)接觸熱阻模板,實(shí)現(xiàn)熱阻公式的計(jì)算;
開(kāi)發(fā)位移傅立葉分解模板,實(shí)現(xiàn)位移的傅立葉分解,并合并各階結(jié)果;
提供豐富的參數(shù)輸入和輸出界面;
結(jié)論:
形成了完整的熱流固耦合分析模板;
模板包括了沸騰修正、接觸熱阻和位移傅立葉分解功能。
Customer Benefit
熱流固耦合分析模板搭建的流程包含了沸騰修正、接觸熱阻和位移傅立葉分解的功能,已經(jīng)直接集成在柴油發(fā)動(dòng)機(jī)仿真分析模板系統(tǒng)中,成為了柴油機(jī)整體仿真方案的一部分。
本文來(lái)自安世亞太微信公號(hào),如果您對(duì)耦合分析有需求或感興趣,歡迎聯(lián)系溝通:
400-6600-388
展開(kāi)