Edit Item 

PDF
agrid.pdf Download (860Kb)  Preview 

Cite this document as:

Two approaches pertaining to modeling largescale ocean circulation on unstructured meshes are described. Both use the finitevolume ideology, unstructured surface triangular mesh and geopotential vertical coordinate, and promise better numerical efficiency than P1–P1 finite element models. The first one is formulated on mediandual control volumes for all variables and presents a finitevolume implementation of P1–P1 finiteelement discretization (Agrid). The second one differs by the cellcentered placement of horizontal velocities (quasiBgrid). Two practical tasks have to be solved to ensure their stable performance in longterm simulations. For triangular Agrids, it is the stabilization against pressure modes triggered by the stepwise bottom topography. The proposed solution preserves volume and tracers by introducing a composite representation for the horizontal velocity (with an elementwise constant velocity correction). The quasiBgrid setup is free of pressure modes but requires efficient filtering and dissipation in the momentum equation because of its too large velocity space. Implementations of momentum advection and viscosity that serve this goal are proposed. Both setups show stable performance and similar numerical efficiency, as exemplified by simulations of a baroclinic channel flow and circulation in the North Atlantic.