Communication-Efficient Property Preservation in Tracer Transport

TitleCommunication-Efficient Property Preservation in Tracer Transport
Publication TypeJournal Article
Year of Publication2019
JournalSIAM Journal on Scientific Computing
Volume41
Number3
PagesC161-C193
Date Published06/2019
Abstract / Summary

Atmospheric tracer transport is a computationally demanding component of the atmospheric dynamical core of weather and climate simulations. Simulations typically have tens to hundreds of tracers. A tracer field is required to preserve several properties, including mass, shape, and tracer consistency. To improve computational efficiency, it is common to apply different spatial and temporal discretizations to the tracer transport equations than to the dynamical equations. Using different discretizations increases the difficulty of preserving properties. This paper provides a unified framework to analyze the property preservation problem and classes of algorithms to solve it. We examine the primary problem and a safety problem; describe three classes of algorithms to solve these; introduce new algorithms in two of these classes; make connections among the algorithms; analyze each algorithm in terms of correctness, bound on its solution magnitude, and its communication efficiency; and study numerical results. A new algorithm, QLT, has the smallest communication volume, and in an important case it redistributes mass approximately locally. These algorithms are only very loosely coupled to the underlying discretizations of the dynamical and tracer transport equations and thus are broadly and efficiently applicable. In addition, they may be applied to remap problems in applications other than tracer transport.

URLhttp://dx.doi.org/10.1137/18m1165414
DOI10.1137/18m1165414
Journal: SIAM Journal on Scientific Computing
Year of Publication: 2019
Volume: 41
Number: 3
Pages: C161-C193
Date Published: 06/2019

Atmospheric tracer transport is a computationally demanding component of the atmospheric dynamical core of weather and climate simulations. Simulations typically have tens to hundreds of tracers. A tracer field is required to preserve several properties, including mass, shape, and tracer consistency. To improve computational efficiency, it is common to apply different spatial and temporal discretizations to the tracer transport equations than to the dynamical equations. Using different discretizations increases the difficulty of preserving properties. This paper provides a unified framework to analyze the property preservation problem and classes of algorithms to solve it. We examine the primary problem and a safety problem; describe three classes of algorithms to solve these; introduce new algorithms in two of these classes; make connections among the algorithms; analyze each algorithm in terms of correctness, bound on its solution magnitude, and its communication efficiency; and study numerical results. A new algorithm, QLT, has the smallest communication volume, and in an important case it redistributes mass approximately locally. These algorithms are only very loosely coupled to the underlying discretizations of the dynamical and tracer transport equations and thus are broadly and efficiently applicable. In addition, they may be applied to remap problems in applications other than tracer transport.

DOI: 10.1137/18m1165414
Citation:
Bradley, A, P Bosler, O Guba, M Taylor, and G Barnett.  2019.  "Communication-Efficient Property Preservation in Tracer Transport."  SIAM Journal on Scientific Computing 41(3): C161-C193.  https://doi.org/10.1137/18m1165414.