BARUCQ, Hélène; DURUFLÉ, Marc; N'DIAYE, Mamadou

doi:10.1002/num.22228

Métadonnées

Licence d’utilisation du document

BARUCQ, Hélène
Advanced 3D Numerical Modeling in Geophysics [Magique 3D]
Université de Pau et des Pays de l'Adour [UPPA]

DURUFLÉ, Marc
Advanced 3D Numerical Modeling in Geophysics [Magique 3D]
Institut de Mathématiques de Bordeaux [IMB]

N'DIAYE, Mamadou
Advanced 3D Numerical Modeling in Geophysics [Magique 3D]
Université de Pau et des Pays de l'Adour [UPPA]

Langue

Article de revue

Ce document a été publié dans

Numerical Methods for Partial Differential Equations. 2018, vol. 34, p. 760-798

Wiley

Résumé en anglais

In this paper we address the problem of constructing high-order implicit time schemes for wave equations. We consider two classes of one-step A-stable schemes adapted to linear Ordinary Differential Equation (ODE). The first class, which is not dissipative is based upon the diagonal Padé approximant of exponential function. In this class, the obtained schemes have the same stability function as Gauss Runge-Kutta (Gauss RK) schemes. They have the advantage to involve solution of smaller linear system at each time step compared to Gauss RK. The second class of schemes are constructed such that they require the inversion of a unique linear system several times at each time step like the Singly Diagonally Runge-Kutta (SDIRK) schemes. While the first class of schemes is constructed for an arbitrary order of accuracy, the second class schemes is given up to order 12. The performance assessment we provide shows a very good level of accuracy with very reduced computational costs for both class of schemes. But diagonal Padé schemes seem to be more accurate and more robust.< Réduire

Métadonnées

Partager cette publication !

Licence d’utilisation du document

High-order Padé and Singly Diagonally Runge-Kutta schemes for linear ODEs, application to wave propagation problems

Langue

Ce document a été publié dans

Résumé en anglais

URI

DOI

Origine

Unités de recherche