Abstract:Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. The traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

EXPODE - Advanced Exponential Time Integration Toolbox for MATLAB, Code Documentation

FESTUNG: A MATLAB /GNU Octave toolbox for the discontinuous Galerkin method. Part IV: Generic problem framework and model-coupling interface

A toolbox of Equation-Free functions in Matlab\Octave for efficient system level simulation

EPIRK-W and EPIRK-K time discretization methods

ODE Test Problems: a MATLAB suite of initial value problems

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

A new approach to constructing efficient stiffly accurate exponential propagation iterative methods of Runge-Kutta type (EPIRK)

Fourier Transform of the Stretched Exponential Function: Analytic Error Bounds, Double Exponential Transform, and Open-Source Implementation libkww

KSSOLV—a MATLAB Toolbox for Solving the Kohn-Sham Equations

Towards a MATLAB Toolbox to compute backstepping kernels using the power series method

A fourth-order exponential time differencing scheme with dimensional splitting for non-linear reaction-diffusion systems

Partitioned Exponential Methods for Coupled Multiphysics Systems

A Matlab Toolbox for Extended Dynamic Mode Decomposition Based on Orthogonal Polynomials and p-q Quasi-Norm Order Reduction

A new algorithm for automatic computation of solitary wave solutions to nonlinear partial differential equations based on the Exp-function method.

Exponential time propagators for elastodynamics

Exasim: Generating Discontinuous Galerkin Codes for Numerical Solutions of Partial Differential Equations on Graphics Processors

Exponential time differencing for matrix-valued dynamical systems

Simple MATLAB and Python scripts for multi‐exponential analysis

A Matlab toolbox for fractional relaxation-oscillation equations

Stochastic Runge-Kutta Software Package for Stochastic Differential Equations