Adaptive time splitting method for multi-scale evolutionary partial differential equations

From MaRDI portal



Abstract: This paper introduces an adaptive time splitting technique for the solution of stiff evolutionary PDEs that guarantees an effective error control of the simulation, independent of the fastest physical time scale for highly unsteady problems. The strategy considers a second order Strang method and another lower order embedded splitting scheme that takes into account potential loss of order due to the stiffness featured by time-space multi-scale phenomena. The scheme is then built upon a precise numerical analysis of the method and a complementary numerical procedure, conceived to overcome classical restrictions of adaptive time stepping schemes based on lower order embedded methods, whenever asymptotic estimates fail to predict the dynamics of the problem. The performance of the method in terms of control of integration errors is evaluated by numerical simulations of stiff propagating waves coming from nonlinear chemical dynamics models as well as highly multi-scale nanosecond repetitively pulsed gas discharges, which allow to illustrate the method capabilities to consistently describe a broad spectrum of time scales and different physical scenarios for consecutive discharge/post-discharge phases.


The authors consider the question of how to estimate the splitting error when solving nonlinear reaction-diffusion equations by applying the Strang (-Marchuk) splitting to partial differential equations. For the discretization in space they refer to an adaptive multiresolution technique. They prove theorems about the error of the original and a shifted (first order) splitting for a scalar reaction-diffusion equation, and, based on this, develop a strategy to compute the actual time step as well as an interval in which the estimates are valid, for more general equations and systems of equations, and for a given accuracy. Their numerical tests include the stiff Belousov-Zhabotinski model and multi-scale gas discharges, applying high-order methods in time in order to make the main error source the splitting method.











This page was built for publication: Adaptive time splitting method for multi-scale evolutionary partial differential equations

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