Efficient fully discrete summation-by-parts schemes for unsteady flow problems

We make an initial investigation into the temporal efficiency of a fully discrete summation-by-parts approach for unsteady flows. As a model problem for the Navier–Stokes equations we consider a two-dimensional advection–diffusion problem with a boundary layer. The problem is discretized in space using finite difference approximations on summation-by-parts form together with weak boundary conditions, leading to optimal stability estimates. For the time integration part we consider various forms of high order summation-by-parts operators and compare with an existing popular fourth order diagonally implicit Runge–Kutta method. To solve the resulting fully discrete equation system, we employ a multi-grid scheme with dual time stepping.