Comparison of ROW, ESDIRK, and BDF2 for Unsteady Flows with the High-Order Flux Reconstruction Formulation

Author/Creator ORCID

Date

2020-05-13

Department

Program

Citation of Original Publication

Wang, L., Yu, M. Comparison of ROW, ESDIRK, and BDF2 for Unsteady Flows with the High-Order Flux Reconstruction Formulation. J Sci Comput 83, 39 (2020). https://doi.org/10.1007/s10915-020-01222-z

Rights

This item is likely protected under Title 17 of the U.S. Copyright Law. Unless on a Creative Commons license, for uses protected by Copyright Law, contact the copyright holder or the author.
This is a post-peer-review, pre-copyedit version of an article published in Journal of Scientific Computing. The final authenticated version is available online at: http://dx.doi.org/10.1007/s10915-020-01222-z.
Access to this item will begin on 5/13/21

Abstract

We conduct a comparative study of the Jacobian-free linearly implicit Rosenbrock–Wanner (ROW) methods, the explicit first stage, singly diagonally implicit Runge–Kutta (ESDIRK) methods, and the second-order backward differentiation formula (BDF2) for the high-order flux reconstruction/correction procedure via reconstruction solution of the unsteady Navier–Stokes equations. Pseudo-transient continuation is employed to solve the nonlinear equation at each stage of ESDIRK (excluding the first stage) and each step of BDF2. A Jacobian-free implementation of the restarted generalized minimal residual method solver is employed with a low storage element-Jacobi preconditioner to solve the linear system at each stage of ROW and each pseudo time iteration of ESDIRK and BDF2. Several numerical experiments, including both laminar and turbulent flow simulations, are conducted to carry out the comparison. We observe that the multistage ROW2 and ESDIRK2 are more efficient than the multistep BDF2, and higher-order implicit time integrators are more efficient than lower-order ones. In general, the ESDIRK method allows a larger physical time step size for unsteady flow simulation than the ROW method when the element-Jacobi preconditioner is employed, especially for wall-bounded flows; and the ROW method can be more efficient than the ESDIRK method when the time step size is refined.