High-order Padé and Singly Diagonally Runge-Kutta schemes for linear ODEs, application to wave propagation problems
BARUCQ, Hélène
Advanced 3D Numerical Modeling in Geophysics [Magique 3D]
Université de Pau et des Pays de l'Adour [UPPA]
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]
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]
Advanced 3D Numerical Modeling in Geophysics [Magique 3D]
Université de Pau et des Pays de l'Adour [UPPA]
BARUCQ, Hélène
Advanced 3D Numerical Modeling in Geophysics [Magique 3D]
Université de Pau et des Pays de l'Adour [UPPA]
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]
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]
< Reduce
Advanced 3D Numerical Modeling in Geophysics [Magique 3D]
Université de Pau et des Pays de l'Adour [UPPA]
Language
en
Article de revue
This item was published in
Numerical Methods for Partial Differential Equations. 2018, vol. 34, p. 760-798
Wiley
English Abstract
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 ...Read more >
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.Read less <
Origin
Hal imported