基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用

龚小权; 贾洪印; 陈江涛; 赵辉; 周桂宇

doi:10.7638/kqdlxxb-2018.0189

欢迎来到能源电力期刊网！

English Version 期刊登录

您当前的位置：

首页 >

文章列表页 >

基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用

专栏——首届中国空气动力学大会青年优秀论文 | 更新时间：2026-05-28

- 基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用
- Applications of GMRES based on exact calculations of Jacobian matrix in discontinuous Galerkin methods
- 空气动力学学报 2019年37卷第1期页码：121-132
- 作者机构：
  
  中国空气动力研究与发展中心,四川,绵阳,621000
- 作者简介：
  
  龚小权(1984-), 男, 重庆潼南人, 助理研究员, 硕士研究生, 研究方向:高阶精度DG方法及亚跨超声速计算空气动力学.E-mail:xiaoquan19492002@163.com
  贾洪印, 男, E-mail:hongyinjia@foxmail.com
- 基金信息：
- DOI：10.7638/kqdlxxb-2018.0189
  中图分类号： V211.3;V211.78
- 纸质出版：2019
- 稿件说明：
移动端阅览
龚小权, 贾洪印, 陈江涛, 赵辉, 周桂宇. 基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用[J]. 空气动力学学报, 2019, 37(1): 121-132.

GONG Xiaoquan, JIA Hongyin, CHENG Jiangtao, et al. Applications of GMRES based on exact calculations of Jacobian matrix in discontinuous Galerkin methods[J]. 2019, 37(1): 121-132.
龚小权, 贾洪印, 陈江涛, 赵辉, 周桂宇. 基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用[J]. 空气动力学学报, 2019, 37(1): 121-132. DOI： 10.7638/kqdlxxb-2018.0189.

GONG Xiaoquan, JIA Hongyin, CHENG Jiangtao, et al. Applications of GMRES based on exact calculations of Jacobian matrix in discontinuous Galerkin methods[J]. 2019, 37(1): 121-132. DOI： 10.7638/kqdlxxb-2018.0189.

摘要

为改善高阶间断Galerkin有限元方法（DG）时间推进效率，在三维非结构网格下针对该方法建立了并行广义最小残差（Generalized Minimal Residual，GMRES）隐式时间迭代方法，GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率，发展了方程组右端项残值雅可比精确计算方法，针对无黏通量Roe格式和黏性通量BR2（Bassi Rebay 2）黏性计算方法，分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响，然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率，同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS（Lower Upper-Symmetric Gauss-Seidel）的计算效率。结果表明，建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL（Courant-Friedrichs-Lewy）数，相比前面提到的其它方法具有更高的计算效率，其收敛速度实现量级以上的提高。

Abstract

To improve iteration efficiency of high-order discontinuous Galerkin finite element (DG) method

an efficient implicit generalized minimal residual (GMRES) method has been applied to DG method based on unstructured grids. The method is implemented by the ksp solver of portable

extensible toolkit for scientific computation (PETSc) library. In order to further improve computational efficiency of GMRES

the exact calculation method of Jacobian matrix is developed. The Jacobian matrix of Roe scheme and bassi rebay 2 (BR2) scheme are exactly obtained for the freedom of conservation variables. The developed method is applied to compute a variety of steady inviscid and viscous flow problems. First

a NACA0012 airfoil is used to study the effect of restarted iteration and convergence parameters on the GMRES convergence. The computational efficiencies of the GMRES based on different calculation methods of Jacobian matrix are compared by the inviscid and viscous examples. The computational efficiencies are also compared for the GMRES and LU-SGS iterative method. The results show that the exact calculation method of Jacobian matrix can greatly increase the CFL number of the GMRES for different DG methods. Compared with other methods

the GMRES based on the exact calculation method of Jacobian matrix has higher computational efficiency

and its convergence speed is improved by more than one order of magnitude.

关键词

Keywords

references

COCKBURN B, SHU C W. The Runge-Kutta discontinuous Galerkin method for conservation laws V:multidimensional systems[J]. Journal of Computational Physics, 1998, 141:99-224.

