JOURNAL OF COMPUTATIONAL PHYSICS | 卷:350 |
Symplectic Hamiltonian HDG methods for wave propagation phenomena | |
Article | |
Sanchez, M. A.1,2,3  Ciuca, C.4  Nguyen, N. C.5  Peraire, J.5  Cockburn, B.1  | |
[1] Univ Minnesota, Sch Math, Minneapolis, MN 55455 USA | |
[2] Pontificia Univ Catolica Chile, Inst Math & Computat Engn, Sch Engn, Santiago, Chile | |
[3] Pontificia Univ Catolica Chile, Fac Math, Santiago, Chile | |
[4] Imperial Coll London, Dept Aeronaut, London SW7 2AZ, England | |
[5] MIT, Dept Aeronaut & Astronaut, Cambridge, MA 02139 USA | |
关键词: Finite element methods; Discontinuous Galerkin methods; Hybrid/mixed methods; Acoustic wave equation; Hamiltonian systems; Symplectic time integrators; Energy conservation; | |
DOI : 10.1016/j.jcp.2017.09.010 | |
来源: Elsevier | |
【 摘 要 】
We devise the first symplectic Hamiltonian hybridizable discontinuous Galerkin (HDG) methods for the acoustic wave equation. We discretize in space by using a Hamiltonian HDG scheme, that is, an HDG method which preserves the Hamiltonian structure of the wave equation, and in time by using symplectic, diagonally implicit and explicit partitioned Runge-Kutta methods. The fundamental feature of the resulting scheme is that the conservation of a discrete energy, which is nothing but a discrete version of the original Hamiltonian, is guaranteed. We present numerical experiments which indicate that the method achieves optimal approximations of order k + 1 in the L-2-norm when polynomials of degree k >= 0 and Runge-Kutta time-marching methods of order k + 1 are used. In addition, by means of post-processing techniques and by increasing the order of the Runge-Kutta method to k+2, we obtain superconvergent approximations of order k+2 in the L-2-norm for the displacement and the velocity. We also present numerical examples that corroborate that the methods conserve energy and that they compare favorably with dissipative HDG schemes, of similar accuracy properties, for long-time simulations. (C) 2017 Elsevier Inc. All rights reserved.
【 授权许可】
Free
【 预 览 】
Files | Size | Format | View |
---|---|---|---|
10_1016_j_jcp_2017_09_010.pdf | 1361KB | download |