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.

Parallel Iteration Algorithms Research on the Continuous Extension for Stiff Systems

On the Convergence of the Accelerated Riccati Iteration Method.

High-Order Schemes of Exponential Time Differencing for Stiff Systems with Nondiagonal Linear Part

On convergence of waveform relaxation for nonlinear systems of ordinary differential equations

High order implicit time integration schemes on multiresolution adaptive grids for stiff PDEs

The Generalization of Diagonally Implicit Runge–Kutta–Nyström Method with Controllable Numerical Dissipation for Structural Dynamics

Generalized high-order iterative methods for solutions of nonlinear systems and their applications

The Existence the Solution of Nonlinear Discrete Schemes and Convergence of a Linearized Iterative Method for time-dependent PNP Equations

Convergence properties of the Runge-Kutta-Chebyshev method

Global convergence of iterative solvers for problems of nonlinear magnetostatics

The Convergence and Fractals of Attractive Basins

Projective Integration Methods in the Runge-Kutta Framework and the Extension to Adaptivity in Time

A note on continuous-stage Runge-Kutta methods

Efficient and Parallel Solution of High-order Continuous Time Galerkin for Dissipative and Wave Propagation Problems

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

Algebraically Stable SDIRK Methods with Controllable Numerical Dissipation for First/second-Order Time-Dependent Problems

New high-order convergence iteration methods without employing derivatives for solving nonlinear equations

Construction of symplectic (partitioned) Runge-Kutta methods with continuous stage

A Successive over Relaxation Implicit Iterative Algorithm for Solving Stochastic Linear Systems with Markov Jumps

Performance analysis of relaxation Runge–Kutta methods