A high order characteristic discontinuous Galerkin scheme for advection on unstructured meshes
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
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. Furthermore, the scheme is being developed for implementation within the Model for Prediction Across Scales (MPAS) Ocean model, an unstructured grid finite volume ocean model.
- Research Organization:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Biological and Environmental Research (BER)
- Grant/Contract Number:
- AC52-06NA25396; LA-UR-16-22694
- OSTI ID:
- 1460624
- Alternate ID(s):
- OSTI ID: 1359304
- Report Number(s):
- LA-UR-16-22694; TRN: US1901880
- Journal Information:
- Journal of Computational Physics, Vol. 324, Issue C; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
A Flux-Form Semi-Lagrangian advection scheme for tracer transport on arbitrary meshes
An Arbitrary Lagrangian-Eulerian Discretization of MHD on 3D Unstructured Grids