建立一種修正的HLLC(Harten-Lax-Van Leer-Contact)格式下穩(wěn)定的Mie-Grüneisen多介質(zhì)混合計(jì)算模型.Mie-Grüneisen混合模型中的通量包括守恒和非守恒兩個(gè)部分,原始的HLLC格式對(duì)守恒部分適用,但是原始的HLLC格式直接用于非守恒部分,很難控制住數(shù)值振蕩的產(chǎn)生.在原始格式中,間斷面的移動(dòng)速度為間斷網(wǎng)格的左側(cè)或右側(cè)速度,修正后替換為網(wǎng)格中的平均速度,經(jīng)過修正后,對(duì)HLLC格式重新進(jìn)行推導(dǎo),并隨之?dāng)U展到二維問題.數(shù)值實(shí)驗(yàn)表明,利用修改后的HLLC格式,Mie-Grüneisen混合模型可以取得較好的穩(wěn)定性和準(zhǔn)確性.

【文章頁數(shù)】：8 頁

【部分圖文】：

圖3120μs時(shí),鉬和MORB相互作用的p、u、Γ計(jì)算結(jié)果

該算例中,壓縮態(tài)的鉬與正常態(tài)的鉬先在0.4m處發(fā)生相互作用,再與MORB發(fā)生相互作用.由于情況復(fù)雜,文獻(xiàn)[11]僅提供了數(shù)值解,未提供準(zhǔn)確解.計(jì)算100網(wǎng)格下的HLLC的數(shù)值解,選取文獻(xiàn)[9]的計(jì)算結(jié)果進(jìn)行對(duì)比,如圖3所示.從圖3中可以看出,未修改的HLLC方法,從0.6m開....

圖4T=50μs,(a)密度,(b)壓力和(c)體積分?jǐn)?shù)分布

該問題在文獻(xiàn)[9]中同樣關(guān)注過,且計(jì)算結(jié)果大體相似.鉬與MORB的作用機(jī)制與上一算例基本相同,受到鉬沖擊波的沖擊后,部分進(jìn)入MORB內(nèi)部進(jìn)行擴(kuò)散,部分反射回鉬中形成稀疏波.圖4為50μs時(shí)刻的云圖,從圖中可以看到,由于MORB為正方形,進(jìn)入MORB中繼續(xù)傳播的激波,以強(qiáng)間斷形式....

圖5T=100μs,(a)密度,(b)壓力和(c)體積分?jǐn)?shù)分布

圖4T=50μs,(a)密度,(b)壓力和(c)體積分?jǐn)?shù)分布3結(jié)論

圖1網(wǎng)格處的左右間斷面

對(duì)于存在m種流體介質(zhì)的流場(chǎng),假設(shè)第I種介質(zhì)的體積分?jǐn)?shù)為ZI,那么體積分?jǐn)?shù)滿足輸送方程,即方程(1)中的最后一個(gè)方程.1.2原始HLLC計(jì)算方法

本文編號(hào)：3908568