本文關鍵詞： 有限差分方法 高階精度差分格式 幾何守恒律 復雜構型 數(shù)值精度分析 網(wǎng)格光滑性 網(wǎng)格質量檢測　出處：《中國空氣動力研究與發(fā)展中心》2015年博士論文　論文類型：學位論文

【摘要】：作為最為古老的數(shù)值離散方法之一,有限差分方法因其簡單易懂而備受初學者的青睞,也因此而成為諸多數(shù)值離散類教科書的必講內容。然而,令許多從事數(shù)值計算尤其是流體動力學數(shù)值仿真工作的研究人員感到困惑的是：對數(shù)值離散方法具有建立基本概念和啟蒙思想作用的有限差分方法在實際工程問題的應用中卻難覓其蹤跡。在固體力學問題的計算應用中有限元方法已經(jīng)相當成熟,而幾乎所有的流體力學商業(yè)計算軟件都是基于有限體積框架的,僅有少數(shù)研究人員基于有限差分方法編寫計算代碼,并主要結合高階精度差分格式應用于繞簡單外形的流動機理研究中。目前有限差分方法難以應用于大規(guī)模的工程實際問題中,其主要原因是有限差分方法在復雜網(wǎng)格中的計算魯棒性較差。不同于有限體積方法,有限差分方法實施的前提是坐標變換(從物理空間到計算空間)及其逆變換(從計算空間到物理空間)的嚴格成立,因此有限差分方法強烈依賴于坐標變換。大量的研究工作表明：在坐標變換的諸多特性中,其幾何守恒律特性對有限差分方法的計算魯棒性具有決定性的影響。坐標變換的幾何守恒律特性在微分意義下是自動成立的,但在差分格式的離散下卻不一定成立。差分離散下幾何守恒律的難以滿足嚴重影響了有限差分方法尤其是高階精度有限差分方法的魯棒性,極大地限制了其在工程實際問題中的應用。本文詳細分析了滿足幾何守恒律的充分條件(CMM),并在CMM條件的指引下依據(jù)各變換系數(shù)離散后的幾何意義給出了其對稱守恒計算形式(SCMM),隨后還對SCMM的優(yōu)越性通過數(shù)值試驗進行了驗證。針對幾何守恒律誤差對數(shù)值離散結果的影響,’本文首先詳細分析了網(wǎng)格光滑性對幾何守恒律誤差數(shù)值收斂精度的影響規(guī)律,隨后分析了幾何守恒律的滿足與否對數(shù)值離散結果數(shù)值精度的影響,并均通過相應的數(shù)值試驗對理論分析的結論進行了驗證。本文的研究結果表明：幾何守恒律的滿足能夠有效減小數(shù)值離散誤差,而當網(wǎng)格不夠光滑時,幾何守恒律的滿足能夠將數(shù)值離散結果的精度提高整整一階。在分析雅可比的不同計算形式所代表的幾何意義時,我們發(fā)現(xiàn)不同計算形式離散后的數(shù)值體積是不一樣的,當計算網(wǎng)格單元為平行六面體時,三種計算形式的數(shù)值體積相等；而當網(wǎng)格單元不是平行六面體時,三個數(shù)值體積一般不相等。一般而言,計算網(wǎng)格單元相對于平行六面體的畸變愈大,則雅可比的三種計算形式離散后的數(shù)值體積差異便愈大,據(jù)此本文提出了基于SCMM的網(wǎng)格檢測方法,可以對結構網(wǎng)格的網(wǎng)格質量進行檢測。與現(xiàn)有商業(yè)軟件的網(wǎng)格檢測方法不同,本文的網(wǎng)格檢測方法與數(shù)值離散格式相關。即使是對于同一套網(wǎng)格,數(shù)值離散格式不同時其檢測結果也是各異的,能較為充分地體現(xiàn)網(wǎng)格質量與數(shù)值離散格式的相關性。最后,在相關研究結果的基礎上,本文初步實現(xiàn)了高階精度有限差分格式對三維復雜外形流動的數(shù)值模擬能力。全文共分為七章,各章內容概述如下：第一章為引言。首先簡要介紹了有限差分方法尤其是高階精度有限差分方法的研究歷史與現(xiàn)狀,并指出幾何守恒律問題的懸而未決直接限制了高階精度有限差分方法的推廣應用。隨后對有限差分方法中關于幾何守恒律問題的研究工作進行了簡單敘述,并分別介紹了體積守恒律和面積守恒律的國內外研究概況。最后圍繞有限差分方法的幾何守恒律問題提出了本文的主要研究內容。第二章為幾何守恒律問題的提出。首先簡要介紹了有限差分方法中的坐標變換關系,并據(jù)此將一般守恒律控制方程由直角坐標系變換到計算坐標系。在計算坐標系下給出有限差分離散算子的一般約定特性后,由自由流守恒條件導出有限差分離散中的幾何守恒律問題。第三章為幾何守恒律問題的理論研究。在簡單分析了國內外在兒何守恒律方面的各種解決途徑的優(yōu)劣后,依據(jù)滿足幾何守恒律的充分條件,提出了有限差分方法中網(wǎng)格變換導數(shù)的對稱守恒計算方法(SCMM)。并在隨機網(wǎng)格中對SCMM的有效性進行了數(shù)值驗證,基于線性和非線性問題的數(shù)值實驗均表明：本文提出的幾何守恒律解決途徑(SCMM)能有效減小有限差分方法在非光滑網(wǎng)格中的離散誤差,尤其是對于非線性問題。第四章為網(wǎng)格光滑性對幾何守恒律誤差的影響研究。對幾何守恒律不滿足時(網(wǎng)格變換導數(shù)的內、外層差分算子不相同)其數(shù)值誤差與網(wǎng)格光滑性的關系進行了理論推導,并通過數(shù)值實驗對理論推導的正確性進行了驗證。另外,對數(shù)值實驗中涉及的不同精度誤差的統(tǒng)計分析結果也給出了詳細的理論證明。第五章為幾何守恒律誤差對數(shù)值離散精度的影響研究,通過數(shù)值實驗和理論分析研究著重考察了在不同光滑程度網(wǎng)格中,幾何守恒律的滿足與否對物理解的數(shù)值離散精度的影響。幾何守恒律的滿足能有效減小物理解的數(shù)值離散誤差,而當計算網(wǎng)格不是足夠光滑時,幾何守恒律的滿足還能將物理解的數(shù)值離散精度提高整整一階。第六章為網(wǎng)格檢測方法研究。針對現(xiàn)有網(wǎng)格檢測方法與數(shù)值離散過程脫節(jié)的問題,提出了與數(shù)值離散方法相關的SCMM網(wǎng)格檢測方法,并基于坐標變換雅可比的不同離散形式之間的差異體現(xiàn)了網(wǎng)格單元與標準平行六面體之間的差異給出了SCMM網(wǎng)格檢測方法的檢測指標。隨后的數(shù)值實驗表明：SCMM網(wǎng)格檢測方法對網(wǎng)格質量較為敏感,且能較好地體現(xiàn)數(shù)值方法對網(wǎng)格質量的依賴性。第七章是本文的結束語。概括了論文的主要創(chuàng)新點,并對未來的工作方向進行了展望。
[Abstract]:As one of the most ancient numerical discretization methods , the finite difference method is very popular with beginners because of its simple and easy - to - understand . However , it is difficult to find trace in the application of finite difference method based on finite difference method . In this paper , the influence of grid smoothness on the precision of geometric conservation laws is analyzed in detail . The results of this paper are as follows :

