Tracer transport within an unstructured grid ocean model using characteristic discontinuous Galerkin advection
Abstract
In a previous article a characteristic discontinuous Galerkin (CDG) advection scheme was presented for tracer transport (Lee et al., 2016). The scheme is conservative, unconditionally stable with respect to time step and scales sub-linearly with the number of tracers being advected. Here in this paper, we present the implementation of the CDG advection scheme for tracer transport within MPAS-Ocean, a Boussinesque unstructured grid ocean model with an arbitrary Lagrangian Eulerian vertical coordinate. The scheme is implemented in both the vertical and horizontal dimensions, and special care is taken to ensure that the scheme remains conservative in the context of moving vertical layers. Consistency is ensured with respect to the dynamics by a renormalization of the fluxes with respect to the volume fluxes derived from the continuity equation. For spherical implementations, the intersection of the flux swept regions and the Eulerian grid are determined for great circle arcs, and the fluxes and element assembly are performed on the plane via a length preserving projection. Lastly, solutions are presented for a suite of test cases and comparisons made to the existing flux corrected transport scheme in MPAS-Ocean.
- Authors:
- Publication Date:
- Research Org.:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Biological and Environmental Research (BER); USDOE National Nuclear Security Administration (NNSA)
- OSTI Identifier:
- 1989728
- Alternate Identifier(s):
- OSTI ID: 1477647; OSTI ID: 1636880
- Report Number(s):
- LA-UR-17-22608
Journal ID: ISSN 0898-1221; S0898122118305364; PII: S0898122118305364
- Grant/Contract Number:
- AC52-06NA25396; LA-UR-17-22608
- Resource Type:
- Published Article
- Journal Name:
- Computers and Mathematics with Applications (Oxford)
- Additional Journal Information:
- Journal Name: Computers and Mathematics with Applications (Oxford) Journal Volume: 78 Journal Issue: 2; Journal ID: ISSN 0898-1221
- Publisher:
- Elsevier
- Country of Publication:
- United Kingdom
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; Discontinuous Galerkin; Semi Lagrangian; Advection equation; Unstructured grid; Arbitrary Lagrangian Eulerian vertical coordinate
Citation Formats
Lee, D., Petersen, M., Lowrie, R., and Ringler, T. Tracer transport within an unstructured grid ocean model using characteristic discontinuous Galerkin advection. United Kingdom: N. p., 2019.
Web. doi:10.1016/j.camwa.2018.09.024.
Lee, D., Petersen, M., Lowrie, R., & Ringler, T. Tracer transport within an unstructured grid ocean model using characteristic discontinuous Galerkin advection. United Kingdom. https://doi.org/10.1016/j.camwa.2018.09.024
Lee, D., Petersen, M., Lowrie, R., and Ringler, T. Mon .
"Tracer transport within an unstructured grid ocean model using characteristic discontinuous Galerkin advection". United Kingdom. https://doi.org/10.1016/j.camwa.2018.09.024.
@article{osti_1989728,
title = {Tracer transport within an unstructured grid ocean model using characteristic discontinuous Galerkin advection},
author = {Lee, D. and Petersen, M. and Lowrie, R. and Ringler, T.},
abstractNote = {In a previous article a characteristic discontinuous Galerkin (CDG) advection scheme was presented for tracer transport (Lee et al., 2016). The scheme is conservative, unconditionally stable with respect to time step and scales sub-linearly with the number of tracers being advected. Here in this paper, we present the implementation of the CDG advection scheme for tracer transport within MPAS-Ocean, a Boussinesque unstructured grid ocean model with an arbitrary Lagrangian Eulerian vertical coordinate. The scheme is implemented in both the vertical and horizontal dimensions, and special care is taken to ensure that the scheme remains conservative in the context of moving vertical layers. Consistency is ensured with respect to the dynamics by a renormalization of the fluxes with respect to the volume fluxes derived from the continuity equation. For spherical implementations, the intersection of the flux swept regions and the Eulerian grid are determined for great circle arcs, and the fluxes and element assembly are performed on the plane via a length preserving projection. Lastly, solutions are presented for a suite of test cases and comparisons made to the existing flux corrected transport scheme in MPAS-Ocean.},
doi = {10.1016/j.camwa.2018.09.024},
journal = {Computers and Mathematics with Applications (Oxford)},
number = 2,
volume = 78,
place = {United Kingdom},
year = {Mon Jul 01 00:00:00 EDT 2019},
month = {Mon Jul 01 00:00:00 EDT 2019}
}
https://doi.org/10.1016/j.camwa.2018.09.024
Web of Science
Works referenced in this record:
Fully multidimensional flux-corrected transport algorithms for fluids
journal, June 1979
- Zalesak, Steven T.
- Journal of Computational Physics, Vol. 31, Issue 3
An Incremental Remapping Transport Scheme on a Spherical Geodesic Grid
journal, August 2005
- Lipscomb, William H.; Ringler, Todd D.
- Monthly Weather Review, Vol. 133, Issue 8
Incremental Remapping as a Transport/Advection Algorithm
journal, May 2000
- Dukowicz, John K.; Baumgardner, John R.
- Journal of Computational Physics, Vol. 160, Issue 1
A multi-resolution approach to global ocean modeling
journal, September 2013
- Ringler, Todd; Petersen, Mark; Higdon, Robert L.
- Ocean Modelling, Vol. 69
Anisotropic slope limiting for discontinuous Galerkin methods
journal, January 2017
- Aizinger, Vadym; Kosík, Adam; Kuzmin, Dmitri
- International Journal for Numerical Methods in Fluids, Vol. 84, Issue 9
An overview of research on Eulerian–Lagrangian localized adjoint methods (ELLAM)
journal, August 2002
- Russell, Thomas F.; Celia, Michael A.
- Advances in Water Resources, Vol. 25, Issue 8-12
An efficient WENO limiter for discontinuous Galerkin transport scheme on the cubed sphere
journal, September 2015
- Guo, Wei; Nair, Ramachandran D.; Zhong, Xinghui
- International Journal for Numerical Methods in Fluids, Vol. 81, Issue 1
Evaluation of the arbitrary Lagrangian–Eulerian vertical coordinate method in the MPAS-Ocean model
journal, February 2015
- Petersen, Mark R.; Jacobsen, Douglas W.; Ringler, Todd D.
- Ocean Modelling, Vol. 86
Spurious dianeutral mixing and the role of momentum closure
journal, January 2012
- Ilıcak, Mehmet; Adcroft, Alistair J.; Griffies, Stephen M.
- Ocean Modelling, Vol. 45-46
A simple weighted essentially nonoscillatory limiter for Runge–Kutta discontinuous Galerkin methods
journal, January 2013
- Zhong, Xinghui; Shu, Chi-Wang
- Journal of Computational Physics, Vol. 232, Issue 1
An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation
journal, December 1990
- Celia, Michael A.; Russell, Thomas F.; Herrera, Ismael
- Advances in Water Resources, Vol. 13, Issue 4
On Mass Conservation in High-Order High-Resolution Rigorous Remapping Schemes on the Sphere
journal, June 2013
- Erath, Christoph; Lauritzen, Peter H.; Tufo, Henry M.
- Monthly Weather Review, Vol. 141, Issue 6
A compatible and conservative spectral element method on unstructured grids
journal, August 2010
- Taylor, Mark A.; Fournier, Aimé
- Journal of Computational Physics, Vol. 229, Issue 17
A Conservative Semi-Lagrangian Discontinuous Galerkin Scheme on the Cubed Sphere
journal, January 2014
- Guo, Wei; Nair, Ramachandran D.; Qiu, Jing-Mei
- Monthly Weather Review, Vol. 142, Issue 1
A standard test case suite for two-dimensional linear transport on the sphere
journal, January 2012
- Lauritzen, P. H.; Skamarock, W. C.; Prather, M. J.
- Geoscientific Model Development, Vol. 5, Issue 3
A conservative semi-Lagrangian multi-tracer transport scheme (CSLAM) on the cubed-sphere grid
journal, March 2010
- Lauritzen, Peter H.; Nair, Ramachandran D.; Ullrich, Paul A.
- Journal of Computational Physics, Vol. 229, Issue 5
A vertex-based hierarchical slope limiter for <mml:math altimg="si27.gif" display="inline" overflow="scroll" xmlns:xocs="http://www.elsevier.com/xml/xocs/dtd" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.elsevier.com/xml/ja/dtd" xmlns:ja="http://www.elsevier.com/xml/ja/dtd" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:tb="http://www.elsevier.com/xml/common/table/dtd" xmlns:sb="http://www.elsevier.com/xml/common/struct-bib/dtd" xmlns:ce="http://www.elsevier.com/xml/common/dtd" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:cals="http://www.elsevier.com/xml/common/cals/dtd"><mml:mi>p</mml:mi></mml:math>-adaptive discontinuous Galerkin methods
journal, April 2010
- Kuzmin, Dmitri
- Journal of Computational and Applied Mathematics, Vol. 233, Issue 12
A high order characteristic discontinuous Galerkin scheme for advection on unstructured meshes
journal, November 2016
- Lee, D.; Lowrie, R.; Petersen, M.
- Journal of Computational Physics, Vol. 324
Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works
journal, January 1973
- Boris, Jay P.; Book, David L.
- Journal of Computational Physics, Vol. 11, Issue 1
A “Vertically Lagrangian” Finite-Volume Dynamical Core for Global Models
journal, October 2004
- Lin, Shian-Jiann
- Monthly Weather Review, Vol. 132, Issue 10