【摘要】：隨著計(jì)算機(jī)技術(shù)的發(fā)展,計(jì)算流體力學(xué)(CFD)在工程領(lǐng)域得到了廣泛應(yīng)用,并發(fā)揮著重要的作用,但是CFD仍然存在一些具有挑戰(zhàn)性的問(wèn)題。當(dāng)運(yùn)用CFD解決現(xiàn)代工業(yè)界所面臨的各種設(shè)計(jì)、計(jì)算以及復(fù)雜幾何形狀下的流動(dòng)等問(wèn)題時(shí),數(shù)值計(jì)算精度和湍流模型仍然是CFD工作者面臨的一個(gè)難題。為此,本文進(jìn)行以下研究:首先基于區(qū)域分裂技術(shù),采用求解偏微分方程并結(jié)合基本的代數(shù)方法,生成繞ONERA M6(M6)機(jī)翼的高質(zhì)量多塊結(jié)構(gòu)化拼接網(wǎng)格�？臻g采用Roe的通量差分分裂(Flux Difference Splitting)迎風(fēng)格式,時(shí)間采用隱式推進(jìn)方法,通過(guò)數(shù)值模擬M6機(jī)翼跨音速流場(chǎng)以探究?jī)煞N限制器:minimum-modulus(min-mod)限制器和modified 3rd-Upwind偏置限制器在跨聲速繞流中的性能�？刂品匠虨槔字Z平均Navier-Stokes(RANS)方程耦合Spalart-Allmaras(SA)一方程湍流模型,為了加速收斂,計(jì)算還采用了多重網(wǎng)格方法。結(jié)果表明:min-mod限制器預(yù)測(cè)到的激波質(zhì)量相對(duì)較好,modified 3rd-Upwind偏置限制器的收斂性和精度好于min-mod限制器。其次針對(duì)ONERA M6機(jī)翼,通過(guò)求解RANS方程分別耦合一方程Spalart-Allmaras(SA)湍流模型、二方程Wilcox’s k-?(Wilcox)和Menter’s SST(SST)湍流模型,考察三種湍流模型在跨聲速繞流中的性能和對(duì)氣動(dòng)特征的影響。為了提高精度、避免激波附近的數(shù)值振蕩,計(jì)算中使用了modified 3rd-Upwind偏置限制器。通過(guò)對(duì)比分析不同湍流模型之間的差異,獲得了對(duì)湍流模型的選取和改進(jìn)具有一定參考價(jià)值的結(jié)果。最后基于“超立方體”概念,提出了一種通用的高質(zhì)量多塊拼接結(jié)構(gòu)化網(wǎng)格構(gòu)建方法,生成了繞細(xì)長(zhǎng)旋成體的計(jì)算網(wǎng)格;通過(guò)求解RANS方程分別耦合SA、Wilcox和SST三種湍流模型,數(shù)值模擬了細(xì)長(zhǎng)旋成體在超聲速大攻角下的分離流流場(chǎng),比較了SA、Wilcox和SST三種湍流模型的差異,獲得了與實(shí)驗(yàn)結(jié)果吻合較好的物面壓力分布。三種湍流模型較好地模擬了主渦強(qiáng)度、渦核位置以及靠近壁面的二次分離渦。本文通過(guò)計(jì)算分別對(duì)比分析了不同限制器、不同湍流模型之間的差異,獲得了比較理想的結(jié)果,驗(yàn)證了本文算法和求解技術(shù)的可行性。本文結(jié)果對(duì)提高CFD的模擬精度、湍流模型的改進(jìn)及限制器的構(gòu)建具有一定的參考價(jià)值。
【圖文】：

n 1nQ Q Q .5 多重網(wǎng)格法多重網(wǎng)格法從工程應(yīng)用角度講是一種多分辨率的算法，多重網(wǎng)格的生成方法一重疊式和聚合式，重疊式對(duì)于幾何結(jié)構(gòu)的某些關(guān)鍵特征的粗網(wǎng)格不能充分表達(dá)，能使粗網(wǎng)格和細(xì)網(wǎng)格間的傳遞關(guān)系變得簡(jiǎn)單，在很大程度上減少工作量，避免某幾何外形不能產(chǎn)生粗網(wǎng)格的情況。在采用多重網(wǎng)格技術(shù)時(shí)，需要運(yùn)用循環(huán)準(zhǔn)則確定某時(shí)刻從一層網(wǎng)格轉(zhuǎn)換到另一，循環(huán)準(zhǔn)則主要分為自適應(yīng)循環(huán)方法和確定的循環(huán)方法。由于用自適應(yīng)循環(huán)方法前網(wǎng)格層高頻誤差的消除較困難，故本文采用確定的循環(huán)方法，，該方法簡(jiǎn)單、穩(wěn)含兩種最常用的循環(huán)：V 循環(huán)和 W 循環(huán)[44]。圖 2.1 為 V 型循環(huán)方式示意圖， W 型循環(huán)方式示意圖。

