An optimization approach to modeling sea ice dynamics. Part 1

Lagrangian framework

Helga S. Huntley, Esteban G. Tabak, Edward H. Suh

Research output: Contribution to journalArticle

Abstract

A new model for the dynamics of sea ice is proposed. The pressure field, instead of being derived from a local rheology as in most existing models, is computed from a global optimization problem. Here the pressure is seen as emerging not from an equation of state but as a Lagrange multiplier that enforces the ice's resistance to compression while allowing divergence. The resulting variational problem is solved by minimizing the pressure globally throughout the domain, constrained by the equations of momentum and mass conservation, as well as the limits on ice concentration (which has to stay between 0 and 1). This formulation has an attractive mathematical elegance while being physically motivated. Moreover, it leads to an analytic formulation that is also easily implemented in a numerical code, which exhibits marked stability and is suited to capturing discontinuities. In order to test the theory, the equations for a one-dimensional model are cast in terms of Lagrangian mass coordinates. The solution to the minimization problem is compared to an exact analytic solution derived using jump conditions in a simple test case. Another case is examined, which is somewhat more complicated but still allows our physical intuition to verify the qualitative results of the model. Good agreement is found. A final validation is performed by a comparison with a particle-based model, which tracks individual ice floes and their inelastic interaction in a one-dimensional domain.

Original languageEnglish (US)
Pages (from-to)543-560
Number of pages18
JournalSIAM Journal on Applied Mathematics
Volume67
Issue number2
DOIs
StatePublished - Dec 2006

Fingerprint

Sea Ice
Sea ice
Optimization
Ice
Modeling
Jump Conditions
Mass Conservation
Rheology
Formulation
One-dimensional Model
Lagrange multipliers
Variational Problem
Analytic Solution
Equation of State
Model
Global Optimization
Minimization Problem
Global optimization
Discontinuity
Divergence

Keywords

  • Ice dynamics
  • Lagrangian fluid dynamics
  • Rheology

ASJC Scopus subject areas

  • Mathematics(all)
  • Applied Mathematics

Cite this

An optimization approach to modeling sea ice dynamics. Part 1 : Lagrangian framework. / Huntley, Helga S.; Tabak, Esteban G.; Suh, Edward H.

In: SIAM Journal on Applied Mathematics, Vol. 67, No. 2, 12.2006, p. 543-560.

Research output: Contribution to journalArticle

@article{865fe76fc5db4227a5e1f832c9526f9e,
title = "An optimization approach to modeling sea ice dynamics. Part 1: Lagrangian framework",
abstract = "A new model for the dynamics of sea ice is proposed. The pressure field, instead of being derived from a local rheology as in most existing models, is computed from a global optimization problem. Here the pressure is seen as emerging not from an equation of state but as a Lagrange multiplier that enforces the ice's resistance to compression while allowing divergence. The resulting variational problem is solved by minimizing the pressure globally throughout the domain, constrained by the equations of momentum and mass conservation, as well as the limits on ice concentration (which has to stay between 0 and 1). This formulation has an attractive mathematical elegance while being physically motivated. Moreover, it leads to an analytic formulation that is also easily implemented in a numerical code, which exhibits marked stability and is suited to capturing discontinuities. In order to test the theory, the equations for a one-dimensional model are cast in terms of Lagrangian mass coordinates. The solution to the minimization problem is compared to an exact analytic solution derived using jump conditions in a simple test case. Another case is examined, which is somewhat more complicated but still allows our physical intuition to verify the qualitative results of the model. Good agreement is found. A final validation is performed by a comparison with a particle-based model, which tracks individual ice floes and their inelastic interaction in a one-dimensional domain.",
keywords = "Ice dynamics, Lagrangian fluid dynamics, Rheology",
author = "Huntley, {Helga S.} and Tabak, {Esteban G.} and Suh, {Edward H.}",
year = "2006",
month = "12",
doi = "10.1137/040621156",
language = "English (US)",
volume = "67",
pages = "543--560",
journal = "SIAM Journal on Applied Mathematics",
issn = "0036-1399",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "2",

}

TY - JOUR

T1 - An optimization approach to modeling sea ice dynamics. Part 1

T2 - Lagrangian framework

AU - Huntley, Helga S.

AU - Tabak, Esteban G.

AU - Suh, Edward H.

PY - 2006/12

Y1 - 2006/12

N2 - A new model for the dynamics of sea ice is proposed. The pressure field, instead of being derived from a local rheology as in most existing models, is computed from a global optimization problem. Here the pressure is seen as emerging not from an equation of state but as a Lagrange multiplier that enforces the ice's resistance to compression while allowing divergence. The resulting variational problem is solved by minimizing the pressure globally throughout the domain, constrained by the equations of momentum and mass conservation, as well as the limits on ice concentration (which has to stay between 0 and 1). This formulation has an attractive mathematical elegance while being physically motivated. Moreover, it leads to an analytic formulation that is also easily implemented in a numerical code, which exhibits marked stability and is suited to capturing discontinuities. In order to test the theory, the equations for a one-dimensional model are cast in terms of Lagrangian mass coordinates. The solution to the minimization problem is compared to an exact analytic solution derived using jump conditions in a simple test case. Another case is examined, which is somewhat more complicated but still allows our physical intuition to verify the qualitative results of the model. Good agreement is found. A final validation is performed by a comparison with a particle-based model, which tracks individual ice floes and their inelastic interaction in a one-dimensional domain.

AB - A new model for the dynamics of sea ice is proposed. The pressure field, instead of being derived from a local rheology as in most existing models, is computed from a global optimization problem. Here the pressure is seen as emerging not from an equation of state but as a Lagrange multiplier that enforces the ice's resistance to compression while allowing divergence. The resulting variational problem is solved by minimizing the pressure globally throughout the domain, constrained by the equations of momentum and mass conservation, as well as the limits on ice concentration (which has to stay between 0 and 1). This formulation has an attractive mathematical elegance while being physically motivated. Moreover, it leads to an analytic formulation that is also easily implemented in a numerical code, which exhibits marked stability and is suited to capturing discontinuities. In order to test the theory, the equations for a one-dimensional model are cast in terms of Lagrangian mass coordinates. The solution to the minimization problem is compared to an exact analytic solution derived using jump conditions in a simple test case. Another case is examined, which is somewhat more complicated but still allows our physical intuition to verify the qualitative results of the model. Good agreement is found. A final validation is performed by a comparison with a particle-based model, which tracks individual ice floes and their inelastic interaction in a one-dimensional domain.

KW - Ice dynamics

KW - Lagrangian fluid dynamics

KW - Rheology

UR - http://www.scopus.com/inward/record.url?scp=34247560620&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=34247560620&partnerID=8YFLogxK

U2 - 10.1137/040621156

DO - 10.1137/040621156

M3 - Article

VL - 67

SP - 543

EP - 560

JO - SIAM Journal on Applied Mathematics

JF - SIAM Journal on Applied Mathematics

SN - 0036-1399

IS - 2

ER -