粘性流動模擬中基于伴隨方程的網(wǎng)格自適應(yīng)
本文關(guān)鍵詞:粘性流動模擬中基于伴隨方程的網(wǎng)格自適應(yīng),,由筆耕文化傳播整理發(fā)布。
【摘要】:本文發(fā)展了一種以提高目標函數(shù)計算值的準確度為目的,基于伴隨方程的網(wǎng)格自適應(yīng)方法。伴隨方程的解表征了目標函數(shù)對于粘性流動方程殘值的靈敏度。利用伴隨解,目標函數(shù)的誤差可以表示為可計算修正項與剩余誤差項兩部分之和,可計算修正項即為伴隨方程解與粘性流動方程殘值的內(nèi)積;谑S嗾`差項構(gòu)建了對于粘性流動模擬具有良好穩(wěn)定性的網(wǎng)格自適應(yīng)探測器。通過減小剩余誤差項,提高可計算修正項的精度,從而間接提升目標函數(shù)值的準確度。采用插值法獲得自適應(yīng)參數(shù)計算中所需要的細網(wǎng)格變量值。采用有限體積法median-dual格點格式求解二維可壓層流方程,對流項的離散采用中心格式,在遠場邊界應(yīng)用無反射條件,在物面應(yīng)用無滑移和恒溫壁邊界條件,時間推進采用顯式五步龍格庫塔格式。伴隨方程的求解采用離散方法,并添加時間導(dǎo)數(shù)項,時間推進同樣采用五步龍格庫塔格式。推導(dǎo)了對于中心格式的對流伴隨通量計算方法和人工耗散項的處理方法。采用薄層近似等方法簡化粘性伴隨通量的計算。強邊界條件使得伴隨解在域內(nèi)和物面具有不同的性質(zhì),采用了先獲得域內(nèi)伴隨解,再對物面網(wǎng)格點進行后處理得到物面伴隨解的方法。討論了計算伴隨方程右端目標函數(shù)線性導(dǎo)數(shù)項的方法。最后選擇具有重要工程應(yīng)用意義的升力系數(shù)、阻力系數(shù)等作為目標函數(shù),將網(wǎng)格自適應(yīng)方法和誤差修正技術(shù)應(yīng)用于圓柱繞流和NACA0012翼型繞流模擬。數(shù)值實驗的結(jié)果表明本文方法能有效地提高目標函數(shù)計算值的準確性。
【關(guān)鍵詞】:網(wǎng)格自適應(yīng) 伴隨方程 粘性流動 目標函數(shù) 誤差估算
【學(xué)位授予單位】:南京航空航天大學(xué)
【學(xué)位級別】:碩士
【學(xué)位授予年份】:2015
【分類號】:O357
【目錄】:
- 摘要4-5
- ABSTRACT5-9
- 注釋表9-10
- 第一章 緒論10-13
- 1.1 研究目的與意義10
- 1.2 國內(nèi)外研究現(xiàn)狀10-11
- 1.3 本文的工作11-13
- 第二章 伴隨自適應(yīng)原理13-20
- 2.1 目標函數(shù)的誤差估算及修正13-14
- 2.2 自適應(yīng)參數(shù)與準則14-16
- 2.3 插值算子16-18
- 2.4 網(wǎng)格自適應(yīng)過程18-19
- 2.5 本章小結(jié)19-20
- 第三章 流動方程的求解20-29
- 3.1 控制方程20-21
- 3.2 數(shù)值離散21-26
- 3.2.1 對流項離散22-23
- 3.2.2 粘性項離散23-24
- 3.2.3 邊界條件24-25
- 3.2.4 時間推進25-26
- 3.3 算例分析26-28
- 3.4 本章小結(jié)28-29
- 第四章 伴隨方程的求解29-41
- 4.1 伴隨方法簡介29
- 4.2 離散伴隨方程29-31
- 4.3 伴隨通量的計算31-35
- 4.3.1 對流伴隨通量的計算32-33
- 4.3.2 粘性伴隨通量的計算33-35
- 4.3.3 人工耗散項的處理35
- 4.4 物面伴隨解的計算35-38
- 4.5 右端項的計算38-39
- 4.6 本章小結(jié)39-41
- 第五章 數(shù)值試驗與分析41-51
- 5.1 圓柱繞流42-43
- 5.2 NACA0012 翼型繞流43-49
- 5.3 雙NACA0012 翼型繞流49-50
- 5.4 本章小結(jié)50-51
- 第六章 總結(jié)與展望51-53
- 6.1 全文工作總結(jié)51
- 6.2 未來的工作展望51-53
- 參考文獻53-56
- 致謝56-57
- 在學(xué)期間的研究成果及學(xué)術(shù)論文57
【相似文獻】
中國期刊全文數(shù)據(jù)庫 前10條
1 李亭鶴,閻超,李躍軍;重疊網(wǎng)格技術(shù)中割補法的研究與改進[J];北京航空航天大學(xué)學(xué)報;2005年04期
2 袁益謙;網(wǎng)格高速攝影[J];高速攝影與光子學(xué);1981年02期
3 袁益謙;網(wǎng)格高速攝影成象特性的研究[J];光學(xué)學(xué)報;1983年07期
4 卓向正;;一個利用微機求網(wǎng)格點上氣象要素值的方案[J];浙江氣象科技;1986年03期
5 周春華;;一種基于預(yù)估-修正策略的非定常流網(wǎng)格自適應(yīng)加密/稀疏方法[J];計算物理;2013年05期
6 葉強強;余燁;張靜;劉曉平;;基于網(wǎng)格的球面虛擬實景空間搭建方法[J];華南師范大學(xué)學(xué)報(自然科學(xué)版);2012年04期
7 劉君;白曉征;王巍;;網(wǎng)格之間流場信息傳遞的高精度方法和驗證[J];空氣動力學(xué)學(xué)報;2010年01期
8 曹建發(fā),陳晏清;多段翼型貼體網(wǎng)格的生成方法[J];北京航空航天大學(xué)學(xué)報;1993年01期
9 張政,何世平,續(xù)伯欽;自動網(wǎng)格形狀測試技術(shù)[J];中國科學(xué)技術(shù)大學(xué)學(xué)報;1997年02期
10 何躍龍;李盾;;非結(jié)構(gòu)粘性直角網(wǎng)格物面相貫線修正[J];計算機與數(shù)字工程;2009年12期
中國重要會議論文全文數(shù)據(jù)庫 前9條
1 王瑞利;付漢清;;流體力學(xué)數(shù)值模擬中的網(wǎng)格自適應(yīng)技術(shù)[A];中國工程物理研究院科技年報(1999)[C];1999年
2 李盾;陳靖;方華;;復(fù)雜外形三維非結(jié)構(gòu)粘性直角網(wǎng)格的生成[A];中國造船工程學(xué)會電子技術(shù)學(xué)術(shù)委員會2006學(xué)術(shù)年會論文集(上冊)[C];2006年
3 葉強強;余燁;張靜;劉曉平;;基于網(wǎng)格的球面虛擬實景空間搭建方法[A];安徽首屆科普產(chǎn)業(yè)博士科技論壇——暨社區(qū)科技傳播體系與平臺建構(gòu)學(xué)術(shù)交流會論文集[C];2012年
4 鄭冠男;楊國偉;;三維非結(jié)構(gòu)動網(wǎng)格計算方法研究[A];第十三屆全國激波與激波管學(xué)術(shù)會議論文集[C];2008年
5 周磊;譚偉偉;;聚合多重網(wǎng)格的體網(wǎng)格可視化技術(shù)[A];探索 創(chuàng)新 交流(第4集)——第四屆中國航空學(xué)會青年科技論壇文集[C];2010年
6 晉長秋;曹雄;;二維拉氏程序中的重分[A];中國工程物理研究院科技年報(2000)[C];2000年
7 陳立崗;梁義;陳建軍;鄭耀;;保特征的曲面網(wǎng)格重生成方法[A];慶祝中國力學(xué)學(xué)會成立50周年暨中國力學(xué)學(xué)會學(xué)術(shù)大會’2007論文摘要集(下)[C];2007年
8 楊永;;基于網(wǎng)格計算的CFD模擬可信度分析[A];計算流體力學(xué)研究進展——第十二屆全國計算流體力學(xué)會議論文集[C];2004年
9 龔也君;;浸入界面法快速生成網(wǎng)格[A];第十七屆全國反應(yīng)堆結(jié)構(gòu)力學(xué)會議論文集[C];2012年
中國博士學(xué)位論文全文數(shù)據(jù)庫 前10條
1 郭e
本文編號:301020
本文鏈接:http://www.sikaile.net/shoufeilunwen/benkebiyelunwen/301020.html