n 1nQ Q Q .5 多重網(wǎng)格法多重網(wǎng)格法從工程應(yīng)用角度講是一種多分辨率的算法，多重網(wǎng)格的生成方法一重疊式和聚合式，重疊式對(duì)于幾何結(jié)構(gòu)的某些關(guān)鍵特征的粗網(wǎng)格不能充分表達(dá)，能使粗網(wǎng)格和細(xì)網(wǎng)格間的傳遞關(guān)系變得簡(jiǎn)單，在很大程度上減少工作量，避免某幾何外形不能產(chǎn)生粗網(wǎng)格的情況。在采用多重網(wǎng)格技術(shù)時(shí)，需要運(yùn)用循環(huán)準(zhǔn)則確定某時(shí)刻從一層網(wǎng)格轉(zhuǎn)換到另一，循環(huán)準(zhǔn)則主要分為自適應(yīng)循環(huán)方法和確定的循環(huán)方法。由于用自適應(yīng)循環(huán)方法前網(wǎng)格層高頻誤差的消除較困難，故本文采用確定的循環(huán)方法，該方法簡(jiǎn)單、穩(wěn)含兩種最常用的循環(huán)：V 循環(huán)和 W 循環(huán)[44]。圖 2.1 為 V 型循環(huán)方式示意圖， W 型循環(huán)方式示意圖。
【學(xué)位授予單位】：長(zhǎng)安大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2018
【分類號(hào)】：O241.82

【參考文獻(xiàn)】

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

1 涂國(guó)華;燕振國(guó);趙曉慧;馬燕凱;毛枚良;;SA和SST湍流模型對(duì)高超聲速邊界層強(qiáng)制轉(zhuǎn)捩的適應(yīng)性[J];航空學(xué)報(bào);2015年05期

2 唐志共;張益榮;陳堅(jiān)強(qiáng);毛枚良;張毅鋒;劉化勇;;更準(zhǔn)確、更精確、更高效——高超聲速流動(dòng)數(shù)值模擬研究進(jìn)展[J];航空學(xué)報(bào);2015年01期

3 張翔;閻超;;高超聲速熱流計(jì)算湍流模型性能評(píng)估[J];北京航空航天大學(xué)學(xué)報(bào);2015年02期

4 周培培;;基于結(jié)構(gòu)網(wǎng)格的柵格翼繞流數(shù)值模擬[J];空氣動(dòng)力學(xué)學(xué)報(bào);2014年03期

5 聶智軍;成水燕;麻蓉;顏洪;;基于RANS方法的民用飛機(jī)網(wǎng)格分析[J];航空計(jì)算技術(shù);2013年05期

6 韓濤;;湍流模型在民機(jī)跨聲速繞流計(jì)算的應(yīng)用[J];航空計(jì)算技術(shù);2013年03期

7 尹翔;馮金富;張俊祥;;一種典型翼型在湍流模型中的數(shù)值仿真分析[J];計(jì)算機(jī)仿真;2013年04期

8 劉景源;;SST湍流模型在高超聲速繞流中的改進(jìn)[J];航空學(xué)報(bào);2012年12期

9 鄭秋亞;王寶園;劉三陽(yáng);;繞翼身組合體高質(zhì)量網(wǎng)格設(shè)計(jì)和阻力計(jì)算[J];應(yīng)用力學(xué)學(xué)報(bào);2009年03期

10 劉仙名,符松;幾種k-ε模式在細(xì)長(zhǎng)體分離流計(jì)算中的應(yīng)用[J];計(jì)算力學(xué)學(xué)報(bào);2004年04期

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

1 鄭秋亞;基于Navier-Stokes方程的復(fù)雜流動(dòng)數(shù)值模擬精度與并行計(jì)算研究[D];西安電子科技大學(xué);2011年

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

1 馬戎;基于動(dòng)態(tài)混合網(wǎng)格的氣動(dòng)/運(yùn)動(dòng)耦合一體化計(jì)算方法研究[D];中國(guó)空氣動(dòng)力研究與發(fā)展中心;2015年

2 陸應(yīng)龍;基于N-S方程的飛行器氣動(dòng)力氣動(dòng)熱數(shù)值研究[D];上海交通大學(xué);2014年

3 李杰;Navier-Stokes方程數(shù)值模擬及湍流模型研究[D];西北工業(yè)大學(xué);2006年

4 王玉玲;歐拉方程的高效數(shù)值方法及三角翼大迎角脫體渦特性的研究[D];西北工業(yè)大學(xué);2005年

本文編號(hào)：2631086