A full‐stress, thermomechanical flow band model using the finite volume method

A thermomechanical ice flow model is formulated using the finite volume method. Separate submodels solve the full, two‐dimensional momentum equations, the advective‐diffusive heat equation, and evolution of the free surface. A unique aspect of the method is the use of a boundary‐fitted, orthogonal, curvilinear coordinate system, which simplifies the implementation of boundary conditions, leads to a straightforward discretization scheme, and results in banded sparse coefficient matrices that can be inverted directly. For simple boundary conditions and geometries, the model compares well with analytical solutions. For more complicated boundary conditions and geometries, the model compares well with full‐stress solutions obtained by previous authors.