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.

Powell’s method-based nonlinear least-squares data fitting for the Mittag-Leffler relaxation function

Parameter Identification in Mixed Brownian–fractional Brownian Motions Using Powell's Optimization Algorithm

Modeling Enthalpy Relaxation Using the Mittag-Leffler Function

Computation of the inverse Mittag–Leffler function and its application to modeling ultraslow dynamics

A STOCHASTIC EQUIVALENT LINEARIZATION FOR ANALYSIS OF MULTIPLE-DEGREE-OF-FREEDOM FLEXIBLE STRUCTURES

Sampling Points-Independent Identification of the Fractional Maxwell Model of Viscoelastic Materials Based on Stress Relaxation Experiment Data

A fast numerical scheme for fractional viscoelastic models of wave propagation

Application of fractional calculus methods to viscoelastic behaviours of solid propellants

Function Based Nonlinear Least Squares and Application to Jelinski--Moranda Software Reliability Model

A mathematical model for fitting and predicting relaxation modulus and simulating viscoelastic responses

A meshless method to solve nonlinear variable-order time fractional 2D reaction–diffusion equation involving Mittag-Leffler kernel

Direct integration of measured viscoelastic relaxation data in time-domain finite element simulations

Optimal a priori error estimate of relaxation-type linear finite element method for nonlinear Schrödinger equation

Damped least squares method for nonlinear mixed additive and multiplicative errors model

An extended search method for identifying optimal parameters of the generalized Maxwell model

Ramp-hold relaxation solutions for the KVFD model applied to soft viscoelastic media

A Coupled Numerical Method for Relaxed Minimizers and Inhomogeneous Microstructures

Sensitivity analysis and determination of free relaxation parameters for the weakly-compressible MRT-LBM schemes

Performance analysis of relaxation Runge–Kutta methods

Variational Multiscale Interpolating Element-Free Galerkin Method for the Nonlinear Darcy–Forchheimer Model

Least-Squares finite element method for the simulation of sea-ice motion