基于Newton/Gauss-Seidel迭代的DGM隐式方法  

AN IMPLICIT ALGORITHM FOR DISCONTINUOUS GALERKIN METHOD BASED ON NEWTON/GAUSS-SEIDEL ITERATIONS

作　　者：刘伟[1] 张来平[1,2] 赫新[2] 贺立新[2] 张涵信[2] 

机构地区：[1]中国空气动力研究与发展中心空气动力学国家重点实验室,绵阳621000 [2]中国空气动力研究与发展中心计算空气动力学研究所,绵阳621000

出　　处：《力学学报》2012年第4期792-796,共5页Chinese Journal of Theoretical and Applied Mechanics

基　　金：国家重点基础研究发展计划(2009CB723802);国家自然科学基金(11028205;91016011;91130029);空气动力学国家重点实验室基金(JBKY11010913)资助项目~~

摘　　要：在Newton迭代方法的基础上,对高阶精度间断Galerkin有限元方法 (DGM)的时间隐式格式进行了研究.Newton迭代法的优势在于收敛效率高效,并且定常和非定常问题能够统一处理,对于非定常问题无需引入双时间步策略.为了避免大型矩阵的求逆,采用一步Gauss-Seidel迭代和Matrix-free技术消去残值Jacobi矩阵的上、下三角矩阵,从而只需计算和存储对角(块)矩阵.对角(块)矩阵采用数值方法计算.空间离散采用Taylor基,其优势在于对于任意形状的网格,基函数的形式是一致的,有利于在混合网格上推广.利用该方法,数值模拟了Bump绕流和NACA0012翼型绕流.计算结果表明,与显式的Runge-Kutta时间格式相比,隐式格式所需的迭代步数和CPU时间均在很大程度上得到减少,计算效率能够提高1～2个量级.An efficient implicit algorithm was developed for high-order discontinuous Galerkin （DGM） based on Newton/Gauss-Seidel iteration approach. The second-order to the forth order DGMs based on Taylor basis functions were employed to carry out the spatial discretization. Newton iteration scheme was used to solve the nonlinear syStem, and the linear system was solved with one-step Gauss-Seidel iteration. In addition, the effects of several parameters in the implicit scheme, such as the CFL number, the Newton sub-iteration steps, and the update frequency of mass-matrix, have been investigated for two-dimensional Euler equations. Two typical cases, including subsonic flows over a bump and a NACA0012 airfoil, were simulated, and compared with the traditional explicit Runge-Kutta scheme. The numerical results demonstrate that the present implicit scheme can accelerate the convergence history evidently.

关 键 词：间断Galerkin有限元 Taylor基函数 NEWTON迭代 Gauss-Seidel迭代 时间隐式方法 

分 类 号：V211.3[航空宇航科学与技术—航空宇航推进理论与工程]