Abstract:The International Journal of High Performance Computing Applications, Ahead of Print. Recently, global and local relaxation Runge–Kutta methods have been developed for guaranteeing the conservation, dissipation, or other solution properties for general convex functionals whose dynamics are crucial for an ordinary differential equation solution. These novel time integration procedures have an application in a wide range of problems that require dynamics-consistent and stable numerical methods. The application of a relaxation scheme involves solving scalar nonlinear algebraic equations to find the relaxation parameter. Even though root-finding may seem to be a problem technically straightforward and computationally insignificant, we address the problem at scale as we solve full-scale industrial problems on a CPU-powered supercomputer and show its cost to be considerable. In particular, we apply the relaxation schemes in the context of the compressible Navier–Stokes equations and use them to enforce the correct entropy evolution. We use seven different algorithms to solve for the global and local relaxation parameters and analyze their strong scalability. As a result of this analysis, within the global relaxation scheme, we recommend using Brent's method for problems with a low polynomial degree and of small sizes for the global relaxation scheme, while secant proves to be the best choice for higher polynomial degree solutions and large problem sizes. For the local relaxation scheme, we recommend secant. Further, we compare the schemes' performance using their most efficient implementations, where we look at their effect on the timestep size, overhead, and weak scalability. We show the global relaxation scheme to be always more expensive than the local approach—typically 1.1–1.5 times the cost. At the same time, we highlight scenarios where the global relaxation scheme might underperform due to its increased communication requirements. Finally, we present an analysis that sets expectations on the computational overhead anticipated based on the system properties.

Uniform accuracy of implicit-explicit Runge-Kutta methods for linear hyperbolic relaxation systems

Uniform accuracy of implicit-explicit Runge-Kutta (IMEX-RK) schemes for hyperbolic systems with relaxation

Uniform accuracy of implicit-explicit backward differentiation formulas (IMEX-BDF) for linear hyperbolic relaxation systems

Implicit-Explicit Runge-Kutta schemes for hyperbolic systems with stiff relaxation and applications

Uniformly Accurate Schemes for Hyperbolic Systems with Relaxation

Implicit-explicit Relaxation Runge-Kutta Methods: Construction, Analysis and Applications to PDEs

Application of a Fourth-Order Relaxation Scheme to Hyperbolic Systems of Conservation Laws

Entropy-Preserving and Entropy-Stable Relaxation IMEX and Multirate Time-Stepping Methods

Energy diminishing implicit-explicit Runge--Kutta methods for gradient flows

Energy diminishing implicit-explicit Runge–Kutta methods for gradient flows

Asymptotic Analysis of IMEX-RK Methods for ES-BGK Model at Navier-Stokes level

A Relaxation Scheme for a Multi-Class Lighthill-Whitham-Richards Traffic Flow Model

Performance analysis of relaxation Runge–Kutta methods

Implicit-explicit time discretization schemes for a class of semilinear wave equations with nonautonomous dampings

Energy Conservative Relaxation-Free Runge-Kutta Schemes

A Second-Order, Global-in-Time Energy Stable Implicit-Explicit Runge–Kutta Scheme for the Phase Field Crystal Equation

On the Stability of IMEX Upwind gSBP Schemes for 1D Linear Advection-Diffusion Equations

Error analysis of an implicit-explicit time discretization scheme for semilinear wave equations with application to multiscale problems

High order strong stability preserving multi-derivative implicit and IMEX Runge--Kutta methods with asymptotic preserving properties

IMEX-RK methods for Landau-Lifshitz equation with arbitrary damping