Multi-resolution unstructured grid-generation for geophysical applications on the sphere

Darren Engwirda

DOI: https://doi.org/10.48550/arXiv.1512.00307

2015-12-01

Abstract:An algorithm for the generation of non-uniform unstructured grids on ellipsoidal geometries is described. This technique is designed to generate high quality triangular and polygonal meshes appropriate for general circulation modelling on the sphere, including applications to atmospheric and ocean simulation, and numerical weather predication. Using a recently developed Frontal-Delaunay-refinement technique, a method for the construction of high-quality unstructured ellipsoidal Delaunay triangulations is introduced. A dual polygonal grid, derived from the associated Voronoi diagram, is also optionally generated as a by-product. Compared to existing techniques, it is shown that the Frontal-Delaunay approach typically produces grids with near-optimal element quality and smooth grading characteristics, while imposing relatively low computational expense. Initial results are presented for a selection of uniform and non-uniform ellipsoidal grids appropriate for large-scale geophysical applications. The use of user-defined mesh-sizing functions to generate smoothly graded, non-uniform grids is discussed.

Atmospheric and Oceanic Physics,Computational Geometry,Numerical Analysis,Computational Physics,Fluid Dynamics

What problem does this paper attempt to address?

This paper aims to solve the problem of generating non - uniform and unstructured grids on spheres in geophysical applications. Specifically, the author proposes a new algorithm for generating high - quality triangular and polygonal grids on ellipsoidal geometries. These grids are particularly suitable for general circulation modeling in fields such as atmospheric and ocean simulations and numerical weather forecasting. By using the recently developed frontal Delaunay refinement technique, this method can construct high - quality unstructured ellipsoidal Delaunay triangulations and can generate dual polygonal grids derived from the associated Voronoi diagrams as a by - product. Compared with existing techniques, the frontal Delaunay method can generally produce near - optimal cell quality and smooth grading characteristics while maintaining a low computational cost. The paper also discusses how to use user - defined grid size functions to generate smoothly graded non - uniform grids.