【學位授予單位】：中國空氣動力研究與發(fā)展中心
【學位級別】：博士
【學位授予年份】：2015
【分類號】：O241.82

【相似文獻】

相關期刊論文 前10條

1 王加玲;孫志忠;;熱方程非線性邊界條件問題的有限差分方法(英文)[J];高等學校計算數(shù)學學報;2008年04期

2 舒阿秀;;時滯非線性拋物型方程時間周期解的有限差分方法[J];合肥工業(yè)大學學報(自然科學版);2008年02期

3 舒阿秀;;一類時滯非線性拋物型方程時間周期解的有限差分方法[J];河南科技大學學報(自然科學版);2009年01期

4 馬宗源;廖紅建;;雙剪統(tǒng)一彈塑性有限差分方法研究[J];計算力學學報;2012年01期

5 楊志峰,王?;高精度有限差分方法研究進展[J];自然科學進展;1999年09期

6 郭南南;馬犀大;;內置偏心電極的空心圓管中有限差分方法的應用[J];科技信息;2006年S5期

7 王海明;;一類變系數(shù)半線性拋物型方程的有限差分方法[J];北華大學學報(自然科學版);2007年03期

8 李明;;有限差分方法求任意形狀勢壘的透射系數(shù)[J];科技信息;2013年21期

