Submitted as: development and technical paper 11 Aug 2021

Submitted as: development and technical paper | 11 Aug 2021

Review status: this preprint is currently under review for the journal GMD.

Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

Adrian K. Turner1, Kara J. Peterson2, and Dan Bolintineanu2 Adrian K. Turner et al.
  • 1Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
  • 2Sandia National Laboratories, Albuquerque, NM, USA

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-gridscale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds-preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of non-uniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

Adrian K. Turner et al.

Status: open (extended)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on gmd-2021-199', Anonymous Referee #1, 10 Sep 2021 reply

Adrian K. Turner et al.

Adrian K. Turner et al.


Total article views: 266 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
222 39 5 266 5 2
  • HTML: 222
  • PDF: 39
  • XML: 5
  • Total: 266
  • BibTeX: 5
  • EndNote: 2
Views and downloads (calculated since 11 Aug 2021)
Cumulative views and downloads (calculated since 11 Aug 2021)

Viewed (geographical distribution)

Total article views: 217 (including HTML, PDF, and XML) Thereof 217 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 21 Oct 2021
Short summary
We developed a technique to remap sea-ice tracer quantities between circular discrete element distributions. This is needed for a global discrete element method sea-ice model being developed jointly by Los Alamos National Laboratory and Sandia National Laboratories that has the potential to better utilize newer supercomputers with graphics processing units and better represent sea-ice dynamics. This new remapping technique ameliorates the effect of element distortion created by sea-ice ridging.