Adaptive time splitting method for multi-scale evolutionary partial differential equations
For computational scientists solving multi-scale PDEs, this method offers improved error control over classical adaptive schemes, but the improvement is incremental.
The paper introduces an adaptive time splitting method for stiff evolutionary PDEs that ensures effective error control independent of the fastest time scale, demonstrated on nonlinear chemical dynamics and nanosecond pulsed gas discharges.
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.