摘要
双曲守恒型方程的高精度、高分辨率计算格式的研究一直是计算流体力学的热点问题。针对原WENO-JS格式分辨率较低和计算量偏大的不足问题,提出利用简单的重构数值通量的方法以提高计算效率,构造了新的简单限制器的5阶迎风型WENO格式。通过MATLAB软件的仿真对Lax-Wendroff WENO-JS格式、Lax-Wendroff简单限制器WENO格式、Runge-Kutta WENO-JS格式、Runge-Kutta简单限制器的WENO格式的实验结果进行了分析,并比较了这四种计算格式的计算效率和计算精度。数值实验表明:新格式Lax-Wendroff简单限制器WENO格式在保持原WENO分辨率的前提下,计算速度有明显提高,减少了20%的计算时间。
The research of high accuracy and high resolution schemes have been a hot topic in computational mathematics. According to low resolution and large amount of calculation of the original WENO-JS scheme,we propose a simple new limiter fifth order upwind WENO scheme to reconstruct the numerical flux of the simple structure to improve the computational efficiency. Compared with other efficient high accuracy schemes such as ENO and WENO,it is shown that the computational cost of this scheme is less than that of WENO-JS in the same accuracy. By use of MATLAB software,we compared and analyzed computational efficiencies and computational accuracies of Lax-Wendroff WENO-JS scheme,Lax-Wendroff simple limiter WENO scheme,Runge-Kutta simple limiter WENO scheme and Runge-Kutta WENO-JS scheme. The numerical results show that the new Lax-Wendroff simple limiter WENO scheme can improve the computing speed and reduce the computing time by 20% while maintaining the original WENO resolution.
引文
[1]陆金甫.偏微分方程数值解法[M].北京:清华大学出版社,2004.
[2]袁光伟,杭旭登,盛志强,等.福射扩散计算方法若干研究进展[M].计算物理,2009,26(4):475-500.
[3]欧阳颀.非线性科学与斑图动力学导论[M].北京:北京大学出版社,2010.
[4]刘儒勋,舒其望.计算流体力学的若干新方法[M].北京:科学出版社,2003.
[5]WARMING R F.Flux Vector Splitting of the Inviscid Gas-Dynamic Equation with Application of Finite Difference Method[J].Journal of Computational,1981,40:263-293.
[6]GODONOV S K.A Difference Schemes for Numerical Computation of Discontinuous Solution of Hydrodynamics Equations[J].Math.Sbornik,1959,47:271-306.
[7]ENGQUIST,OSHER S.Stable and Entropy Satisfying Approximation for Transonic Flow Calculation[J].Mathematics of Computation,1980,34:45-75.
[8]COCKBURN B,LIN S Y,SHU C W.TVB Runge-Kutta Local Projection Discontinuous Galerkin Finite Element Method for Conservation Laws III:One-Dimensional Systems[J].J.Comput.Phys.,1989,84:90-113.
[9]COCKBURN B,HOU S,SHU C W.TVB Runge-Kutta Local Projection Discontinuous Galerkin Finite element method for conservation laws IV[J].Math.Comput.,1990,54:545-581.
[10]LIU H,XU H,GONG H.Modeling the Asymmetry in Traffic Flow(b):Macroscopic Approach[J].Appl.Math.Model,2013,37(22):9441-9450.
[11]ZHU J,QIU J X.A New Fifth Order Finite Difference WENO Scheme for Solving Hyperbolic Conservation Laws[J].J.Comput.Phys.,2016,318:110-121.
[12]LIU X,ZHANG S,ZHANG H,SHU C W.A New Class of Central Compact Schemes with Spectral-like Resolution II:Hybrid Weighted Nonlinear Schemes[J].J.Comput.Phys.,2015,284:133-154.
[13]PIROZZOLI S.Conservative Hybrid Compact-WENO Schemes for Shock-Turbulence Interaction[J].J.Comput.Phys.,2002,178(1):81-117.
[14]QIU J X,DUMBSER M,SHU C W.The Discontinuous Galerkin Method with Lax-Wendroff type Time Discretizations[J].Comput.Methods Appl.Mech.Engrg.2005,194:4528-4543.
[15]LAX P D.Hyperbolic System of Conservation Laws II[J].Comm.Pure and Applied Mathematics.1957,10:537-566.
[16]ZHANG X,SHU C W.On Maximum-Principle-Satisfying High Order Schemes for Scalar Conservation Laws[J].J.Comput.Phys.,2010,229(9):3091-3120.
[17]ZHANG X,SHU C W.Maximum-Principle-Satisfying and Positivity-Preserving High-Order Schemes for Conservation Laws:Survey and New Developments Proc.R.Soc.A[J].Math.Phys.Eng.Sci.,2011,467(2):2752-2776.
[18]LIU Y J,SHU C W,TADMOR E,ZHANG M P.Central Discontinuous Galerkin Methods on Overlapping Cells with A Non-Oscillatory Hierarchical Reconstruction[J].SIAM J.Numer.Anal.,2007,45:2442-2467.
[19]姚慧丽,宋晓秋,李兴华.一类半线性微分方程的渐进概自守温和解[J].哈尔滨理工大学学报,2012,17(1):72-78.
[20]姚慧丽,卜宪江,宋晓秋.一类微分方程的指数增长的温和渐进概自守解[J].哈尔滨理工大学学报,2014,19(5):23-26.
