Special Session 140: 

Resolving large-scale geophysical flows over unstructured meshes

Qingshan Chen
Clemson University
USA
Co-Author(s):    Lili Ju and Roger Temam
Abstract:
Unstructured meshes have been gaining popularity in recent years, because they are almost free of polar singularities, and remain highly scalable even at eddy resolving resolutions. However, to unleash the full potential of these meshes, new schemes are needed. The classical C-grid scheme, which is widely popular on structured meshes, has serious issues concerning the reconstruction of the tangential velocity component. This talk presents new numerical schemes based on an old idea, namely the collocated vorticity-divergence formulation (so-called Z-grid), for large-scale geophysical flows on unstructured centroidal Voronoi meshes. Using the finite-volume discretization technique, the schemes conserve the mass and the absolute vorticity locally, and the potential enstrophy globally. It is also shown that, in an area-averaged sense, the schemes reproduce the Lagrangian transport property for potential vorticity, which is fundamental to the understanding of the dynamics of large-scale geophysical flows. A major challenge of vorticity-divergence based numerical schemes is the specification of the boundary conditions for the PDEs. This project adopts a hybrid approach that combines explicit and implicit implementation of the boundary conditions on the streamfunction and velocity potential. This talk will go over the analytical and practical aspects of the schemes, and finish with some high-resolution numerical results.