A high order characteristic discontinuous Galerkin scheme for advection on unstructured meshes

A new characteristic discontinuous Galerkin (CDG) advection scheme is presented. In contrast to standard discontinuous Galerkin schemes, the test functions themselves follow characteristics in order to ensure conservation and the edges of each element are also traced backwards along characteristics...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2016-11, Vol.324 (C), p.289-302
Main Authors: Lee, D., Lowrie, R., Petersen, M., Ringler, T., Hecht, M.
Format: Article
Language:English
Subjects:
Advection equation
CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS
Discontinuous Galerkin
High order advection
Lagrangian characteristics
Mathematics
Unstructured grid
Unstructured mesh
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:A new characteristic discontinuous Galerkin (CDG) advection scheme is presented. In contrast to standard discontinuous Galerkin schemes, the test functions themselves follow characteristics in order to ensure conservation and the edges of each element are also traced backwards along characteristics in order to create a swept region, which is integrated in order to determine the mass flux across the edge. Both the accuracy and performance of the scheme are greatly improved by the use of large Courant–Friedrichs–Lewy numbers for a shear flow test case and the scheme is shown to scale sublinearly with the number of tracers being advected, outperforming a standard flux corrected transport scheme for 10 or more tracers with a linear basis. Moreover the CDG scheme may be run to arbitrarily high order spatial accuracy and on unstructured grids, and is shown to give the correct order of error convergence for piecewise linear and quadratic bases on regular quadrilateral and hexahedral planar grids. Using a modal Taylor series basis, the scheme may be made monotone while preserving conservation with the use of a standard slope limiter, although this reduces the formal accuracy of the scheme to first order. The second order scheme is roughly as accurate as the incremental remap scheme with nonlocal gradient reconstruction at half the horizontal resolution. The scheme is being developed for implementation within the Model for Prediction Across Scales (MPAS) Ocean model, an unstructured grid finite volume ocean model.
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2016.08.010