Analysis for centered nonoscillatory scheme of third order
-
摘要: 三阶中心无振荡格式Cn3(Centered nonoscillatory scheme of third order)使用对称模板构造具有三阶精度的插值公式.利用单调区域与精确区域修正原始插值参数,该格式能够获得间断附近无振荡、光滑区域高精度的计算结果.通过一维和二维典型算例,将Cn3格式与三阶和五阶WENO格式(Weighted Essentially Non-Oscillatory schemes)进行比较,重点分析了3种格式的间断分辨率、计算稳定性和数值耗散性.分析可见,Cn3格式能够精确、稳定地捕捉激波和接触间断,同时对光滑流动区域的小尺度流动结构保持较低的耗散,值得进一步研究及推广应用.Abstract: Choosing symmetric stencils, a centered nonoscillatory scheme of third order scheme (Cn3) constructs a formula which achieves third-order accuracy. A modification of the original estimates, in the reconstruction procedure, was undertaken by employing a monotonicity region as well as an accuracy one. After the very modification, the Cn3 scheme not only obtained nonoscillatory results near discontinuities but also achieved high-accuracy calculation in smooth regions. By using several typical one- and two-dimensional test cases, Cn3 scheme was compared with the weighted essentially non-oscillatory (WENO) schemes of third and fifth order. The properties of capturing discontinuities, stability/robustness and numerical dissipation were significantly considered. The results of the numerical experiments confirm that Cn3 scheme has the ability to suppress spurious numerical oscillations near shocks and contact discontinuities when sharply capturing them, which indicates its characteristics of stability and accuracy, and has low dissipation in smooth region at the same time. Cn3 scheme is worth of further study and application.
-
Key words:
- high-accuracy scheme /
- WENO schemes /
- discontinuities capturing /
- stability /
- numerical dissipation
-
[1] Kroll N.ADIGMA-A European project on the development of adaptive higher-order variational methods for aerospace applications[C]//47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition.Reston,VA:AIAA,2009 [2] 阎超,于剑,徐晶磊,等.CFD模拟方法的发展成就与展望[J].力学进展,2011,41(5):562-589 Yan Chao,Yu Jian,Xu Jinglei,et al.On the achievements and prospects of the methods for computational fluid dynamics[J].Advances in Mechanics,2011,41(5):562-589(in Chinese) [3] Harten A,Engquist B,Osher S,et al.Uniformly high order accurate essentially non-oscillatory schemes,III[J].Journal of Computational Physics,1987,71(2):231-303 [4] Shu C W,Osher S.Efficient implementation of essentially non-oscillatory shock-capturing schemes[J].Journal of Computational Physics,1988,77(2):439-471 [5] Osher S,Shu C.Efficient implementation of essentially non-oscillatory shock-capturing schemes,II[J].Journal of Computational Physics,1989,83(1):32-78 [6] Liu X D,Osher S,Chan T.Weighted essentially non-oscillatory schemes[J].Journal of Computational Physics,1994,115(1):200-212 [7] Jiang G S,Shu C W.Efficient implementation of weighted ENO schemes[J].Journal of Computational Physics,1996,126(1):202-228 [8] Shu C W.High order weighted essentially nonoscillatory schemes for convection dominated problems[J].SIAM Review,2009,51(1):82-126 [9] Ivan L,Groth C P T.High-order solution-adaptive central essentially non-oscillatory (CENO) method for viscous flows[J].AIAA Paper,2011:2011-367 [10] Huynh H T.Accurate upwind schemes for the Euler equations[C]//12th Computational Fluid Dynamics Conference.San Diego,CA:AIAA,1995:1022-1039 [11] Suresh A.Centered nonoscillatory schemes of third order[C]//12th Computational Fluid Dynamics Conference.San Diego,CA:AIAA,1995:1292-1305 [12] Suresh A,Huynh H T.Accurate monotonicity-preserving schemes with Runge-Kutta time stepping[J].Journal of Computational Physics,1997,136(1):83-99 [13] Sod G A.A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws[J].Journal of Computational Physics,1978,27(1):1-31 [14] Nichols R H,Tramel R W,Buning P G.Evaluation of two high-order weighted essentially nonoscillatory schemes[J].AIAA Journal,2008,46(12):3090-3102 [15] Woodward P,Colella P.The numerical simulation of two-dimensional fluid flow with strong shocks[J].Journal of Computational Physics,1984,54(1):115-173
点击查看大图
计量
- 文章访问数: 831
- HTML全文浏览量: 97
- PDF下载量: 652
- 被引次数: 0