REED W H, HILL T R. Triangular mesh methods for the neutron transport equation[R]. Los Alamos Scientific Laboratory Report, LAUR-73-479, 1973.

MAM S, GONG X Q, DENG Y Q, et al. An implicit LU-SGS scheme for the discontinuous Galerkin method on unstructured grids[J]. Journal of Northwestern Polytechnical University, 2016, 34(4):754-760. (in Chinese)

GUO Y H, YANG Y, ZHANG Q. A high efficient implicit discontinuous Galerkin method[J]. Acta Aerodynamica Sinica, 2012, 30(2):230-233. (in Chinese)

TAKANORI H, KEISUKE S, WANG Z J. Animplicit LU-SGS scheme for the spectral volume method on unstructured tetrahedral grids[J]. Communications in Computational Physics, 2009, 6(5):978-996.

LIU W, ZHANG L P, HE X, et al. An implicit algorithm for discontinuous Galerkin method based on Newton/Gauss-seidel iterations[J]. Chinese Journal of Theoretical and Applied Mechanics, 2012, 44(4):792-796. (in Chinese)

PERSSON P O, PERAIRE J. Newton-GMRES preconditioning for discontinuous Galerkin discretizations of the Navier-Stokes equations[J]. SIAM Journal on Scientific Computing, 2008, 30(6):2709-2733.

CRIVELLINI A, BASSI F. An implicit matrix-free discontinuous Galerkin solver for viscous and turbulent aerodynamic simulations[J]. Computers Fluids, 2011, 50:81-93.

LASLO T D, DAVID L D. Discontinuous Galerkin solutions of the Navier-Stokes equations using linear multigrid preconditioning[R]. AIAA 2007-3942.

CHEN S Q, ZHANG T T. Krylov implicit integration factor methods for spatial discretization on high dimensional unstructured meshes:application to discontinuous Galerkin methods[J]. Journal of Computational Physics, 2011, 230:4336-4352.

LANDMANN B. A parallel discontinuous Galerkin code for the Navier-Stokes and Reynolds-averaged Navier-Stokes equations[D]. Stuttgart, university of Stuttgart, 2008.

DANG Y B, LIU K L, SUN Y F, et al. Implicit high order discontinuous Galerkin method for compressible laminar and turbulent flow simulation[J]. Acta Aerodynamica Sinica, 2018, 36(3):535-541. (in Chinese)

YANG X Q, CHENG J, WANG C J, et al. A fast, implicit discontinuous Galerkin method based on analytical Jacobians for the compressible Navier-Stokes equations[R]. AIAA 2016-1326.

CHENG J, LIU X D, YANG X Q, et al. A direct discontinuous Galerkin method for the compressible Navier-Stokes equations on arbitrary grids[R]. AIAA 2016-3334.

BALAY S, ABHYANKAR S, ADAMS M, et al. PETSc users manual[OL ] . http://www.mcs.anl.gov/petsc/documentation/index.html .

BLAZEK J. Computational fluid dynamics: principles and applications[M]. London Elsevier Science, 2001.

ROE P L. Approximate riemann solvers, parameter vectors, and difference schemes[J]. Journal of Computational Physics, 1981, 43:357-372.

BASSI F, REBAY S, MARIOTTI G, et al. A high-order accurate discontinuous finite element method for inviscid and viscous turbomachinery flows[C]//In Proceedings of the 2nd European Conference on Turbomachinery Fluid Dynamics and Thermodynamics, March 5-71997, 99-108.

KRIST S L, BIEDRON R T, RUMSEY C L. CFL3D User's Manual (Version 5.0)[EB]. NASA/TM-1998-208444.

SCHMITT V, CHARPIN F. Pressure distributions on the ONERA-M6 wing at transonic Mach numbers[R]. Experimental Data Base for Computer Program Assessment, AGARD AR138, 1979.

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

适用于非定常流模拟的分布式并行GMRES方法

基于Ac-准则的最优输入设计方法

水滴撞击运动超疏水表面特性及其调控

离心通风机降噪研究进展

热防护材料动态烧蚀形貌结构光视觉测量技术