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

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 (CAA) 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.

[1]  M. Dumbser,et al.  Arbitrary high order discontinuous Galerkin schemes , 2005 .

[2]  Juhani Pitkäranta,et al.  An analysis of the discontinuous Galerkin method for a scalar hyperbolic equation , 1986 .

[3]  J. Butcher.  The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods , 1987 .

[4]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[5]  Chi-Wang Shu,et al.  A numerical study for the performance of the Runge-Kutta discontinuous Galerkin method based on different numerical fluxes , 2006, J. Comput. Phys..

[6]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[7]  Boleslaw K. Szymanski,et al.  Adaptive Local Refinement with Octree Load Balancing for the Parallel Solution of Three-Dimensional Conservation Laws , 1997, J. Parallel Distributed Comput..

[8]  M. Clemens,et al.  Local timestepping discontinuous Galerkin methods for electromagnetic RF field problems , 2009, 2009 3rd European Conference on Antennas and Propagation.

[9]  J. C. Hardin,et al.  Second Computational Aeroacoustics (CAA) Workshop on Benchmark Problems , 1997 .

[10]  M. Y. Hussaini,et al.  An Analysis of the Discontinuous Galerkin Method for Wave Propagation Problems , 1999 .

[11]  Claes Johnson.  Numerical solution of partial differential equations by the finite element method , 1988 .

[12]  Claus-Dieter Munz,et al.  A Discontinuous Galerkin Scheme based on a Space-Time Expansion II. Viscous Flow Equations in Multi Dimensions , 2008, J. Sci. Comput..

[13]  Eleuterio F. Toro,et al.  ADER schemes for three-dimensional non-linear hyperbolic systems , 2005 .

[14]  Christophe Eric Corre,et al.  Multi-Size-Mesh, Multi-Time-Step Algorithm for Noise Computation Around an Airfoil in Curvilinear Meshes , 2007 .

[15]  P. Rentrop,et al.  Multirate Partitioned Runge-Kutta Methods , 2001 .

[16]  Min Jiang,et al.  A MULTI-TIME-STEP STRATEGY BASED ON AN OPTIMIZED TIME INTERPOLATION SCHEME FOR OVERSET GRIDS , 2010 .

[17]  Willem Hundsdorfer,et al.  A multirate time stepping strategy for stiff ordinary differential equations , 2007 .

[18]  P. Raviart,et al.  On a Finite Element Method for Solving the Neutron Transport Equation , 1974 .

[19]  Michael Dumbser,et al.  Building Blocks for Arbitrary High Order Discontinuous Galerkin Schemes , 2006, J. Sci. Comput..

[20]  Harold L. Atkins,et al.  Two-dimensional Wave Analysis of the Discontinuous Galerkin Method with Non-Uniform Grids and Boundary Conditions , 2002 .

[21]  Claus-Dieter Munz,et al.  A Discontinuous Galerkin Scheme Based on a Space–Time Expansion. I. Inviscid Compressible Flow in One Space Dimension , 2007, J. Sci. Comput..

[22]  Jan S. Hesthaven,et al.  Application of implicit-explicit high order Runge-Kutta methods to discontinuous-Galerkin schemes , 2007, J. Comput. Phys..

[23]  M. Dumbser,et al.  Heterogeneous Domain Decomposition for Computational Aeroacoustics , 2006 .

[24]  Manuel Calvo,et al.  Short note: a new minimum storage Runge-Kutta scheme for computational acoustics , 2004 .

[25]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[26]  Michael Dumbser,et al.  Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems , 2007, J. Comput. Phys..

[27]  Christopher K. W. Tam,et al.  Multi-size-mesh Multi-time-step Dispersion-relation-preserving Scheme for Multiple-scales Aeroacoustics Problems , 2002 .

[28]  Christopher K. W. Tam,et al.  Multi-size-mesh Multi-time-step Dispersion-relation-preserving Scheme for Multiple-scales Aeroacoustics Problems , 2003 .

[29]  C. W. Gear,et al.  Multirate linear multistep methods , 1984 .

[30]  E. Toro,et al.  An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes - V. Local time stepping and p-adaptivity , 2007 .

[31]  Michael Dumbser,et al.  The discontinuous Galerkin method with Lax-Wendroff type time discretizations , 2005 .

[32]  E. Sonnendrücker,et al.  Numerical Methods for Hyperbolic and Kinetic Problems , 2005 .

[33]  Holger Waubke,et al.  A SIMPLE METHOD FOR COMPUTING THE GREEN FUNCTIONS OF THE HELMHOLTZ EQUATION IN THE TWO-DIMENSIONAL IMPEDANCE HALF SPACE , 2010 .

[34]  M. Carpenter,et al.  Fourth-order 2N-storage Runge-Kutta schemes , 1994 .

[35]  Fang Q. Hu,et al.  A Perfectly Matched Layer absorbing boundary condition for linearized Euler equations with a non-uniform mean flow , 2005 .

[36]  Claus-Dieter Munz,et al.  Direct aeroacoustic simulations based on domain decompositions , 2008 .

[37]  Todd E. Peterson,et al.  A note on the convergence of the discontinuous Galerkin method for a scalar hyperbolic equation , 1991 .

[38]  Mark Ainsworth,et al.  Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods , 2004 .

[39]  Harold L. Atkins,et al.  Eigensolution analysis of the discontinuous Galerkin method with non-uniform grids , 2001 .

[40]  W. Habashi,et al.  2N-Storage Low Dissipation and Dispersion Runge-Kutta Schemes for Computational Acoustics , 1998 .

[41]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[42]  C. Bogey,et al.  Low-dissipation and low-dispersion fourth-order Runge–Kutta algorithm , 2006 .

[43]  M.Y. Hussaini,et al.  Low-Dissipation and Low-Dispersion Runge-Kutta Schemes for Computational Acoustics , 1994 .

[44]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[45]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..