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

TL;DR

This paper presents an adaptive time splitting method for stiff evolutionary PDEs that ensures effective error control across multiple physical time scales, improving simulation accuracy for complex multi-scale phenomena.

Contribution

It introduces a novel adaptive time splitting scheme combining second order and embedded lower order methods with enhanced error control for multi-scale PDEs.

Findings

Effective error control across diverse time scales

Successful simulation of nonlinear chemical dynamics

Accurate modeling of nanosecond pulsed gas discharges

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.