Skip to Main content Skip to Navigation
Journal articles

High-order locally A-stable implicit schemes for linear ODEs

Hélène Barucq 1 Marc Duruflé 1 Mamadou N'Diaye 2, 1
1 Magique 3D - Advanced 3D Numerical Modeling in Geophysics
LMAP - Laboratoire de Mathématiques et de leurs Applications [Pau], Inria Bordeaux - Sud-Ouest
Abstract : Accurate simulations of wave propagation in complex media like Earth subsurface can be performed with a reasonable computational burden by using hybrid meshes stuffing fine and coarse cells. Locally implicit time discretiza-tions are then of great interest. They indeed allow using unconditionally stable schemes in the regions of computational domain covered by small cells. The receivable values of the time step are then increased which reduces the computational costs while limiting the dispersion effects. In this work we construct a method that combines optimized explicit schemes and implicit schemes to form locally implicit schemes for linear ODEs, including in particular semi-discretized wave problems that are considered herein for numerical experiments. Both the explicit and implicit schemes used are one-step methods constructed using their stability function. The stability function of the explicit schemes are computed by maximizing the time step that can be chosen. The implicit schemes used are unconditionally stable and do not necessary require the same number of stages as the explicit schemes. The performance assessment we provide shows a very good level of accuracy for locally implicit schemes. It also shows that a locally implicit scheme is a good compromise between purely explicit and purely implicit schemes in terms of computational time and memory usage.
Complete list of metadata
Contributor : Marc Duruflé Connect in order to contact the contributor
Submitted on : Friday, February 5, 2021 - 4:19:13 PM
Last modification on : Friday, February 26, 2021 - 11:58:01 AM


Files produced by the author(s)




Hélène Barucq, Marc Duruflé, Mamadou N'Diaye. High-order locally A-stable implicit schemes for linear ODEs. Journal of Scientific Computing, Springer Verlag, 2020, 85, ⟨10.1007/s10915-020-01313-x⟩. ⟨hal-02419543v3⟩



Les métriques sont temporairement indisponibles