A 1-D elastic-plastic sea-ice model solved with an implicit Eulerian-Lagrangian method

Research output: Contribution to journalArticlepeer-review


A physical model for an elastic-plastic rheology is developed and implemented in a numerical sea-ice model. The rheology describes sea ice as behaving as an elastic material for relatively small deformations and as a plastic material for larger ones. The model equations are solved using an Eulerian-Lagrangian method in which the displacement of granular aggregates from an original Eulerian grid is computed in a Lagrangian sense and the resulting mass distribution is mapped back onto the Eulerian grid. The equations are integrated in time using Krylov solvers in a fully-implicit framework. The model distinguishes itself from previous sea-ice models in a combination of attributes: the absence of a viscous dependence within the rheology, the employment of an Eulerian-Lagrangian grid, and in the use of a fully-implicit time-stepping scheme that allows for a large time-step. Model results and efficiency are presented from a one-dimensional simulation; plans for extension of the physics and numerics to two-dimensions are outlined.

Original languageEnglish (US)
Pages (from-to)1-27
Number of pages27
JournalOcean Modelling
Issue number1
StatePublished - 2007


  • Elastic-plastic
  • Eulerian-Lagrangian
  • Sea-ice dynamics

ASJC Scopus subject areas

  • Computer Science (miscellaneous)
  • Oceanography
  • Geotechnical Engineering and Engineering Geology
  • Atmospheric Science


Dive into the research topics of 'A 1-D elastic-plastic sea-ice model solved with an implicit Eulerian-Lagrangian method'. Together they form a unique fingerprint.

Cite this