使用高格式數(shù)值模擬腔內(nèi)亞跨聲速流
摘要：本文把一階迎風(fēng)有限體積格式的高算法重構(gòu)格式——高的迎風(fēng)有限體積格式(GUVS)與不可壓縮流和可壓縮流SIMPLE 算法相結(jié)合，數(shù)值模擬了各種寬深比的亞跨聲速腔內(nèi)流動(dòng)。利用六階高迎風(fēng)有限體積格式(GUVS)在結(jié)構(gòu)和非結(jié)構(gòu)同位網(wǎng)格上模擬不同寬深比亞跨聲速腔體流動(dòng)，獲得較好的數(shù)值結(jié)果。例如亞聲速方腔流動(dòng)的數(shù)值結(jié)果與Ghia 等用多重網(wǎng)格技術(shù)獲得的Benchmark 解吻合得非常好;對(duì)不同寬深比腔內(nèi)流動(dòng)，數(shù)值模擬了大尺度旋渦、腔角小尺度旋渦的形成、分岔、破碎和邊界層分離等演化情況，表明六階精度GUVS 具有分辨率高、穩(wěn)定性好的良好性能。
關(guān)鍵詞：計(jì)算流體力學(xué);腔內(nèi)流動(dòng);有限體積格式
0、引言
腔體流動(dòng)是存在于自然對(duì)流傳熱、內(nèi)埋式武器艙、浮橋、壁面空腔設(shè)計(jì)等領(lǐng)域中非常廣泛而典型的現(xiàn)象。因此應(yīng)用計(jì)算流體力學(xué)方法數(shù)值模擬腔體流動(dòng)具有非常重要的意義。計(jì)算流體力學(xué)方法通常分為基于壓力和基于密度兩種算法。SIMPLE 系列算法是經(jīng)典的基于壓力的'算法，經(jīng)過(guò)科學(xué)家多年的研究和發(fā)展，SIMPLE 系列算法已廣泛用于不可壓縮和可壓縮各種流動(dòng)。SIMPLE 算法中對(duì)流項(xiàng)的離散通常應(yīng)用迎風(fēng)或者中心格式，通常情況下，任何高于二階的迎風(fēng)格式都是條件穩(wěn)定且在網(wǎng)格較粗時(shí)容易振蕩和發(fā)散，另外，任意高階迎風(fēng)格式均需要多于二階中心有限體積格式的節(jié)點(diǎn)數(shù)。高智提出了有限體積數(shù)值攝動(dòng)算法(NumericalPerturbational Algorithm)，由此導(dǎo)出任意階精度高迎風(fēng)有限體積格式(GUVS)。具有六階精度的高迎風(fēng)有限體積格式，在任何雷諾數(shù)下均穩(wěn)定且對(duì)粗網(wǎng)格也不振蕩和發(fā)散。而且只需要和二階中心有限體積格式相同多的節(jié)點(diǎn)，另外，高迎風(fēng)有限體積格式(GUVS)適用于任意網(wǎng)格，且變步長(zhǎng)與等步長(zhǎng)的GUVS 公式完全一致。這是高迎風(fēng)有限體積格式對(duì)比其他高階迎風(fēng)有限體積格式最突出的優(yōu)點(diǎn)。
本文將高迎風(fēng)有限體積格式與SIMPLE 系列算法相結(jié)合，數(shù)值模擬不可壓縮和可壓縮各種腔體流動(dòng)。
1、對(duì)流擴(kuò)散積分方程高迎風(fēng)格式的具體構(gòu)造。
2、高迎風(fēng)有限體積格式(GUVS)與基于壓力的不可壓和可壓縮SIMPLE 算法相耦合�；趬毫λ惴ㄊ谴蠖鄶�(shù)商用CFD 軟件的基礎(chǔ)，其中壓力是主變量，需要專門求解關(guān)于壓力的方程，密度則由狀態(tài)方程求出，特別適合不可壓縮流動(dòng)和低馬赫數(shù)流動(dòng)。壓力修正是基于壓力的修正速度場(chǎng)和壓力場(chǎng)的分步算法，Chorin 和Temamzuo 作了開(kāi)拓性的理論工作。
壓力修正在有限體積中的具體應(yīng)用就是Pantankar 著名的SIMPLE(Semi-Implicit PressureLinked Equations)系列。后來(lái),Karki 和Pantankar 發(fā)展了可壓縮SIMPLER 方法，廣泛用于各種流動(dòng)。
2.1 不可壓縮 SIMPLE 算法計(jì)算步驟簡(jiǎn)介如下①假定一個(gè)速度分布，記為u0、v0，以此計(jì)算動(dòng)量離散方程中的系數(shù)及常數(shù)項(xiàng);②假定一個(gè)壓力場(chǎng)p*;③依次求解兩個(gè)動(dòng)量方程，得到u *，v*;④求解壓力修正值方程，得到p ';⑤根據(jù) p '改進(jìn)速度值;⑥利用改進(jìn)后的速度場(chǎng)求解那些通過(guò)源項(xiàng)、物性等與速度場(chǎng)耦合的φ 變量，如果φ 與動(dòng)量方程不耦合，即不影響流場(chǎng)，則應(yīng)在速度場(chǎng)收斂后再求解;⑦重復(fù)③~⑥步，直到所有方程收斂。
在 SIMPLE 基礎(chǔ)上衍生、改進(jìn)的系列算法，與上述步驟類似。下面具體介紹GUVS 與不可壓縮SIMPLE 算法的耦合及求解。
本文使用200× 200結(jié)構(gòu)網(wǎng)格和六階精度高迎風(fēng)格式(6GUVS)計(jì)算了Re =1000的頂蓋驅(qū)動(dòng)方腔流動(dòng)，為比較起見(jiàn)，也使用200× 200結(jié)構(gòu)網(wǎng)格和一階迎風(fēng)格式(1US)計(jì)算了同一流動(dòng)。并把數(shù)值結(jié)果與Ghia 的解進(jìn)行了比較、分析。圖3 和圖4 是Benchmark 解和1US、6GUVS 的計(jì)算結(jié)果。從圖中可以看出6GUS 的精度明顯高于1US 的精度。
作為應(yīng)用，本文利用結(jié)構(gòu)網(wǎng)格和六階精度高迎風(fēng)格式(6GUVS)計(jì)算了Re =1000長(zhǎng)寬比分別為2:1和1: 2的頂蓋驅(qū)動(dòng)方腔流動(dòng)，圖5 和圖6 是長(zhǎng)寬比分別為2:1和1: 2頂蓋驅(qū)動(dòng)方腔流動(dòng)的流線圖。
2.2 可壓縮 SIMPLE 算法簡(jiǎn)述如下①假定一個(gè)速度和密度分布，記為u0、v0、ρ 0 以此計(jì)算動(dòng)量離散方程中的系數(shù)及常數(shù)項(xiàng);②假定一個(gè)壓力場(chǎng)p*;③依次求解兩個(gè)動(dòng)量方程，得到u *，v*;④求解壓力修正值方程，得到p ';⑤根據(jù) p '改進(jìn)速度場(chǎng)、壓力場(chǎng)和密度場(chǎng);⑥求解能量方程，更新密度場(chǎng)⑦重復(fù)③~⑥步，直到所有方程收斂。
下面具體介紹高迎風(fēng)格式與可壓縮SIMPLE 算法的耦合及求解。
3、結(jié)論
將高迎風(fēng)有限體積格式(GUVS)與SIMPLE 系列算法相結(jié)合，數(shù)值模擬不可壓縮和可壓縮不同寬深比腔內(nèi)流動(dòng),獲得較好的數(shù)值結(jié)果,表明六階精度GUVS 具有分辨率高、穩(wěn)定性好的良好性能。繼續(xù)將該格式應(yīng)用于復(fù)雜流場(chǎng)的數(shù)值模擬是下一步將進(jìn)行的工作。
[參考文獻(xiàn)] (References)
[1] 高智. 對(duì)流擴(kuò)散方程的攝動(dòng)有限體積方法[A]，第十一屆全國(guó)計(jì)算流體力學(xué)會(huì)議論文集，洛陽(yáng)，2002，38～45.
[2] 高智，柏威. 對(duì)流擴(kuò)散方程的攝動(dòng)有限體積(PFV)方法及討論[J]，力學(xué)學(xué)報(bào)，2004，36(1)：88-9.
[3] 高智，向華，申義慶. 攝動(dòng)有限體積法重構(gòu)近似高精度的意義[J]，計(jì)算物理，2004，21(2)：131-136.
[4] 代民果，高智. 同位網(wǎng)格攝動(dòng)有限體積格式求解浮力驅(qū)動(dòng)方腔流[J]，力學(xué)學(xué)報(bào)，2006，38(6)：733-740.
[5] 代民果. 有限體積迎風(fēng)高格式及基于非結(jié)構(gòu)混合網(wǎng)格的復(fù)雜流動(dòng)數(shù)值模擬[D]，北京：中國(guó)科學(xué)院研究生部，2008.
【使用高格式數(shù)值模擬腔內(nèi)亞跨聲速流】相關(guān)文章：
1.熔鋁爐內(nèi)交變磁場(chǎng)對(duì)于電磁力場(chǎng)的數(shù)值模擬研究
2.臨床執(zhí)業(yè)醫(yī)師助理考點(diǎn)：胸膜腔和胸膜腔內(nèi)壓
3.基于FDS的沙發(fā)火災(zāi)數(shù)值模擬