Open AccessA full‐stress, thermomechanical flow band model using the finite volume method
Author(s) -
Price S. F.,
Waddington E. D.,
Conway H.
Publication year - 2007
Publication title -
journal of geophysical research: earth surface
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.67
H-Index - 298
eISSN - 2156-2202
pISSN - 0148-0227
DOI - 10.1029/2006jf000724
Subject(s) - curvilinear coordinates , finite volume method , discretization , boundary value problem , flow (mathematics) , coordinate system , mechanics , boundary (topology) , momentum (technical analysis) , mathematical analysis , mathematics , geometry , physics , finance , economics
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.