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

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.

[1]  I. Epstein,et al.  An Introduction to Nonlinear Chemical Dynamics , 1998 .

[2]  N. Babaeva,et al.  Two-dimensional modelling of positive streamer dynamics in non-uniform electric fields in air , 1996 .

[3]  Stéphane Descombes,et al.  On the local and global errors of splitting approximations of reaction–diffusion equations with high spatial gradients , 2007, Int. J. Comput. Math..

[4]  D. Veynante,et al.  Stabilization of a Turbulent Premixed Flame Using a Nanosecond Repetitively Pulsed Plasma , 2006, IEEE Transactions on Plasma Science.

[5]  Andrew G. Glen,et al.  APPL , 2001 .

[6]  Marc Massot,et al.  Operator splitting for nonlinear reaction-diffusion systems with an entropic structure : singular perturbation and order reduction , 2004, Numerische Mathematik.

[7]  O. Koch,et al.  Embedded Split-Step Formulae for the Time Integration of Nonlinear Evolution Equations , 2010 .

[8]  I. A. Kossyi,et al.  Kinetic scheme of the non-equilibrium discharge in nitrogen-oxygen mixtures , 1992 .

[9]  Habib N. Najm,et al.  Regular Article: A Semi-implicit Numerical Scheme for Reacting Flow , 1999 .

[10]  Karline Soetaert,et al.  Solving Ordinary Differential Equations in R , 2012 .

[11]  Mechthild Thalhammer,et al.  The Lie–Trotter splitting method for nonlinear evolutionary problems involving critical parameters. An exact local error representation and application to nonlinear Schrödinger equations in the semi-classical regime , 2010 .

[12]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[13]  Mechthild Thalhammer,et al.  The Lie–Trotter splitting for nonlinear evolutionary problems with critical parameters: a compact local error representation and application to nonlinear Schrödinger equations in the semiclassical regime , 2013 .

[14]  Assyr Abdulle,et al.  Fourth Order Chebyshev Methods with Recurrence Relation , 2001, SIAM J. Sci. Comput..

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

[16]  Chemical Oscillations and Instabilities. Non-linear Chemical Kinetics. : Clarendon Press, Oxford, 1990, xvi + 453 pp., £50.00 , 1992 .

[17]  Marc Massot,et al.  New Resolution Strategy for Multi-scale Reaction Waves using Time Operator Splitting and Space Adaptive Multiresolution: Application to Human Ischemic Stroke , 2011 .

[18]  Marc Massot,et al.  Simulation of human ischemic stroke in realistic 3D geometry , 2010, Commun. Nonlinear Sci. Numer. Simul..

[19]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[20]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[21]  Marc Massot,et al.  New Resolution Strategy for Multiscale Reaction Waves using Time Operator Splitting, Space Adaptive Multiresolution, and Dedicated High Order Implicit/Explicit Time Integrators , 2012, SIAM J. Sci. Comput..

[22]  Habib N. Najm,et al.  Modeling unsteady reacting flow with operator splitting and ISAT , 2006 .

[23]  J. Lowke,et al.  Streamer propagation in air , 1997 .

[24]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[25]  Gilbert Strang,et al.  Accurate partial difference methods , 1964 .

[26]  Marc Massot,et al.  New Resolution Strategies for Multi-scale Reaction Waves: Optimal Time Operator Splitting and Space Adaptive Multiresolution , 2011, CLEI Electron. J..

[27]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[28]  P. S. Wyckoff,et al.  A Semi-implicit Numerical Scheme for Reacting Flow , 1998 .

[29]  G. Strang.  Accurate partial difference methods I: Linear cauchy problems , 1963 .

[30]  A. Kulikovsky.  Positive streamer between parallel plate electrodes in atmospheric pressure air , 1997 .