DeepAI
Log In Sign Up

Efficient implementation of partitioned stiff exponential Runge-Kutta methods

12/02/2019
by   Mahesh Narayanamurthi, et al.
0

Multiphysics systems are driven by multiple processes acting simultaneously, and their simulation leads to partitioned systems of differential equations. This paper studies the solution of partitioned systems of differential equations using exponential Runge-Kutta methods. We propose specific multiphysics implementations of exponential Runge-Kutta methods satisfying stiff order conditions that were developed in [Hochbruck et al., SISC, 1998] and [Luan and Osterman, JCAM, 2014]. We reformulate stiffly–accurate exponential Runge–Kutta methods in a way that naturally allows of the structure of multiphysics systems, and discuss their application to both component and additively partitioned systems. The resulting partitioned exponential methods only compute matrix functions of the Jacobians of individual components, rather than the Jacobian of the full, coupled system. We derive modified formulations of particular methods of order two, three and four, and apply them to solve a partitioned reaction-diffusion problem. The proposed methods retain full order for several partitionings of the discretized problem, including by components and by physical processes.

READ FULL TEXT

page 1

page 2

page 3

page 4

08/26/2019

Partitioned Exponential Methods for Coupled Multiphysics Systems

Multiphysics problems involving two or more coupled physical phenomena a...
09/03/2019

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned ...
03/06/2021

Symplectic GARK methods for Hamiltonian systems

Generalized Additive Runge-Kutta schemes have shown to be a suitable too...
08/19/2021

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

We propose entropy-preserving and entropy-stable partitioned Runge-Kutta...
10/10/2019

Lanczos-like algorithm for the time-ordered exponential: The ∗-inverse problem

The time-ordered exponential of a time-dependent matrix A(t) is defined ...