DeepAI AI Chat
Log In Sign Up

Design of DIRK Schemes with High Weak Stage Order

by   Abhijit Biswas, et al.

Runge-Kutta (RK) methods may exhibit order reduction when applied to certain stiff problems. While fully implicit RK schemes exist that avoid order reduction via high-stage order, DIRK (diagonally implicit Runge-Kutta) schemes are practically important due to their structural simplicity; however, these cannot possess high stage order. The concept of weak stage order (WSO) can also overcome order reduction, and it is compatible with the DIRK structure. DIRK schemes of WSO up to 3 have been proposed in the past, however, based on a simplified framework that cannot be extended beyond WSO 3. In this work a general theory of WSO is employed to overcome the prior WSO barrier and to construct practically useful high-order DIRK schemes with WSO 4 and above. The resulting DIRK schemes are stiffly accurate, L-stable, have optimized error coefficients, and are demonstrated to perform well on a portfolio of relevant ODE and PDE test problems.


Algebraic Structure of the Weak Stage Order Conditions for Runge-Kutta Methods

Runge-Kutta (RK) methods may exhibit order reduction when applied to sti...

Very High-Order A-stable Stiffly Accurate Diagonally Implicit Runge-Kutta Methods

A numerical search approach is used to design high-order diagonally impl...

High order, semi-implicit, energy stable schemes for gradient flows

We introduce a class of high order accurate, semi-implicit Runge-Kutta s...

Eliminating Order Reduction on Linear, Time-Dependent ODEs with GARK Methods

When applied to stiff, linear differential equations with time-dependent...

Fast and Stable Schemes for Phase Fields Models

We propose and analyse new stabilized time marching schemes for Phase Fi...