Nonuniform time-step Runge-Kutta discontinuous Galerkin method for Computational Aeroacoustics

被引:40
作者
Liu, Li [1 ]
Li, Xiaodong [1 ]
Hu, Fang Q. [2 ]
机构
[1] Beihang Univ, Sch Jet Propuls, Beijing 100191, Peoples R China
[2] Old Dominion Univ, Dept Math & Stat, Norfolk, VA 23529 USA
基金
美国国家科学基金会; 中国国家自然科学基金;
关键词
Runge-Kutta method; Discontinuous Galerkin method; Nonuniform time-step size; Computational Aeroacoustics; SCALAR HYPERBOLIC EQUATION; LOW-DISSIPATION; UNSTRUCTURED MESHES; CONSERVATION-LAWS; SCHEMES; SYSTEMS; FLOW; ACOUSTICS; EXPANSION; DIMENSION;
D O I
10.1016/j.jcp.2010.05.028
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
With many superior features, Runge-Kutta discontinuous Galerkin method (RKDG), which adopts Discontinuous Galerkin method (DG) for space discretization and Runge-Kutta method (RK) for time integration, has been an attractive alternative to the finite difference based high-order Computational Aeroacoustics (CM) approaches. However, when it comes to complex physical problems, especially the ones involving irregular geometries, the time step size of an explicit RK scheme is limited by the smallest grid size in the computational domain, demanding a high computational cost for obtaining time accurate numerical solutions in CAA. For computational efficiency, high-order RK method with nonuniform time step sizes on nonuniform meshes is developed in this paper. In order to ensure correct communication of solutions on the interfaces of grids with different time step sizes, the values at intermediate-stages of the Runge-Kutta time integration on the elements neighboring such interfaces are coupled with minimal dissipation and dispersion errors. Based upon the general form of an explicit p-stage RK scheme, a linear coupling procedure is proposed, with details on the coefficient matrices and execution steps at common time-levels and intermediate time-levels. Applications of the coupling procedures to Runge-Kutta schemes frequently used in simulation of fluid flow and acoustics are given, including the third-order TVD scheme, and low-storage low dissipation and low dispersion (LDDRK) schemes. In addition, an analysis on the stability of coupling procedures on a nonuniform grid is carried out. For validation, numerical experiments on one-dimensional and two-dimensional problems are presented to illustrate the stability and accuracy of proposed nonuniform time-step RKDG scheme, as well as the computational benefits it brings. Application to a one-dimensional nonlinear problem is also investigated. (C) 2010 Elsevier Inc. All rights reserved.
引用
收藏
页码:6874 / 6897
页数:24
相关论文
共 39 条
  • [1] Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods
    Ainsworth, M
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 2004, 198 (01) : 106 - 130
  • [2] [Anonymous], 1997, 2 COMP AER CAA WORKS
  • [3] [Anonymous], 1974, PUBLICATIONS MATH IN
  • [4] Low-dissipation and low-dispersion fourth-order Runge-Kutta algorithm
    Berland, Julien
    Bogey, Christophe
    Bailly, Christophe
    [J]. COMPUTERS & FLUIDS, 2006, 35 (10) : 1459 - 1463
  • [5] Butcher JC., 1987, The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods
  • [6] A new minimum storage Runge-Kutta scheme for computational acoustics
    Calvo, M
    Franco, JM
    Rández, L
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 2004, 201 (01) : 1 - 12
  • [7] Carpenter M. H., 1994, 109112 NASA
  • [8] Runge-Kutta discontinuous Galerkin methods for convection-dominated problems
    Cockburn, Bernardo
    Shu, Chi-Wang
    [J]. Journal of Scientific Computing, 2001, 16 (03) : 173 - 261
  • [9] The Runge-Kutta discontinuous Galerkin method for conservation laws V - Multidimensional systems
    Cockburn, B
    Shu, CW
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 141 (02) : 199 - 224
  • [10] Arbitrary high order discontinuous Galerkin schemes
    Dumbser, M
    Munz, CD
    [J]. NUMERICAL METHODS FOR HYPERBOLIC AND KINETIC PROBLEMS, 2005, 7 : 295 - 333