《基于简单weno-间断galerkin的euler方程自适应计算-吴泽艳.pdf》由会员分享，可在线阅读，更多相关《基于简单weno-间断galerkin的euler方程自适应计算-吴泽艳.pdf（9页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、2016年 4月第42卷第4期北京航空航天大学学报Journal of Beijing University of Aeronautics and AstronauticsApril 2016V0142 No4http：7f bhxbbuaaeduCYI jbuaabuaaeducnDOI：1013700jbh1001596520150237基于简单WENO-间断Galerkin的Euler方程自适应计算吴泽艳1，(1三峡大学水利与环境学院，宜昌443002王立峰2一，武哲22北京航空航天大学航空科学与工程学院，北京100083)摘 要：为了得到Euler方程的高精度、高分辨率数值解，介绍了间
2、断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器一间断Galerkin方法应用到凸边四边形单元上，通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号，实现了基于“问题单元”的局部网格加密自适应计算。对若干典型问题进行编程计算，结果表明，简单WENO限制器可以应用到曲边四边形单元上，且可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。关 键 词：间断Galerkin方法；简单WENO限制器；Euler方程；自适应计算；曲边四边形单元中图分类号：0242文献标识码：A 文章编号：100159
3、65(2016)04-0806-09间断Galerkin(DG)方法是由Reed和Hill1于1973年为了解决线性输运问题而提出，此后该方法得到了不断发展。20世纪90年代，Cockburn和Shu旧1提出的RungeKutta间断Galerkin(RKDG)方法尤其引人注目，在气动力学、水动力学、电磁场以及气动声学等领域的数值模拟中得到了广泛应用。1998年前后，Bassi和Rebay1、Cockburn和Shu1又将DG方法推广到了黏性流动的计算。近期，Shu。对DG方法作了较详细的综述。DG方法是有限元方法与有限体积方法的结合，在基函数和检验函数的运用方面与有限元方法相似，而在对方程性
4、质的近似方面与有限体积方法更接近，因此有若干优点。首先，与有限元方法一样，间断有限元方法可以方便地提高空间和时间方向的计算精度。在空间上，DG方法可以方便地实施自适应计算，其有2个方面的含义：对网格单元的自适应加密。DG方法允许单元间出现“悬挂节点”，因此网格加密非常简单。DG方法允许不同的单元采用不同阶的多项式。这2个方面分别称作h自适应和P自适应。其次，由DG方法形成的常微分方程组，其质量矩阵是分块对角的，甚至可以是单位矩阵，只要在各个单元内选择合适的正交多项式基函数。因为分块对角矩阵求逆相当于各个块矩阵的求逆，所以质量矩阵的求逆计算量极小。分块对角矩阵的逆仍然是分块对角矩阵，所以存储量也
5、极小。最后，间断有限元方法特别适于并行计算。从已有文献可知，间断有限元方法的并行计算效率在80以上旧1。与用谱方法求解间断问题会出现Gibbs现象类似，如果解在单元内出现间断，则用高阶多项式作为形函数会导致近似解的振荡。为了获得稳定收稿日期：2015-0418；录用日期：2015-07-17；网络出版时间：201510-14 15：18网络出版地址：WWWcnkinetkcmsdetail112625V，201510141518003html基金项目：三峡大学人才科研启动基金(KJ20148031)通讯作者：Tel：01082339560 Email：wanglifengase，buaaedu
6、cn引用格式：吴泽艳，王立峰，武哲基于简单WENO-间断Galerkin的Euler方程自适应计算J北京航空航天大学学撮，2016，42f4)：806-814WU Z Y，WANC L F，WU ZAdaptive simple WENO limiter-discontinuous Galerkin method for Euler equations【l 3Journal of Beijing University of Aeronautics and Astronautics，2016，42(4)：806814(in Chinese)万方数据第4期 吴泽艳，等：基于简单WENO一间断Gal
7、erkin的Euler方程自适应计算807的近似解，需要引入限制器对解作修正。早期，研究人员开发了一些经典的限制器，如基于minmod函数的TVB限制器。、Biswas等。81的moment限制器、Burbeau等1的改进moment限制器、保单调限制器0。和改进保单调限制器-等。之后，类似从差分法中借鉴数值通量概念，间断有限元方法从差分法中借鉴了WENO的思想1，通过对本单元及相邻若干单元的变量的单元平均作WENO重构，形成了一种新的限制器，得到了高分辨率的结果，成为近年的研究热点。在构造一个合适的WENO限制器前，需要确定哪些单元需要被限制，即识别哪些单元是“问题单元”。早期研究人员按照文
8、献7提出的限制器来识别问题单元，即如果某单元需要用到该限制器，则该单元是“问题单元”，否则不是。文献13的研究表明，DG方法在单元的出流边界具有很强的超收敛特性，基于此特性构造了一种间断侦测器。在国内，张来平等4。联合应用该间断侦测器与经典限制器进行了激波捕捉，取得了较好的效果。最近，Zhong和Shu1 5。提出了一种简单WENO限制器，取得了较好的间断捕捉效果，而且因为这种限制器只需要用到“问题单元”及其相邻单元的信息，所以编程简单。本文的计算采用这种简单WENO限制器。1 数学模型考虑Euler方程的守恒形式：iOu 4-兰粤4-掣堕：o (1)一一一=I，a az ay式中：pp“pv
9、EF(U)=G(U)=P珏pu。+PPM(E+P)MpvPM”pv+P(E+P)以及P=(71)p8 (4)式中：s为单位质量内能；7为比热比，对于空气，取y=14。2空间及时间离散21 DG空间离散的一般理论用DG方法对式(1)进行空间离散。设，。为对区域n的一个有限剖分，单元K E，。，I，(K)为单元K上的局部有限元空间，为P(k0)次多项式的集合。在时间0，T内，在如下间断有限元空间中寻找近似解玑(X，t)，这里X=(戈，Y)。Vh=口L。(n)：口I KV(K)，V KF运用加权余量法，在单元K上用权函数W乘以式(1)，并在单元K上积分，有瓤u(X,t)埘(x)do=一J V厂(u(
10、X，t)W(X)d仃 (5)式中：，(U(X，t)=(F(U(X，t)，G(U(X，t)用近似解巩(X，t)K代替式(5)中的精确解u(X，t)，用w。Vh代替权函数W，利用Green公式有瓤Uh(砷)Wh(X)dO=一f，(U。(x，f)n e,K埘。(x)d厂+I f(U。(x，)V埘。(X)d门 (6)式中：e为单元K的边界；Ile,K为单元K的边界外法线单位向量。对式(6)用数值积分代替解析积分有J，(u。(X，)咒。，。W(x)d厂一601，(U。(x)。n e,KWn(x。，)(2) f，(Uh(x，f)Vwh(x)d力。JK(-Om-厂(Uh(x)。Vw。(x。)I K其中：队u
11、、”和P分别为气体的密度、戈方向速度分量、Y方向速度分量和压力；E为单位体积的总能量。有E=pc+p(M2+u2)2 (3)式中：I e l为单元边长；l K I为单元面积；为高斯积分系数。间断有限元方法从差分法中吸收了数值通量的概念，将边界的通量，(巩(x，t)n e用数值,K通量函数h“(X，t)代替，有瓤州X,t)Wh(X)do=万方数据808 北京航空航天大学学报 2016年一w zh“(x，)Wh(x。z)+P E JK l=l盯(U。(x。，#)V埘一(x。)l K (7)m=1数值通量函数h5(X，t)在单元边界OK上的值依赖于近似解在边界两边的值，即由单元K的内部的值乩(x：，
12、t)与单元K的外部的值U。(x：2，t)决定。矿“(X，t)=h“(U。(x：，t)，U。(xo；，t) (8)数值通量要求满足相容性、单调性、全局Lipsehitz连续性以及守恒性。满足这些条件的数值通量常见类型有Godunov型、EngquistOsher型、LaxFriedrichs型、Local LaxFriedrichs型以及Roe型等。一般来说，Roe型与Local LaxFriedrichs型计算效果接近，随着逼近多项式阶次的提高，它们对近似解的影响差别很小。对于线性情形，这些数值通量与迎风型数值通量一致；对于非线性问题，因为Godunov型的数值通量产生的人工黏性最小而成为最佳
13、选择。特别地，LocalLaxFriedrichs型数值通量比Godunov型数值通量产生更多的人工黏性。实际计算中，大多采用简化的LaxFriedrichs型数值通量，而在处理高维双曲系统时，选择Local LaxFriedrichs型数值通量更为方便。本文计算中均采用Local LaxFriedrichs型数值通量。在单元中选取基函数为咖，咖：，咖，这里基函数是关于坐标咒和Y的多项式，近似解可表示为u。(x，)=q()咖j(x)代人式(6)，且令权函数W。=咖，则式(6)可以写为常微分方程组的形式：JM芋uh(t)=L(U(t) (9)U式中：M为质量矩阵，是由各个单元的质量矩阵L。=R。
14、=Mx 2；上咖i(x)咖，(x)d门集成的分块对角矩阵。22时间离散式(9)可写为面d un(t)=M一1L(Uh()为描述方便，可写作面dU=，(U) (10)使用如下TVD型显式RungeKutta方法求解式(10)：f U=U”+A，(U”)=驴)+了1渺心，(1)I旷川=了1旷+j2。U(2)+护，(扩)为满足稳定性要求，时间步长取为At力一(2k+1)ftax(u；栉。，。f十c，珏。r露。，。+C r)do-JaK(12)式中：力。为单元面积；k为形函数的次数；上标z和r分别为本单元和相邻单元；C为当地声速。3简单WENO限制器简单WENO限制器由Zhong和Shu 015提出。
15、简单WENO限制器重构模板小，仅需要用到“问题单元”及其相邻单元的信息，所以特别适合于h自适应计算。这里以三角形单元为例对简单WENO限制器作详细介绍。设已用某种方法判断出哪些单元为“问题单元”，记“问题单元”的第i条边的外法向为n=(rb。n。，)1，Jacobian矩阵为(，(U)，g 7(H)n。，则该Jacobian矩阵的左右特征向量矩阵为B2+(urt。，+秽n。y)c Bl M+nixc Bl割+n。yc B1 一12 2 2凡。yMnix勘 一凡iy n珏01一B2 Bl M Bl 一B1B2一(un。，+vn。y)c BIunixc Bl uIgiyc Bl丁一1。 一2 2
16、J1 0 l 1Ucrbi。 一rti， M u+Cnl移一CRix niz 口 +cniyHc(un。+虢，)一耗；，配+半嚣+c(黼：，+移niy)(13)万方数据第4期 吴泽艳，等：基于简单WENO一间断Galerkin的Euler方程自适应计算809式中：c=厮；BI=(71)c2；B：=B，(M2+F2)2；H=(E+P)p。设需要对“问题单元”。的向量多项式p。(z，Y)进行重构，并保持它在单元内的平均值H(t)。设选择的WENO限制器重构模板为S=。，。，：，各个单元上的DG解分别为po(戈，Y)，Pl(戈，y)，P2(z，Y)，P 3(戈，Y)。首先，在单元。的每一个方向n构造
17、一个新的多项式(P。)?”，分为如下6个步骤：1)修改每一个p，为F。，使得F。在“问题单元”。上具有与P。(戈，Y)相同的单元平均。引)2Pt()一亩协)以dy+衙J。(x,y)也dy扛12，3(14)为了方便，记玩(z，Y)=P。(戈，Y)。将瓦(石，Y)，歹，(戈，Y)，苁(戈，Y)，F，(，y)作投影F：=L。Fi，z=0，1，2，3，每一个F：都是一个七次向量多项式。WENO限制器重构多项式将是4个多项式歹，(戈，y)(=0，1，2，3)的凸组合。2)记线性加权系数为彘、亭。、孝：和邑，因为在光滑区域，4个多项式F。(并，)，)(2-0，1，2，3)都具有矗+1阶精度，所以对于加权系
18、数，首先要求满足亭。十宇。+f：+亭，=1。另外，通过调节加权系数f。、f。、f：和手，来满足解的精度与激波处的本质无振荡的平衡。一般地，取孝。=0997，孝。=孝2=孝3=0001。3)计算光滑指示子卢。，f=0，1，2，3。光滑指示子显式了函数萝，(戈，Y)在“问题单元”。上的光滑度。取光滑指示子为局=f。s|：15 lFi(互，)，)a艽51 a矿2式中：s=(sl，s2)。4)计算基于光滑指示子的非线性加权系数f，Z=0，1，2，3。h=面面，l=0【五。=孝i(8+p；)26=101。(16)小的正数，一般取5)非线性WENO限制器重构多项式(歹。)?“(z，Y)定义为F。(z，Y)
19、(i=0，1，2，3)的凸组合。(瓦)?“(z，Y)=。Fo(z，y)+叫，。(x，Y)+2反(z，y)+甜3F，(x，y) (17)6)将(万。)?“(戈，Y)投影回物理空间得到(p。)?“(z，y)=R。(贰)?“(z，Y) (18)则“问题单元”。上的重构多项式为3 3p飘戈，y)=(p。)：lewJ。J)J (19)最后用重构多项式p：“(x，Y)修改变量即(f)。4 基于自适应网格加密的激波捕捉Euler方程是强非线性方程，即使初值连续，其解也可能出现间断。为了捕捉间断，一般会在计算中引入人工黏性，或者用自带耗散效应的格式(如迎风格式、WENO重构等)来抑制振荡。这时，一般只能在15
20、个网格内捕捉到激波，计算的分辨率比较低。所以，要提高间断的分辨率，就需要非常细的网格。如果在全流场采用一致的细网格，无疑会增加内存消耗和计算量，也意味着增加了计算时间，降低了计算效率。为了有效解决这个问题，依据解的变化来动态调整网格疏密程度是必要的，这就是通常的h自适应方法。这里简单介绍本文采用的h自适应方法。41单元间的数据交换由DG方法的原理可知，单元问的信息是通过单元边界交换传递的，具体表现在单元的边界积分的数值通量需要用到相邻单元的边界信息上。对于不做自适应计算的DG程序而言，只需要通过网格节点标号确定相邻单元的单元编号即可。对于h自适应的DG程序，可以通过单元编号的二叉树方法确定相邻
21、单元编号。本文不采用二又树方法，而是通过边界高斯积分点的坐标来搜索相邻单元，得到相邻单元的单元编号。这种搜索只在网格细化或粗化后做一次，然后保存该结果。因为需要搜索的点总不在计算域的边界上，而非计算域边界上的单元边总可以修改为直线，所以这种搜索是很方便的。如当单元编某边界高斯积分点的坐标为P(石。，y。)时，首先由点P的坐标判断出其可能属于单元门。哦上，然后若(菇。一名。)2+(Yoy。)2+(石。一菇b)2+(YoY6)2一(咒。一戈6)2一(Y。一)，6)210“(20)且(石。，Y。)e。力l哦I(戈6，Y6)e。力l门o 万方数据810 北京航空航天大学学报 2016年则表明单元编上与
22、P点相邻的单元为力。这里(戈。，Y。)和(戈。，Y。)是单元n。120某条边上的端点坐标。即可确定所有单元上全部高斯积分点所对应的相邻单元，从而方便地确定信息传递对象。42加密粗化单元的标识一般而言，h自适应方法根据解的性质来确定局部单元的加密或粗化，即在解急剧变化的地方细化网格，而在解变化较小的地方粗化网格，或者根据解的先验误差估计确定哪些单元细化或粗化。本文旨在提高激波捕捉的分辨率，所以只对“问题单元”进行局部单元加密，对曾是“问题单元”而后来判断不是“问题单元”的单元进行粗化。计算实践表明，这种加密单元的方式是保守的，因为一般情况下根据“问题单元”的判断准则，最终“问题单元”的集合总包含
23、了全部激波附近的单元和和部分非激波附近的单元。43加密粗化过程对四边形单元而言，当确定某单元需要加密时，首先将该单元均分为4个子单元，并计算出各子单元上的正交多项式，然后将母单元上变量的表达式映射到子单元上，最后按第41节所述方法，得到子单元的各个相邻单元。当确定某单元需要粗化时，先检验该单元是否可以粗化，因为只有由细化步骤得到的单元才可以粗化，这只需要检验本单元与相邻单元是否可以构成新的直边四边形即可。如果该单元可以粗化，并得到了向量的3个单元，则由这4个单元就可以粗化为1个新的单元。这时，需要计算出新单元上的正交多项式，并将4个单元上的解用最小二乘法投影到新单元上。5 数值算例51双马赫反
24、射计算域为0，4X0，1，计算域的底部是固壁。初始条件为一右行的激波，其马赫数为10，始于戈=16，Y=0并与戈轴正方向成600夹角。在底部戈0，16的范围内采用激波条件，其他区域采用反射边界条件，计算域的顶部用来流条件。计算网格大小为h=l100。计算时问取为t=02，用3阶精度DG方法计算该问题。图1为不同时刻的密度等值线。不同时刻密度等值线万方数据第4期 吴泽艳，等：基于简单WENO一间断Galerkin的Euler方程自适应计算811从图1可以看出，使用WENO限制器能清晰 14。原始网格及局部加密一次后的网格分别如地捕捉到流场的间断而无明显的数值振荡产生，能够清晰地捕捉到马赫杆结构。
25、52 NACA0012翼型跨声速绕流来流马赫数为08，攻角为125。取来流密度P=10，速度u=08，钟=00，压力P=10O 5 0一O 5一O 5 O 5图2和图3所示。原始网格下八节点四边形单元数为5 000，局部网格加密一次后的单元数为5 597。图4和图5分别为网格局部加密前后的马赫数等值线。图6和图7分别为网格局部加密前后表面压力系数C。曲线。可以看出，简单WENO0 5x0一O 5-0 5 O 5刳2 NACA0012翼型跨声速绕流原始网格(局部) 网3 NACA0012翼型跨声速绕流局部加密网格Fig2 Local、7iew of original mesh for trans
26、(N、(：A001 2 flt)网4 NACA0012翼型跨声速绕流原始网格下的马赫数等值线Fig4 Maeh number contours on original mesh fortransonic NACA0012 flow图6 NACA0012翼型跨声速绕流原始网格下的表面压力系数分布Fig6 Surface pressure coefficient distribution onoriginal mesh for transonic NACA0012 flowFig3 Local view of adaptive mesh tolIransonNACAOOl2 flow图5 NACA
27、0012翼型跨声速绕流网格局部加密后的马赫数等值线Fig5 Mach number contours on adaptive mesh fortransonic NACA0012 flow图7 NACA0012翼型跨声速绕流网格局部加密后的表面压力系数分布Fig7 Surface pressure coefficient distribution onadaptive mesh for transonic NACA0012 flow如丝H魄帅跎鳃鸵拍博m万方数据812 北京航空航天大学学报 2016年限制器能够无波动地捕捉到间断；另外，经过一次局部网格加密后，机翼上下表面激波处马赫数等值线更密
28、集，压力系数曲线斜率更大。53 4凸起计算域及原始网格如图8所示。其中凸起部分为一段圆弧，圆心坐标为(15，一3105)，半径为r=3145。来流马赫数为14，取来流密度P=10，速度“=14，3=00，压力P=1014。计算域的左边为人VI边界条件，右边为出口边界条件，上下为反射边界条件。单元采用八节点四边形单元。局部加密一次后的网格如图9所示。原始网格下单元数为7 500，局部网格加密一次后的单元数为10476。图10和图11分别为网格局部罔8 4凸起原始网格pig8 Full view of OIiginal Icsh foi subsonic n(1past 4bump configm
29、altonJ图9 4凸起局部加密网格Iig9 Full view of adaptive mesh folsubsonic f1()past 4huInp confignialton0 2 j一图10 4凸起原始网格下的密度等值线Fig1 0 Density contours OII original nlesh for subsonic flowpast 41)LllUl)configurationP00 8060 40 20 2 3J图1 1 4凸起局部网格加密后的密度等值线Fig1 1 Density contours on adaptive mesh forsubsonic flow
30、past 4bump configuration加密前后密度等值线。可以看出，简单WENO限制器能够很好地捕捉到间断；另外，经过一次局部单元加密，密度分辨率有了一定的提高。54 2段NACA0012翼型超声速绕流取2个并行错位的NACA0012翼型超声速绕流，翼型的攻角均为0。，来流马赫数为20，取来流密度P=10，速度“=20，口=00，压力P=1014。本算例同时存在内外双重的流动特征，由于翼型错位和超声速来流，翼型之间的干扰所产生的流场波系比较复杂。计算域为一6，65X一65，65。计算域的左边为人口边界条件，上、下、右边为远场边界条件，翼型壁面为反射边界条件。单元采用八节点四边形单元。
31、初始局部网格如图12所示。局部加密一次后的网格如图13所示。初始网格下单元数为41 237，局部网格加密一次后的单元数为45 233。图14和图15分别为网格局部加密前后马赫数等值线。可以看出，下面的翼型未受到干扰，在翼型前缘产生上下对称激波，而上面的激波与下面的翼型产生的激波相交；上面翼型产生的激波的下罔12 2段NACA0012翼型原始网格(局部)Fig1 2 Local view of original mesh for superson ic flowthrough two staggeled NACA0012 airfoils图13 2段NACA0012翼型局部加密网格Fig1 3
32、Local view of adaptive mesh for supersonic flowthrough two staggered NACA00 1 2 airfoils万方数据第4期 吴泽艳，等：基于简单WENO一间断Galerkin的Euler方程白适应计算813图14 2段NACA0012翼型原始网格下的马赫数等值线Fig14 Ma ch number contours on original mesh forsupersonic flow through two staggered NACA0012 airfoils4O 8O2-0 41 O一0 5 0 0 5 1 O 1 5
33、2 0图15 2段NACA0012翼型局部网格加密后的马赫数等值线Fig15 Mach number contours on adaptive mesh forsupersonic flow through two staggered NACA00 1 2 airfoils半支在下面翼型的上表面相交、反射，形成马赫杆结构。另外，简单WENO限制器能够很好地捕捉到间断，经过一次局部单元加密，间断分辨率有了一定的提高。6 结 论本文介绍了间断Galerkin方法、简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法，最后将简单WENO限制器应用到曲边四边形单元上，实现了局部网格加密的自
34、适应计算。通过对若干典型算例进行编程计算，得到如下结论：1)简单WENO限制器可以推广应用到曲边四边形单元上。2)简单WENO限制器可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。参考文献(References)1REED W H，HILL T RTriangular mesh methods for the neutron transport equation：LAUR-73-479RLos Alamos：Scientific Laboratory，19732COCKBURN B，SHU C WThe RungeKutta discontinuousGalerkin meth
35、od for conservation laws V：Multidimensional sys-temsJJournal of Computational Physics1998，141(2)：1992243BASSI F，REBAY SA high order accurate discontinuous finiteelement method for the numerical solution of the compressibleNavier Stokes equationsJJournal of Computational Physics1997，131(2)：267-2794CO
36、CKBURN BSHU C WThe local diseontinuous Galerkinmethod for time-dependent convection diffusion systemsJSIAM Journal on Numerical Analysis，1 998，35(6)：244024635SHU C WA brief survey on discontinuous Galerkin methods incomputationalfluid dynamicsJ力学进展，2013，43(6)：541553SHU C WA brief survey on discontin
37、uous Galerkin methods incomputational fluid dynamicsJAdvances in Mechanics2013，43(6)：541-553(in Chinese)6COCKBURN B，KARNIADAKIS G。SHU c WThe developmenof discontinuous Galerkin methodsMCOCKBURN B，KARNIADAKIS G，SHU C WDiscontinuous Galerkin methods：Theory，computation and applicationsNew York：Springer
38、，2000：1-507COCKBURN B，SHU C WTVB Runge-Kutta local projectiondiscontinuous Galerkin finite element method for eonservationlaws lI：General frameworkJMathematics of Computation1989，52(186)：411-4358BISWAS R，DEVINE K D，FLAHERTY JParalleladaptivefinite element methods for conservation lawsJApplied Nu-mer
39、ical Mathematics，1994，14(sl-3)：255-2839BURBEAU A，SAGAUT P。BRUNEAU C HA problem independent limiter for high order Runge Kutta discontinuous Galerkin methodsJJournal of Computational Physics，200 1，169(1)：11115010SURESH A，HUYNH H TAccurate monotonicity preservingschemes with Runge Kutta time steppingJ
40、Journal of Compu-tational Physics，1997，136(1)：83-99I 1RIDER W J，MARGOLIN L GSimple modifications of monotonicity preserving limitersJJournal of Computational Physics。2001，174(1)：473-4881 2QIU J X，SHU C WRungeKutta discontinuous Galerkin method using WENO limitersJSIAM Journal on Scientific Computing
41、，2005，26(3)：907-92913ADJERID S，DEVINE K，FLAHERTY J，et a1A posteriori er-ror estimation for discontinuous Galerkin solutions of hyperbolicproblemsJComputer Methods in Applied Mechanics andEngineering，2002，191(1112)：10971112万方数据814 北京航空航天大学学报 2016年张来平，刘伟，贺立新，等一种新的间断侦测器及其在DGM中的应用J空气动力学学报，2011，29(4)：401
42、-406ZHANG L P，L1U W，HE L X，et a1A shock detection methodand applications in DGM for hyperbolic conservation laws onunstructured gridsJActa Aerodynamica Sinica，201 1，29(4)：401-406(in Chinese)ZHONG X，SHU C WA simple weighted essentially nonoscillatory limiter for Runge-Kutta discontinuous Galerkin met
43、hodsJJournal of Computational Physics，2013，232(1)：397415作者简介：吴泽艳男，博士，讲师。主要研究方向：计算流体力学。E-mail：wuzeyan2000163corn王立峰男，博士，副教授。主要研究方向：飞行器设计、数值模拟。Tel：010-82339560Email：wanglifengasebuaaeduen武哲男，博士，教授，博士生导师。主要研究方向：飞行器设计、数值模拟。Adaptive simple WENO limiterdiscontinuous Galerkin method forEuler equationsWU Ze
44、yanl，WANG Lifen92一，WU Zhe2(1College of Hydraulic and Environmental Engineering，China Three Gorges University，Yichang 443002，China；2School of Aeronautic Science and Engineering，Beijing University of Aeronautics and Astronautics，Beijing 1 00083，China)Abstract：To achieve high precision and high resolut
45、ion numerical result of Euler equations，the basicprinciple of discontinuous Galerkin method，the simple WENO limiter on triangular meshes and shock capturing method based on adaptive mesh refinement were introducedThe simple WENO limiterdiscontinuousGalerkin method was applied to the curved quadrilat
46、eral element，and the adjacent elements of every elementwith the same coordinates of the Gauss integral points on the boundaries were foundThe adaptive computationbased on“trouble element”refinement was accomplishedSeveral benchmark test cases were computedThenumerical results show that the simple WE
47、NO limiter is appropriate for the curvilinear boundary quadrilateralelement and for the shock capturing based on unstructured grids with hanging nodesKey words：discontinuous Galerkin method；simple WENO limiter；Euler equations；adaptive computation；curvilinear boundary quadrilateral elementReceived：2015-04-18；Accepted：2015-0717：Published online：201510-14 15：18URL：WWWcnkinetkemsdetail112625V201510141518003htmlFoundation item：Talents Scientific Research Starting Foundation of China Three Gorges University(KJ2014803 1)COrresponding authorTel：010-82339560 E-mail：wanglifengasebuaaeduc