【摘要】：計(jì)算流體力學(xué)的持續(xù)發(fā)展對(duì)算法的求解精度提出了更高的要求,但是目前多數(shù)流場(chǎng)求解器都只能達(dá)到二階精度,難以滿足CFD工作者解決復(fù)雜流體工程問(wèn)題和更深入地研究流動(dòng)本質(zhì)的需要。在此背景下,間斷伽遼金有限元方法作為一種高精度數(shù)值方法近年來(lái)受到廣泛關(guān)注。它通過(guò)增加單元內(nèi)解函數(shù)多項(xiàng)式項(xiàng)數(shù)和階數(shù)的方法提高求解精度。該方法兼具傳統(tǒng)有限元法和有限體積法的優(yōu)點(diǎn),并且融合了數(shù)值通量、近似黎曼解、TVD和限制器等內(nèi)容。相比較于其它幾種計(jì)算流體力學(xué)傳統(tǒng)方法,間斷伽遼金有限元法的優(yōu)勢(shì)相當(dāng)顯著。鑒于間斷伽遼金有限元法的巨大優(yōu)勢(shì)和課題組研究的需要,本文開(kāi)展了龍格庫(kù)塔間斷伽遼金有限元法(RKDGM)的復(fù)現(xiàn)工作,驗(yàn)證該方法的計(jì)算表現(xiàn),并摸索方法中時(shí)間推進(jìn)方法、通量分裂格式、限制器、不同類型網(wǎng)格的具體實(shí)現(xiàn)。研究基于歐拉方程進(jìn)行,首先詳細(xì)闡述了一維和二維下RKDGM以及多種通量分裂格式、限制器的理論框架,并詳細(xì)介紹了不同網(wǎng)格類型對(duì)應(yīng)的基函數(shù)、坐標(biāo)變換和數(shù)值求積方法的選擇。然后編寫了一維RKDGM求解器,在探究通量分裂格式、限制器等算法模塊具體程序?qū)崿F(xiàn)細(xì)節(jié)的同時(shí),探究了基函數(shù)初始系數(shù)、格式空間離散精度、限制器對(duì)計(jì)算結(jié)果的影響,并得出一些結(jié)論。以Shu激波管問(wèn)題和光滑流動(dòng)問(wèn)題測(cè)試了基函數(shù)初始系數(shù)對(duì)計(jì)算精度的影響。通過(guò)比較分析發(fā)現(xiàn):即使初始守恒變量初始分布是精確的,基函數(shù)初始系數(shù)的不精確也會(huì)導(dǎo)致計(jì)算精度的損失,且該問(wèn)題在光滑流動(dòng)問(wèn)題中更為突出。進(jìn)而提出一種可以保證基函數(shù)初始系數(shù)精確性的賦值方法,且對(duì)計(jì)算精度的改善效果十分顯著。以激波管問(wèn)題為算例測(cè)試了提高格式精度和加密網(wǎng)格對(duì)計(jì)算精度的改善程度。數(shù)值實(shí)驗(yàn)結(jié)果顯示,兩種方法對(duì)計(jì)算結(jié)果都有改善作用,其中提高格式精度對(duì)計(jì)算精度的改善更為突出。但是這個(gè)改善是有限度的,當(dāng)格式精度提高到一定階數(shù),網(wǎng)格達(dá)到一定密度,進(jìn)一步提高得空間極小。出于節(jié)約計(jì)算資源和提高計(jì)算精度的綜合考慮,給出了非必須情況下,計(jì)算時(shí)應(yīng)選擇稍高格式階數(shù)和稍小網(wǎng)格量的建議。以Lax激波管問(wèn)題為算例測(cè)試minmod限制器和矩限制器的性能表現(xiàn)。數(shù)值實(shí)驗(yàn)結(jié)果顯示,minmod限制器耗散較小,容易出現(xiàn)振蕩;矩限制器耗散相對(duì)較大,對(duì)間斷有些許抹平。基于該發(fā)現(xiàn),提出間斷探測(cè)器與矩限制器搭配使用的策略,通過(guò)設(shè)置適當(dāng)?shù)奶綔y(cè)閾值,有效提高了矩限制器的性能表現(xiàn)。
【圖文】：

