Ice-Sheet / Ocean Interaction Model for Greenland Fjords Using High-Order Discontinuous Galerkin Methods (Final Report)
This report summarizes the development and validation of the Non-hydrostatic Unified Model of the Ocean (NUMO). The overarching science goal of the project was to advance Earth System science by proposing a more realistic representation of physical processes in a Greenland fjord. To achieve this objective, we proposed to build a new, separate, high-resolution ocean model to realistically represent the complex three-dimensional physics of the fjord, its bathymetry, coastlines using continuous/discontinuous Galerkin (CG/DG) methods. The report describes the technical details of model components, i.e., choice of equations set, time integration scheme, boundary conditions, spatial discretization methods, mesh generation. We show a suite of validation results for test cases ranging from the Kovasznay flow, through the lock-exchange, simplified two-dimensional fjord, to a three-dimensional density current on a slope. Results are compared with the literature. The report is concluded with a discussion of efforts taken to perform a simulation of Sermilik Fjord with realistic bathymetry data.