Optimized Runge-Kutta (LDDRK) timestepping schemes for non-constant-amplitude oscillations

by   Aldaïr Petronilia, et al.

Finite differences and Runge-Kutta time stepping schemes used in Computational AeroAcoustics simulations are often optimized for low dispersion and dissipation (e.g. DRP or LDDRK schemes) when applied to linear problems in order to accurately simulate waves with the least computational cost. Here, the performance of optimized Runge-Kutta time stepping schemes for linear time-invariant problems with non-constant-amplitude oscillations is considered. This is in part motivated by the recent suggestion that optimized spatial derivatives perform poorly for growing and decaying waves, as their optimization implicitly assumes real wavenumbers. To our knowledge, this is the first time the time-stepping of non-constant-amplitude oscillations has been considered. It is found that current optimized Runge-Kutta schemes perform poorly in comparison with their maximal order equivalents for non-constant-amplitude oscillations. Moreover, significantly more accurate results can be achieved for the same computation cost by replacing a two-step scheme such as LDDRK56 with a single step higher-order scheme with a longer time step. Attempts are made at finding optimized schemes that perform well for non-constant-amplitude oscillations, and three such examples are provided. However, the traditional maximal order Runge-Kutta time stepping schemes are still found to be preferable for general problems with broadband excitation. These theoretical predictions are illustrated using a realistic 1D wave-propagation example.


page 4

page 5

page 6

page 7

page 8

page 10

page 11


Efficient Explicit Time Stepping of High Order Discontinuous Galerkin Schemes for Waves

This work presents algorithms for the efficient implementation of discon...

Spline-oriented inter/extrapolation-based multirate schemes of higher order

Multirate integration uses different time step sizes for different compo...

Rosenbrock-Wanner Time Integration in Atmospheric Modelling

Non-hydrostatic atmospheric models often use semi-implicit temporal disc...

A complete and consistent second-order hydrodynamic model for floating structures with large horizontal motions

Floating offshore structures often exhibit low-frequency oscillatory mot...

Transmission operators for the non-overlapping Schwarz method for solving Helmholtz problems in rectangular cavities

In this paper we discuss different transmission operators for the non-ov...

Stable explicit schemes for simulation of nonlinear moisture transfer in porous materials

Implicit schemes have been extensively used in building physics to compu...

Convergence of optimized non-overlapping Schwarz method for Helmholtz problems in closed domains

In this paper we discuss the convergence of state-of-the-art optimized S...

Please sign up or login with your details

Forgot password? Click here to reset