A robust implicit high-order discontinuous Galerkin method for solving compressible Navier-Stokes equations on arbitrary grids
严佳 1杨小权 1翁培奋2
扫码查看
点击上方二维码区域，可以放大扫码查看
作者信息
1. Shanghai Institute of Applied Mathematics and Mechanics,Shanghai Key Laboratory of Mechanics in Energy Engineering,School of Mechanics and Engineering Science,Shanghai University,Shanghai 200072,China
2. Shanghai Institute of Applied Mathematics and Mechanics,Shanghai Key Laboratory of Mechanics in Energy Engineering,School of Mechanics and Engineering Science,Shanghai University,Shanghai 200072,China;College of Energy and Mechanical Engineering,Shanghai University of Electric Power Shanghai 200090,China
The primary impediments impeding the implementation of high-order methods in simulating viscous flow over complex configurations are robustness and convergence.These challenges impose significant constraints on computational efficiency,particularly in the domain of engineering applications.To address these concems,this paper proposes a robust implicit high-order discontinuous Galerkin(DG)method for solving compressible Navier-Stokes(NS)equations on arbitrary grids.The method achieves a favorable equilibrium between computational stability and efficiency.To solve the linear system,an exact Jacobian matrix solving strategy is employed for preconditioning and matrix-vector generation in the generalized minimal residual(GMRES)method.This approach mitigates numerical errors in Jacobian solution during implicit calculations and facilitates the implementation of an adaptive Courant-Friedrichs-Lewy(CFL)number increasing strategy,with the aim of improving convergence and robustness.To further enhance the applicability of the proposed method for intricate grid distortions,all simulations are performed in the reference domain.This practice significantly improves the reversibility of the mass matrix in implicit calculations.A comprehensive analysis of various parameters influencing computational stability and efficiency is conducted,including CFL number,Krylov subspace size,and GMRES convergence criteria.The computed results from a series of numerical test cases demonstrate the promising results achieved by combining the DG method,GMRES solver,exact Jacobian matrix,adaptive CFL number,and reference domain calculations in terms of robustness,convergence,and accuracy.These analysis results can serve as a reference for implicit computation in high-order calcula-tions.
万方数据