發(fā)布時間：2018-12-07 16:33

【摘要】：發(fā)展了一套基于伴隨方法的Euler方程自適應(yīng)網(wǎng)格DG求解方法,采用當?shù)鼐W(wǎng)格上的輸出變量的誤差估計作為離散誤差指示器驅(qū)動網(wǎng)格自適應(yīng).在離散誤差指示器的構(gòu)造中,采用GMRES方法求解伴隨方程獲得伴隨變量,將p階流場變量和伴隨變量映射到p+1階函數(shù)空間,并進行有限次塊對角Jacobi迭代得到近似p+1階細網(wǎng)格解.壁面邊界通量及輸出變量計算僅依賴邊界外側(cè)值以保證伴隨相容性.首先采用NACA0012翼型不同精度的阻力系數(shù)結(jié)果驗證了誤差估計的可靠性,然后對NACA0012翼型亞聲速、跨聲速流場和圓柱6馬赫高超聲速繞流進行了自適應(yīng)加密計算研究.研究結(jié)果表明,基于伴隨方法的網(wǎng)格自適應(yīng)可有效提高阻力系數(shù)等輸出變量的模擬準確性,在亞聲速NACA0012翼型計算中使用約17%的全局加密自由度獲得了與全局加密精度相當?shù)淖枇ο禂?shù).
[Abstract]:A set of adaptive mesh DG solutions for Euler equation based on adjoint method is developed. The error estimation of the output variables on the local grid is used as the discrete error indicator to drive the mesh adaptation. In the construction of discrete error indicator, the adjoint variable is obtained by solving adjoint equation by GMRES method, and the p order flow field variable and adjoint variable are mapped to the function space of p 1 order. The finite block diagonal Jacobi iteration is used to obtain the approximate p1 order fine grid solution. The calculation of wall boundary flux and output variables only depends on the lateral boundary value to ensure adjoint compatibility. The reliability of the error estimation is verified by the results of the drag coefficients of the NACA0012 airfoil with different accuracy. Then, the adaptive encryption calculation of the subsonic velocity, transonic flow field and the hypersonic flow around the cylinder 6 Mach is carried out for the NACA0012 airfoil. The results show that the mesh adaptation based on adjoint method can effectively improve the simulation accuracy of output variables such as drag coefficient. In the calculation of subsonic NACA0012 airfoil, about 17% of the degree of freedom of global encryption is used to obtain the resistance coefficient equivalent to the precision of global encryption.
【作者單位】： 中國航天空氣動力技術(shù)研究院;
【分類號】：O241.8

【相似文獻】

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

1 丁吉豫;關(guān)于人口預(yù)報的一種自適應(yīng)方法(摘要)[J];齊齊哈爾師范學(xué)院學(xué)報(自然科學(xué)版);1982年02期

2 黃季墀,夏人偉,何慶芝;一種具有自適應(yīng)能力的結(jié)構(gòu)綜合方法[J];航空學(xué)報;1988年03期

3 鄒青松;;求解一般橢圓形方程不標記震蕩項的自適應(yīng)方法[J];高等學(xué)校計算數(shù)學(xué)學(xué)報;2009年04期

4 周俊;;曲線最小二乘擬合的自適應(yīng)方法[J];太原師范學(xué)院學(xué)報(自然科學(xué)版);2011年03期

5 蔡顯新;王文凱;蔣燕英;江萍;;一種有效的網(wǎng)格自適應(yīng)方法[J];計算力學(xué)學(xué)報;2007年02期

6 于春肖;彈性問題的自適應(yīng)數(shù)值解析研究[J];燕山大學(xué)學(xué)報;2005年01期

7 竇海勇;研究稀有事件的自適應(yīng)弦方法[J];華南理工大學(xué)學(xué)報(自然科學(xué)版);2005年04期

8 謝承迪;自適應(yīng)運籌濾波方法[J];數(shù)值計算與計算機應(yīng)用;1994年01期

9 周宏宇;王愛民;;非線性偏微分方程數(shù)值求解的自適應(yīng)方法研究[J];計算機工程與應(yīng)用;2011年20期

10 梅樹立,陸啟韶,張森文;求解非線性偏微分方程的自適應(yīng)小波精細積分法[J];計算物理;2004年06期

相關(guān)會議論文 前4條

1 袁駟;和雪峰;;一個高效的一維有限元自適應(yīng)求解的新方案——第十三屆全國結(jié)構(gòu)工程學(xué)術(shù)大會特邀報告[A];第十三屆全國結(jié)構(gòu)工程學(xué)術(shù)會議論文集（第Ⅰ冊）[C];2004年

2 汪新;;自適應(yīng)邊界元方法[A];計算力學(xué)研究與進展——中國力學(xué)學(xué)會青年工作委員會第三屆學(xué)術(shù)年會論文集[C];1999年

3 袁駟;和雪峰;;一個高效的一維有限元自適應(yīng)求解的新方案[A];工程力學(xué)學(xué)術(shù)研討會論文集[C];2004年

4 徐旭;郭崇慧;;稀疏成分分析的簡要綜述[A];中國運籌學(xué)會第十屆學(xué)術(shù)交流會論文集[C];2010年

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

1 杜炎;基于EEP法的一維非線性有限元自適應(yīng)分析[D];清華大學(xué);2012年

2 肖嘉;基于EEP法的線法二階常微分方程組有限元自適應(yīng)分析[D];清華大學(xué);2009年

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

1 陳星;帶齊次混合邊界特征值問題的一種基于多尺度離散的有限元自適應(yīng)算法[D];貴州師范大學(xué);2015年

2 余媛媛;基于移位反迭代的非協(xié)調(diào)Crouzeix-Raviart有限元自適應(yīng)方法求Laplace特征值問題[D];貴州師范大學(xué);2015年

3 王彪;弱不連續(xù)問題的p型自適應(yīng)有限元及其快速求解方法[D];湘潭大學(xué);2015年

4 白思林;h-，p-，，hp-自適應(yīng)邊界元方法研究[D];燕山大學(xué);2009年

5 舒冬;二維自適應(yīng)有限元靜力分析方法研究[D];中南大學(xué);2012年

6 楊銀;奇異攝動問題的自適應(yīng)方法[D];湘潭大學(xué);2006年

7 鄭國賢;Level Set方程的求解與應(yīng)用[D];浙江大學(xué);2006年

8 竇海勇;研究復(fù)雜系統(tǒng)中稀有事件的轉(zhuǎn)化路徑的數(shù)值方法[D];清華大學(xué);2004年

9 楊繼明;奇異攝動問題自適應(yīng)移動網(wǎng)格迭代算法研究[D];湘潭大學(xué);2004年

10 鄭雙超;多尺度波動問題的耦合求解方法[D];清華大學(xué);2011年

本文編號：2367467