Splitting methods with complex times for parabolic equations

Using composition procedures, we build up high order splitting methods to solve evolution equations posed in finite or infinite dimensional spaces. Since high-order splitting methods with real time are known to involve large and/or negative time steps, which destabilizes the overall procedure, the key point of our analysis is, we develop splitting methods that use complex time steps having positive real part: going to the complex plane allows to considerably increase the accuracy, while keeping small time steps; on the other hand, restricting our attention to time steps with positive real part makes our methods more stable, and in particular well adapted in the case when the considered evolution equation involves unbounded operators in infinite dimensional spaces, like parabolic (diffusion) equations.We provide a thorough analysis in the case of linear equations posed in general Banach spaces. We also numerically investigate the nonlinear situation. We illustrate our results in the case of (linear and nonlinear) parabolic equations.

[1]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[2]  Stéphane Descombes,et al.  Convergence of a splitting method of high order for reaction-diffusion systems , 2001, Math. Comput..

[3]  Robert I. McLachlan,et al.  Composition methods in the presence of small parameters , 1995 .

[4]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[5]  S. Fauve,et al.  Localized structures generated by subcritical instabilities , 1988 .

[6]  M. Suzuki,et al.  Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simulations , 1990 .

[7]  Michelle Schatzman,et al.  Numerical Integration of Reaction–Diffusion Systems , 2002, Numerical Algorithms.

[8]  Etienne Forest.  Canonical integrators as tracking codes (or how to integrate perturbation theory with tracking) , 1989 .

[9]  G. Sussman,et al.  Chaotic Evolution of the Solar System , 1992, Science.

[10]  Stefan Vandewalle,et al.  Multigrid Waveform Relaxation for Anisotropic Partial Differential Equations , 2002, Numerical Algorithms.

[11]  Alexander Ostermann,et al.  Exponential splitting for unbounded operators , 2009, Math. Comput..

[12]  Q. Sheng.  Solving Linear Partial Differential Equations by Exponential Splitting , 1989 .

[13]  Stéphane Descombes,et al.  Convergence of the Peaceman-Rachford approximation for reaction-diffusion systems , 2003, Numerische Mathematik.

[14]  John E. Chambers,et al.  Symplectic Integrators with Complex Time Steps , 2003 .

[15]  André D. Bandrauk,et al.  Complex integration steps in decomposition of quantum exponential evolution operators , 2006 .

[16]  Fernando Casas,et al.  On the necessity of negative coefficients for operator splitting schemes of order higher than two , 2005 .

[17]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[18]  Boun Oumar Dia.  Méthodes de directions alternées d'ordre élevé en temps , 1996 .

[19]  Zurab Gegechkori,et al.  The fourth order accuracy decomposition scheme for an evolution problem , 2004 .

[20]  Creutz,et al.  Higher-order hybrid Monte Carlo algorithms. , 1989, Physical review letters.

[21]  A. Ostermann,et al.  High order splitting methods for analytic semigroups exist , 2009 .

[22]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[23]  Zurab Gegechkori,et al.  High degree precision decomposition method for the evolution problem with an operator under a split form , 2002 .

[24]  M. Suzuki,et al.  General theory of fractal path integrals with applications to many‐body theories and statistical physics , 1991 .

[25]  H. H. Rosenbrock,et al.  Some general implicit processes for the numerical solution of differential equations , 1963, Comput. J..

[26]  Tasso J. Kaper,et al.  N th-order operator splitting schemes and nonreversible systems , 1996 .