Abstract
A third-order numerical scheme was developed for 2D irregular hexagonal meshes for the advection problems in this study. The scheme is based on a multi-moment constrained finite-volume method (MCV) in Cartesian coordinates and entails the introduction of a general integration method over a hexagonal cell. Unlike in the conventional finite-volume method, various discrete moments, that is, point value and volume-integrated average, are adopted as computational constraints to achieve high-order computation. The high-order spatial reconstruction can therefore be built in a local space, which considerably reduces the stencil length. The numerical scheme is tested using various idealized experiments. Compared with the existing schemes, this scheme is demonstrated to be flexible for application in irregular hexagonal meshes without increasing cost or compromising on accuracy. The general integration formulation based on a third-order polynomial helps to expand the application to arbitrary hexagons that does not require the use of centroids as computational points or Voronoi tessellation. It is also convenient to define the orthogonal wind components in the Cartesian system to directly drive the atmospheric transport.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have