Assessment of high-order IMEX methods for incompressible flow

12/08/2021
by   Montadhar Guesmi, et al.
0

This paper investigates the competitiveness of semi-implicit Runge-Kutta (RK) and spectral deferred correction (SDC) time-integration methods up to order six for incompressible Navier-Stokes problems in conjunction with a high-order discontinuous Galerkin method for space discretization. It is proposed to harness the implicit and explicit RK parts as a partitioned scheme, which provides a natural basis for the underlying projection scheme and yields a straight-forward approach for accommodating nonlinear viscosity. Numerical experiments on laminar flow, variable viscosity and transition to turbulence are carried out to assess accuracy, convergence and computational efficiency. Although the methods of order 3 or higher are susceptible to order reduction due to time-dependent boundary conditions, two third-order RK methods are identified that perform well in all test cases and clearly surpass all second-order schemes including the popular extrapolated backward difference method. The considered SDC methods are more accurate than the RK methods, but become competitive only for relative errors smaller than ca 10^-5.

READ FULL TEXT
research
01/31/2020

A spectral deferred correction method for incompressible flow with variable viscosity

This paper presents a semi-implicit spectral deferred correction (SDC) m...
research
04/28/2022

High-Order Discretization of Backward Anisotropic Diffusion and Application to Image Processing

Anisotropic diffusion is a well recognized tool in digital image process...
research
11/14/2019

Fast multigrid solution of high-order accurate multi-phase Stokes problems

A fast multigrid solver is presented for high-order accurate Stokes prob...
research
12/31/2017

GPU Acceleration of a High-Order Discontinuous Galerkin Incompressible Flow Solver

We present a GPU-accelerated version of a high-order discontinuous Galer...

Please sign up or login with your details

Forgot password? Click here to reset