Section: New Results
Hybrid time discretizations of high-order
Construction and analysis of a fourth order, energy preserving, explicit time discretization for dissipative linear wave equations.
Participants : Juliette Chabassier, Julien Diaz, Anh-Tuan Ha.
We submitted a paper to M2AN. This paper deals with the construction of a fourth order, energy preserving, explicit time discretization for dissipative linear wave equations. This scheme is obtained by replacing the inversion of a matrix, that comes naturally after using the technique of the Modified Equation on the second order Leap Frog scheme applied to dissipative linear wave equations, by an explicit approximation of its inverse. The stability of the scheme is studied first using an energy analysis, then an eigenvalue analysis. Numerical results in 1D illustrate the good behavior regarding space/time convergence and the efficiency of the newly derived scheme compared to more classical time discretizations. A loss of accuracy is observed for non smooth profiles of dissipation, and we propose an extension of the method that fixes this issue. Finally, we assess the good performance of the scheme for a realistic dissipation phenomenon in Lorentz's materials. This work has been done in collaboration with Sébastien Imperiale (Inria Project-Team M3DISIM) and Alain Anh-Tuan Ha (Internship at Magique 3D in 2016).
Higher-order optimized explicit Runge-Kutta schemes for linear ODEs
Participants : Hélène Barucq, Marc Duruflé, Mamadou N'Diaye.
In this work, we have constructed optimized explicit Runge-Kutta schemes for linear ODEs that we called Linear-ERK. Theses schemes can be applied to the following ODE
where
This work is a chapter of the thesis defended by Mamadou N'diaye on December 8, 2017, under the joint supervision of Hélène Barucq and Marc Duruflé.
High-order locally implicit time schemes for linear ODEs
Participants : Hélène Barucq, Marc Duruflé, Mamadou N'Diaye.
In this work we have proposed a method that combines optimized explicit schemes and implicit schemes to form locally implicit schemes for linear ODEs, including in particular ODEs coming from the space discretization of wave propagation phenomena. This method can be applied to the following ODE
Like in the local time-stepping developed by Grote and co-workers, the computational domain is split into a fine region and a coarse region. The matrix
where
This work has been presented at the Mathias annual Total seminar and is a chapter of the thesis defended by Mamadou N'diaye on December 8, 2017, under the joint supervision of Hélène Barucq and Marc Duruflé.