English

Beyond Strang: A practical assessment of some second-order 3-splitting methods

Numerical Analysis 2023-08-14 v2 Numerical Analysis Analysis of PDEs Computational Physics Fluid Dynamics Plasma Physics

Abstract

Operator splitting is a popular divide-and-conquer strategy for solving differential equations. Typically, the right-hand side of the differential equation is split into a number of parts that are then integrated separately. Many methods are known that split the right-hand side into two parts. This approach is limiting, however, and there are situations when 3-splitting is more natural and ultimately more advantageous. The second-order Strang operator-splitting method readily generalizes to a right-hand side splitting into any number of operators. It is arguably the most popular method for 3-splitting because of its efficiency, ease of implementation, and intuitive nature. Other 3-splitting methods exist, but they are less well-known, and \rev{analysis and} evaluation of their performance in practice are scarce. We demonstrate the effectiveness of some alternative 3-split, second-order methods to Strang splitting on two problems: the reaction-diffusion Brusselator, which can be split into three parts that each have closed-form solutions, and the kinetic Vlasov--Poisson equations that is used in semi-Lagrangian plasma simulations. We find alternative second-order 3-operator-splitting methods that realize efficiency gains of 10\%--20\% over traditional Strang splitting. Our analysis for the practical assessment of efficiency of operator-splitting methods includes the computational cost of the integrators and can be used in method design.

Keywords

Cite

@article{arxiv.2302.08034,
  title  = {Beyond Strang: A practical assessment of some second-order 3-splitting methods},
  author = {Raymond J. Spiteri and Arash Tavassoli and Siqi Wei and Andrei Smolyakov},
  journal= {arXiv preprint arXiv:2302.08034},
  year   = {2023}
}
R2 v1 2026-06-28T08:41:23.500Z