Iterative operator-splitting methods with higher-order time integration methods and applications for parabolic partial differential equations (Q929923)

Computational simulations of bioremediation, heat transfer, transport and chemical reaction problems require accurate numerical methods. One of the widely spread methods is to transform time and space dependent problems described by one or a system of partial differential equations to a system of ordinary differential equations (ODEs) using a certain space discretization method like the finite element method and then solve the system of ODEs. The presented method is a combination of two efficient methods: an iterative and an operator-splitting method by using the advantages of both of them. Splitting methods have the benefit in decoupling multiphysics problems into simpler physical problems and solve this reduced problems in each time step. The iterative method is used in each time step for solving reduced problems iteratively. An analysis of the convergence and the rate of convergence of this iterative operator-splitting method is presented. To obtain high-order methods Runge-Kutta and backward differentiation formula methods are used. For both commutative and noncommutative operator of the continuous equations the A-stability of the method is proved under certain conditions. Various numerical experiments for systems of ODEs or parabolic equations also with nonlinearities are presented. Comparisons of various methods of higher order together with various splitting methods are shown.