5.4 四邊形域的坐標(biāo)變換與數(shù)值積分將真實(shí)物理空間( x, y)內(nèi)的矩形域e變換到積分空間( , )中的標(biāo)準(zhǔn)參考單元 E 如圖5.2 示。圖 5.2 四邊形域坐標(biāo)變換示意圖Fig. 5.2 The illustration of the coordinate transformation of a quadrilateral domaing123feabdcyx231ex231E44

5.4 四邊形域的坐標(biāo)變換與數(shù)值積分將真實(shí)物理空間( x, y)內(nèi)的矩形域e變換到積分空間( , )中的標(biāo)準(zhǔn)參考單元 E 如圖5.2 示。圖 5.2 四邊形域坐標(biāo)變換示意圖Fig. 5.2 The illustration of the coordinate transformation of a quadrilateral domaing123feabdcyx231ex231E44
【學(xué)位授予單位】：大連理工大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2018
【分類號(hào)】：O241.82;O35

【參考文獻(xiàn)】

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

1 孫強(qiáng);呂宏強(qiáng);伍貽兆;;基于高階物面近似的自適應(yīng)間斷有限元法歐拉方程數(shù)值模擬[J];空氣動(dòng)力學(xué)學(xué)報(bào);2015年04期

2 劉君;鄒東陽(yáng);徐春光;;基于非結(jié)構(gòu)動(dòng)網(wǎng)格的非定常激波裝配法[J];空氣動(dòng)力學(xué)學(xué)報(bào);2015年01期

3 張來(lái)平;李明;劉偉;赫新;張涵信;;基于非結(jié)構(gòu)/混合網(wǎng)格的高階精度DG/FV混合方法研究進(jìn)展[J];空氣動(dòng)力學(xué)學(xué)報(bào);2014年06期

4 秦望龍;呂宏強(qiáng);伍貽兆;;彎曲網(wǎng)格上的間斷有限元湍流數(shù)值解法研究[J];空氣動(dòng)力學(xué)學(xué)報(bào);2014年05期

5 李珍珍;蔚喜軍;趙國(guó)忠;馮濤;;RKDG有限元法求解一維拉格朗日形式的Euler方程[J];計(jì)算物理;2014年01期

6 秦望龍;呂宏強(qiáng);伍貽兆;;基于混合網(wǎng)格的高階間斷有限元黏流數(shù)值解法[J];力學(xué)學(xué)報(bào);2013年06期

7 劉偉;張來(lái)平;赫新;賀立新;張涵信;;基于Newton/Gauss-Seidel迭代的DGM隱式方法[J];力學(xué)學(xué)報(bào);2012年04期

8 李喜樂(lè);楊永;郝海兵;焦瑾;;求解RANS方程的高階間斷Galerkin方法研究[J];西北工業(yè)大學(xué)學(xué)報(bào);2012年03期

9 夏軼棟;伍貽兆;呂宏強(qiáng);宋江勇;;高階間斷有限元法的并行計(jì)算研究[J];空氣動(dòng)力學(xué)學(xué)報(bào);2011年05期

10 張來(lái)平;劉偉;賀立新;赫新;鄧小剛;;基于靜動(dòng)態(tài)重構(gòu)的高階DG/FV混合格式在二維非結(jié)構(gòu)網(wǎng)格中的推廣[J];計(jì)算物理;2011年02期

相關(guān)博士學(xué)位論文 前1條

1 李明;基于混合網(wǎng)格的高階精度DG/FV混合算法研究[D];中國(guó)空氣動(dòng)力研究與發(fā)展中心;2013年

相關(guān)碩士學(xué)位論文 前1條

1 劉偉;基于混合網(wǎng)格的高階間斷Galerkin/有限體積混合算法的研究[D];中國(guó)空氣動(dòng)力研究與發(fā)展中心;2010年

，

本文編號(hào)：2600808