Multi-resolution unstructured grid-generation for geophysical applications on the sphere
It addresses the need for efficient, high-quality mesh generation for atmospheric and ocean simulation on the sphere, but the improvement is incremental over existing Delaunay-based methods.
This paper presents a Frontal-Delaunay-refinement algorithm for generating high-quality unstructured triangular and polygonal meshes on ellipsoidal geometries, achieving near-optimal element quality and smooth grading with low computational cost, suitable for geophysical modeling.
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.