9 張寶琳;并行求解初邊值問題的有限差分方法研究[J];應用數(shù)學與計算數(shù)學學報;1992年01期

10 高智;攝動有限差分方法研究進展[J];力學進展;2000年02期

相關會議論文 前9條

1 張德良;胡宗民;;攝動有限差分方法及其應用[A];第十屆全國激波與激波管學術討論會論文集[C];2002年

2 唐少強;;晶體多尺度計算的有限差分方法[A];慶祝中國力學學會成立50周年暨中國力學學會學術大會’2007論文摘要集（下）[C];2007年

3 朱力立;張德良;高智;;不可壓流動攝動有限差分方法及其應用[A];計算流體力學研究進展——第十一屆全國計算流體力學會議論文集[C];2002年

4 周逸杰;陳效群;鄧建松;;B忮zier曲線求值方法的比較與選擇[A];中國幾何設計與計算新進展2007——第三屆中國幾何設計與計算大會論文集[C];2007年

5 張偉;陳曉非;;有限差分方法中存在復雜形狀地表時的自由表面條件[A];中國地球物理.2003——中國地球物理學會第十九屆年會論文集[C];2003年

6 唐少強;;晶體多尺度計算的有限差分方法[A];2006年全國固體力學青年學者研討會論文摘要文集[C];2006年

7 唐少強;;晶體多尺度計算的有限差分方法[A];全國計算物理學會第六屆年會和學術交流會論文摘要集[C];2007年

8 陳振興;陳龍偉;;基于有限差分方法的金屬桿與混凝土靶板侵徹的數(shù)值模擬[A];第五屆全國計算爆炸力學會議論文摘要[C];2012年

9 尹兆華;李巧紅;常磊;;不變形雙滴的熱毛細遷移及相互作用[A];中國力學大會——2013論文摘要集[C];2013年

相關博士學位論文 前6條

1 秦雯娣;幾類偏微分方程的動力學相容的非標準有限差分方法[D];哈爾濱工業(yè)大學;2015年

2 閔耀兵;高階精度有限差分方法幾何守恒律研究[D];中國空氣動力研究與發(fā)展中心;2015年

3 鄭罡;頻域有限差分方法的改進研究[D];電子科技大學;2010年

4 廖翠萃;多辛變分積分子及非標準有限差分方法[D];哈爾濱工業(yè)大學;2013年

5 張守慧;拋物型方程的幾種可并行的有限差分方法[D];山東大學;2009年

6 羅志強;非線性水波數(shù)值模擬[D];南開大學;2011年

相關碩士學位論文 前10條

1 蔡成健;有限差分方法與隨機偏微分方程的數(shù)值模擬[D];廣州大學;2009年

2 魏瑋;非線性拋物型方程組時間周期解的有限差分方法[D];華東師范大學;2005年

3 舒阿秀;時滯反應擴散方程組的時間周期解的有限差分方法[D];華東師范大學;2009年

4 吳魏;BBM-Burgers方程的有限差分方法[D];南京航空航天大學;2007年

5 鄒年容;一維無界域上薛定諤方程的有限差分方法[D];湘潭大學;2007年

6 張洪波;非線性反應擴散方程組的高階緊有限差分方法[D];華東師范大學;2008年

7 陳卓;低擴散分數(shù)階微分方程的有限差分方法[D];南京航空航天大學;2011年

8 王迪;一類發(fā)展方程的有限差分方法[D];信陽師范學院;2015年

9 姜海云;高階Lidstone邊值問題的四階緊有限差分方法[D];華東師范大學;2007年

10 柏琰;對稱正則長波（SRLW）方程的有限差分方法[D];南京航空航天大學;2005年

，

本文編號：1537487