Hamiltonian discontinuous Galerkin FEM for linear, stratified (in)compressible Euler equations: internal gravity waves

JOURNAL OF COMPUTATIONAL PHYSICS

卷:330

Article

van Oers, Alexander M.^1,2  Maas, Leo R. M.^1,2,3  Bokhove, Onno⁴ 

[1] NIOZ Royal Netherlands Inst Sea Res, POB 59, NL-1790 AB Texel, Netherlands

[2] Univ Utrecht, POB 59, NL-1790 AB Texel, Netherlands

[3] Univ Utrecht, Inst Marine & Atmospher Res Utrecht, Imau, Princetonpl 5, NL-3584 CC Utrecht, Netherlands

[4] Univ Leeds, Sch Math, Leeds LS2 9JT, W Yorkshire, England

关键词: Linear stratified Euler equations;    Hamiltonian structure;    Discontinuous Galerkin method;    Internal gravity waves;    Wave attractors;   

DOI  :  10.1016/j.jcp.2016.10.032

来源: Elsevier

【 摘 要 】

The linear equations governing internal gravity waves in a stratified ideal fluid possess a Hamiltonian structure. A discontinuous Galerkin finite element method has been developed in which this Hamiltonian structure is discretized, resulting in conservation of discrete analogs of phase space and energy. This required (i) the discretization of the Hamiltonian structure using alternating flux functions and symplectic time integration, (ii) the discretization of a divergence-free velocity field using Dirac's theory of constraints and (iii) the handling of large-scale computational demands due to the 3-dimensional nature of internal gravity waves and, in confined, symmetry-breaking fluid domains, possibly its narrow zones of attraction. (C) 2016 The Authors. Published by Elsevier Inc.

【 授权许可】

Free