A discontinuous finite element baroclinic marine model on unstructured prismatic meshes

Abstract

We describe the space discretization of a three-dimensional baroclinic finite element model, based upon a discontinuous Galerkin method, while the companion paper (Comblen et al. 2010a) describes the discretization in time. We solve the hy- drostatic Boussinesq equations governing marine flows on a mesh made up of triangles extruded from the surface toward the seabed to obtain prismatic threedimensional elements. Diffusion is implemented using the symmetric interior penalty method. The tracer equation is consistent with the continuity equation. A Lax-Friedrichs flux is used to take into account internal wave propagation. By way of illustration, a flow exhibiting internal waves in the lee of an isolated seamount on the sphere is simulated. This enables us to show the advantages of using an unstructured mesh, where the resolution is higher in areas where the flow varies rapidly in space, the mesh being coarser far from the region of interest. The solution exhibits the expected wave structure. Linear and quadratic shape functions are used, and the extension to higher-order discretization is straightforward.

Improving the Upper‐Ocean Temperature in an Ocean Climate Model (FESOM 1.4): Shortwave Penetration Versus Mixing Induced by Nonbreaking Surface Waves

“…To adequately resolve these features and key regions, mesh resolutions are required to be much higher than that are currently affordable in long‐term simulations of global climate models. This promoted the development of unstructured‐mesh ocean models in the past decades (e.g., Blaise et al, ; C. Chen et al, ; Danilov et al, ; Ford et al, ; Q. Wang et al, ; White et al, ). The Finite Element Sea ice‐Ocean Model (FESOM), unlike most of the unstructured‐mesh ocean models that are intended for coastal and regional applications, is the first mature global unstructured‐mesh ocean model that was developed mainly for the purpose of climate research (Danilov et al, ; Sidorenko et al, ; Timmermann et al, ; Q. Wang et al, ).…”

Assessment of numerical schemes for solving the advection–diffusion equation on unstructured grids: case study of the Guaíba River, Brazil

“…Theoretically, the pressure gradient force can be divided into barotropic and baroclinic components by using the Leibniz integration rule (Blaise et al, 2010). And, since the baroclinic components are driven by vertical gradients of density and/or temperature over the water column, IPH-UnTRIM2D includes only barotropic components of pressure in the calculations as part of its depth-integrated approach (i.e.…”

Discontinuous Galerkin unsteady discrete adjoint method for real-time efficient tsunami simulations

“…It is based on the discontinuous Galerkin method, whose potential has already been demonstrated for marine modeling [2][3][4][5] or atmospheric flows [6][7][8]. The main advantages of this method, in addition to the use of unstructured grids, consist of good stability properties thanks to the finite volume fluxes, efficient parallel scaling due the high locality of the method and high-order accuracy.…”