loading page

New constitutive equations for the teardrop and parabolic lens yield curves in viscous-plastic sea ice models
  • Damien Ringeisen,
  • Martin Losch,
  • Bruno Tremblay
Damien Ringeisen
McGill University

Corresponding Author:damien.ringeisen@mcgill.ca

Author Profile
Martin Losch
Alfred Wegener Institute for Polar and Marine Research
Author Profile
Bruno Tremblay
McGill University
Author Profile


Most viscous-plastic sea ice models use the elliptical yield curve. This yield curve has the fundamental flaw that it excludes acute angles between deformation features at high resolution. Conceptually, the teardrop and parabolic lens yield curves offer an attractive alternative. These yield curves feature a non-symmetrical shape, a Coulombic behavior for the low-medium compressive stress, and a continuous transition to the ridging-dominant mode. We show that the current formulation of the teardrop and parabolic lens viscous-plastic yield curves with normal flow rules results in negative or zero bulk and shear viscosities and consequently poor numerical convergence and representation of stress states on or within the yield curve. These issues are mainly linked to the assumption that the constitutive equation applicable for the elliptical yield curve also applies to non-symmetrical yield curves and yield curves with tensile strength. We present a new constitutive relation for the teardrop and parabolic lens yield curves that solves the numerical convergence issues naturally. Results from simple uni-axial loading experiments show that we can reduce the residual norm of the numerical solver with a smaller number of total solver iterations, resulting in significant improvements in numerical efficiency and representation of the stress and deformation field.
07 Dec 2022Submitted to ESS Open Archive
07 Dec 2022Published in ESS Open Archive