Partitioned exponential methods for coupled multiphysics systems

From MaRDI portal
Publication:2227738



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.


New partitioned exponential integrators are developed for coupled multiphysics problems by changing the information between components of the coupling terms. Two different approaches are considered for the construction and analysis of these methods, one based on splitting the component functions into linear and nonlinear terms (split-RHS methods), and the other is based on approximating the Jacobians of individual components (W-type methods). Two new formulations of partitioned exponential methods are proposed: partitioned methods of Rosenbrock-exponential type (PEXPW), and partitioned Runge-Kutta style exponential methods (PEPIRKW). The matrix-exponential-like functions of the proposed partitioned methods are evaluated on the Jacobians of the individual component functions, whereas these functions are applied to the full (coupled) Jacobian in an unpartitioned method. If the individual Jacobians have computationally favorable structures, then the computational expense of evaluating matrix-exponential-like functions is greatly reduced. In reaction-diffusion systems with two or more species, the diffusion Jacobian is block-diagonal, the matrix-exponential-like functions can be evaluated on individual blocks in parallel. The numerical experiments show that using this strategy the partitioned exponential methods are at least twice as fast as the unpartitioned counterpart. One drawback of the partitioned exponential methods is the reduced stability. The numerical tests illustrate that the partitioned methods can exhibit more stable behavior than an unpartitioned method in some stiffness regimes, wheras in some very stiff regimes, partitioned methods can fail to obtain a solution.



Cites work



Describes a project that uses

Uses Software






This page was built for publication: Partitioned exponential methods for coupled multiphysics systems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2227738)