alex-huth / gimpm_ssa_damage

A Generalized Interpolation Material Point Method for the Shallow Shelf Approximation of ice flow with Damage

GIMPM-SSA-Damage

The Generalized Interpolation Material Point Method (GIMPM) for the Shallow Shelf Approximation (SSA) of ice flow with Damage

Developer: Alex Huth ([email protected] ; [email protected])

This repository contains the GIMPM-SSA-Damage model detailed in:

Huth, A., Duddu, R., & Smith, B. (2021a). A generalized interpolation material point method for shallow ice shelves. 1: Shallow shelf approximation and ice thickness evolution. Journal of Advances in Modeling Earth Systems, 13, e2020MS002277. https://doi.org/10.1029/2020MS002277

Huth, A., Duddu, R., & Smith, B. (2021b). A generalized interpolation material point method for shallow ice shelves. 2: Anisotropic nonlocal damage mechanics and rift propagation. Journal of Advances in Modeling Earth Systems, 13, e2020MS002292. https://doi.org/10.1029/2020MS002292

In addition to the GIMPM, this code also includes the standard Material Point Method (sMPM)

Damage models included (see Huth et al., 2021b for details):

• SSA creep damage (Huth et al., 2021b)
• SSA "zero-stress" damage (Sun et al., 2017)
• The SSA "zero-stress" damage model + a modification to include necking/mass balance effects (Bassis & Ma, 2015)

Compilation

See README in PROG, which contains the main source code.

Test cases

The examples from Huth et al., 2021a are found in test1d, test2d, and mismip/steady The examples from Huth et al., 2021b are found in mismip/damage Each directory contains a README with instructions for running the examples.

Notes

• An installation and knowledge of Elmer FEM and Elmer/Ice is required
  • https://github.com/ElmerCSC/elmerfem
  • https://elmerfem.org/elmerice/wiki/
• Tested with Elmer v 9.0