EEMCS

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader 1193 Kb

Official URL: http://www.math.utwente.nl/publications

Abstract

A discontinuous Galerkin finite element method (DGFEM) has been developed and tested for linear, three-dimensional, rotating incompressible Euler equations. These equations admit complicated wave solutions.
The numerical challenges concern: (i) discretisation of a divergence-free velocity field; (ii) discretisation of geostrophic boundary conditions combined with no-normal flow at solid walls; (iii) discretisation of the conserved, Hamiltonian dynamics of the inertial-waves; and, (iv) large-scale computational demands owing to the three-dimensional nature of inertial-wave dynamics and possibly its narrow zones of chaotic attraction. These issues have been resolved: (i) by employing Dirac’s method of constrained Hamiltonian dynamics to our DGFEM for linear, compressible flows, thus enforcing the incompressibility constraints; (ii) by enforcing no-normal flow at solid walls in a weak form and geostrophic tangential flow —along the wall; (iii) by applying a symplectic time discretisation; and, (iv) by combining PETSc’s linear algebra routines with our high-level software.
We compared our simulations with exact solutions of three-dimensional compressible and incompressible flows, in (non)rotating periodic and partly periodic cuboids (Poincar´e waves). Additional verifications concerned semi-analytical eigenmode solutions in rotating cuboids with solid walls.

Export this item as:

• BibTeX
• XML

To correct this item please ask your editor

Repository Staff Only: edit this item