Runge-Kutta間斷Galerkin(RKDG)方法是數(shù)值求解雙曲守恒律方程的主流方法之一,它具有精度高、易于處理復(fù)雜區(qū)域和邊界條件、易于自適應(yīng)和并行效率高等眾多優(yōu)點(diǎn)。很多雙曲守恒律問題的規(guī)模特別巨大,數(shù)值計算求解時需要耗費(fèi)大量的時間和計算機(jī)存儲空間。這讓自適應(yīng)方法成為RKDG方法的一個重要研究方向。我們首先設(shè)計了一個基于壞單元指示子的p自適應(yīng)RKDG算法。這個算法在每一時間步利用壞單元指示子探測并標(biāo)記壞單元,然后對壞單元使用低次逼近多項(xiàng)式,而其它單元使用高次逼近多項(xiàng)式,由此省掉間斷附近不必要的自由度,節(jié)省計算存儲量。我們將這個算法應(yīng)用到了一維Burgers問題、Riemann問題和Buckley-Leverett問題上,得到的數(shù)值結(jié)果表明該算法不僅能夠保持原RKDG方法在光滑區(qū)域的高次逼近,而且能夠有效降低存儲空間。隨后我們在p自適應(yīng)算法的基礎(chǔ)上,通過對壞單元進(jìn)行加密、對被加密的好單元對進(jìn)行合并,又加入了網(wǎng)格的h自適應(yīng),形成了一種求解一維雙曲守恒律方程的hp自適應(yīng)RKDG方法。數(shù)值測試表明新的算法既有p自適應(yīng)算法節(jié)省計算存儲的優(yōu)點(diǎn),還有h自適應(yīng)算法提高間斷處模擬效果的優(yōu)勢。跟同等單元個數(shù)的普通RKDG方法相比,兩個方法都能有效控制數(shù)值振蕩,但hp自適應(yīng)方法得到的數(shù)值解更加優(yōu)質(zhì)。最后,我們嘗試對hp自適應(yīng)RKDG算法進(jìn)行改進(jìn),引入了Q(Q1)個時間步只進(jìn)行一次自適應(yīng)的做法。數(shù)值試驗(yàn)表明,在間隔步數(shù)Q取得比較恰當(dāng)?shù)那闆r下,改進(jìn)的算法在不明顯降低數(shù)值解質(zhì)量的同時,能夠節(jié)省一定的計算時間。
【學(xué)位單位】：南京郵電大學(xué)
【學(xué)位級別】：碩士
【學(xué)位年份】：2018
【中圖分類】：O241.82

【參考文獻(xiàn)】

相關(guān)期刊論文 前10條

1 倪建;朱君;;帶TVB限制器的RKDG方法與浸入邊界方法在Euler方程中的應(yīng)用[J];南京信息工程大學(xué)學(xué)報(自然科學(xué)版);2014年01期

2 陳浩;蘇劍;王尚錦;朱明雷;;雙曲型守恒律方程的自適應(yīng)間斷Galerkin方法[J];工程數(shù)學(xué)學(xué)報;2013年02期

3 徐云;蔚喜軍;;自適應(yīng)間斷有限元方法求解雙曲守恒律方程[J];計算物理;2009年02期

4 金杉,莊達(dá)民,張向陽;計算流體力學(xué)在現(xiàn)代建筑消防設(shè)計中的應(yīng)用[J];消防科學(xué)與技術(shù);2003年03期

5 傅曉英,劉俊,許劍峰,羅麟;計算流體力學(xué)在城市規(guī)劃設(shè)計中的應(yīng)用研究[J];四川大學(xué)學(xué)報(工程科學(xué)版);2002年06期

6 劉寶杰,鄒正平,嚴(yán)明,劉火星,寧方飛,張永新,徐力平,蔣浩康,陳懋章;葉輪機(jī)計算流體動力學(xué)技術(shù)現(xiàn)狀與發(fā)展趨勢[J];航空學(xué)報;2002年05期

7 周連第;船舶與海洋工程計算流體力學(xué)的研究進(jìn)展與應(yīng)用[J];空氣動力學(xué)學(xué)報;1998年01期

8 常思勤,扈圣剛;計算流體力學(xué)進(jìn)展及其在汽車設(shè)計中的應(yīng)用[J];武漢汽車工業(yè)大學(xué)學(xué)報;1997年04期

9 肖柯則，夏藝;計算流體力學(xué)在鑄造過程中的應(yīng)用[J];內(nèi)蒙古工業(yè)大學(xué)學(xué)報(自然科學(xué)版);1995年03期

10 吳正;低速混合型城市交通的流體力學(xué)模型[J];力學(xué)學(xué)報;1994年02期

本文編號：2826243