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.

Convergence properties of the Runge-Kutta-Chebyshev method

Convergence of Newton'S Method for Systems of Equations with Constant Rank Derivatives

The Stability and Convergence of Computing Long-Time Behaviour

Convergence of a fully discrete finite difference scheme for the Korteweg-de Vries equation

Explicit stabilized multirate method for stiff differential equations

Stochastic C-stability and B-consistency of explicit and implicit Euler-type schemes

A Unifying Theory for Runge--Kutta-like Time Integrators: Convergence and Stability

ANALYSIS OF THE IMPLICIT-EXPLICIT ULTRA-WEAK DISCONTINUOUS GALERKIN METHOD FOR CONVECTION-DIFFUSION PROBLEMS

Necessary and sufficient conditions for strong stability of explicit Runge-Kutta methods

Explicit Runge Kutta Methods that Alleviate Order Reduction

Stability and Convergence analysis of a Crank-Nicolson Galerkin scheme for the fractional Korteweg-de Vries equation

The Runge--Kutta discontinuous Galerkin method with stage-dependent polynomial spaces for hyperbolic conservation laws

Convergence and Superconvergence Analysis of a Nonconforming Finite Element Variable-Time-Step BDF2 Implicit Scheme for Linear Reaction-Diffusion Equations

Enforcing strong stability of explicit Runge--Kutta methods with superviscosity

Performance analysis of relaxation Runge–Kutta methods

Electrical polarization of nuclear spins in a breakdown regime of quantum Hall effect

Many-Stage Optimal Stabilized Runge–Kutta Methods for Hyperbolic Partial Differential Equations

Strong Stability Preserving Integrating Factor Runge-Kutta Methods

Conservative iterative methods for implicit discretizations of conservation laws

High Resolution Schemes for Conservation Laws and Convection-Diffusion Equations with Varying Time and Space Grids

Superconvergence of Discontinuous Galerkin and Local Discontinuous Galerkin Schemes for Linear Hyperbolic and Convection-Diffusion Equations in One Space Dimension