A comparison of Rosenbrock and ESDIRK methods combined with iterative solvers for unsteady compressible flows

A comparison of Rosenbrock and ESDIRK methods combined with iterative solvers for unsteady compressible flows

Blom, D.S. (TU Delft Aerodynamics)
Birken, Philipp (Lund University)
Bijl, H. (TU Delft Aerodynamics)
Kessels, F.J.M. (TU Delft Science Centre & Programmering)
Meister, Andreas (University of Kassel)
van Zuijlen, A.H. (TU Delft Aerodynamics)

2016-12-01

In this article, we endeavour to find a fast solver for finite volume discretizations for compressible unsteady viscous flows. Thereby, we concentrate on comparing the efficiency of important classes of time integration schemes, namely time adaptive Rosenbrock, singly diagonally implicit (SDIRK) and explicit first stage singly diagonally implicit Runge-Kutta (ESDIRK) methods. To make the comparison fair, efficient equation system solvers need to be chosen and a smart choice of tolerances is needed. This is determined from the tolerance TOL that steers time adaptivity. For implicit Runge-Kutta methods, the solver is given by preconditioned inexact Jacobian-free Newton-Krylov (JFNK) and for Rosenbrock, it is preconditioned Jacobian-free GMRES. To specify the tolerances in there, we suggest a simple strategy of using TOL/100 that is a good compromise between stability and computational effort. Numerical experiments for different test cases show that the fourth order Rosenbrock method RODASP and the fourth order ESDIRK method ESDIRK4 are best for fine tolerances, with RODASP being the most robust scheme.

ESDIRK
Jacobian-free Newton-Krylov
Navier-Stokes equations
Rosenbrock methods
Time adaptivity
Unsteady flows

http://resolver.tudelft.nl/uuid:8d8e6670-0502-4866-a49b-711e64518da9

https://doi.org/10.1007/s10444-016-9468-x

1019-7168

Advances in Computational Mathematics, 42 (6), 1401-1426

Institutional Repository

journal article

© 2016 D.S. Blom, Philipp Birken, H. Bijl, F.J.M. Kessels, Andreas Meister, A.H. van Zuijlen