High-order discretizations for the wave equation based on the modified equation technique - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Communication Dans Un Congrès Année : 2010

High-order discretizations for the wave equation based on the modified equation technique

Résumé

Highly accurate solution of the full wave equation implies very high computa- tional burdens. Indeed, to improve the accuracy of the numerical solution, one must considerably reduce the space step. Obviously this will result in increasing the number of unknowns of the discrete problem. Besides, the time step is linked to the space step through the CFL (Courant-Friedrichs-Lewy) condition. In the three- dimensional case the problem can have more than ten million unknowns, which must be evaluated at each time-iteration. However, high-order numerical methods can be used for computing accurate solutions with larger space and time steps. Recently, Joly and Gilbert (cf. [2]) have optimized the Modified Equation Technique (MET), which was proposed in [1] for solving the wave equation, and it seems to be very promising given some improvements. In this work, we apply this technique in a new way. Normally, most of the study devoted to the solution of the wave equation consider first the space discretization of the system before addressing the question of the time discretization. We intend here to invert the discretization process by applying first the time discretization using the MET and then to consider the space discretization. The time discretization causes high-order operators to appear (such as p-Laplacian) and we have therefore to consider appropriate methods to discretize them. The Discontinuous Galerkin Methods are well adapted to this discretization, since they allow to consider piecewise discontinuous functions. In particular, using the Interior Penalty Discontinuous Galerkin (IPDG) method (cf. [3, 4]), one can en- force through the elements high-order transmission conditions, which are adapted to the high order operators to be discretized. References [1] Shubin, G. R. and Bell, J. B. A modified equation approach to constructing fourth- order methods for acoustic wave propagation.. SIAM J. Sci. Statist. Comput., vol. 8, pp.135–151, 1987. [2] Gilbert, J.-C. and Joly, P. Higher order time stepping for second order hyperbolic problems and optimal cfl conditions. Numerical Analysis and Scientific Computing for PDE's and their Challenging Applications, 2006. [3] Grote, M. J. and Schneebeli, A. and Sch ̈tzau, D. Dicontinuous Galerkin finite o element method for the wave equation. SIAM J. on Numerical Analysis, vol.44, pp.2408–2431, 2006. [4] Mozolevski, I. and Suli, E. hp-version interior penalty DGFEMs for the biharmonic equation. Oxford University Computing Laboratory, 2004.
Fichier non déposé

Dates et versions

inria-00527136 , version 1 (18-10-2010)

Identifiants

  • HAL Id : inria-00527136 , version 1

Citer

Cyril Agut, Julien Diaz. High-order discretizations for the wave equation based on the modified equation technique. XIV Spanish-French Jacques-Louis Lions School, Sep 2010, A Coruna, Spain. ⟨inria-00527136⟩
224 Consultations
0 Téléchargements

Partager

Gmail Facebook X LinkedIn More