針對數(shù)值求解滲流方程時,使用標(biāo)準(zhǔn)有限體積法出現(xiàn)數(shù)值界面不能有效向前傳播的"數(shù)值熱障"現(xiàn)象,提出一種修正的有限體積法,該方法擴(kuò)散系數(shù)的取值采用密度變量在兩個相鄰單元的代數(shù)平均值.數(shù)值實驗結(jié)果表明,新格式可有效避免"數(shù)值熱障"現(xiàn)象.

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

【部分圖文】：

圖2中心型有限體積法的數(shù)值結(jié)果

圖1當(dāng)t=0.1,1時問題（1）的Barenblatt解2修正的單元中心有限體積法

圖3修正的中心型有限體積法數(shù)值結(jié)果

考慮Barenblatt解問題,使用上述推導(dǎo)的修正的中心型有限體積法求解,數(shù)值結(jié)果如圖3所示.由圖3可見,新方法不僅沒有“數(shù)值熱障”現(xiàn)象,還可以很好地捕捉波陣面且數(shù)值解保正.表1～表3列出了本文提出的新格式與文獻(xiàn)[6-8]中的算法在求解Barenblatt問題時數(shù)值結(jié)果....

圖1當(dāng)t=0.1,1時問題（1）的Barenblatt解

其中u+=max{0,u}.因此初值條件可取為g(x)=Bm(x,t0).圖1為Barenblatt解的圖像.使用1.1中推導(dǎo)的中心型有限體積法求解該問題的數(shù)值結(jié)果如圖2所示.由圖2可見,計算終止時刻數(shù)值波陣面未離開初始位置,產(chǎn)生了數(shù)值振蕩.因此,標(biāo)準(zhǔn)單元中心型....

本文編號：3974052