Steady 2D Stokes Model of a Glacier

FENICS_GLACIER is a FENICS example which defines the steady 2D Stokes equations for a model of the motion of ice in a glacier, by William Mitchell.

Since this example was published, FENICS has changed, and several commands used in the example no longer work. We are actively trying to figure out how to modify this example so that it works with the current version of FENICS.

Related Data and Programs:

DOLFIN-CONVERT, a Python program which can convert a mesh file from a given format to an XML format suitable for use by DOLFIN or FENICS, by Anders Logg.

DOLFIN_XML, a data directory which contains examples of XML files that describe 3D finite element meshes as used by DOLFIN and FENICS.

FENICS, examples which illustrate the use of FENICS, a collection of free software with an extensive list of features for automated, efficient solution of partial differential equations (PDE's).

FENICS_TO_FEM, examples which illustrate how a mesh or scalar function computed by the FENICS program can be written to FEM files, which can then be used to create images, or as input to meshing programs or other analysis tools.

MITCHELL_FENICS, examples which illustrate the implementation of the Mitchell elliptic test problems using FENICS.


  1. Anders Logg, Kent-Andre Mardal, Garth Wells,
    Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book,
    Lecture Notes in Computational Science and Engineering,
    Springer, 2011,
    ISBN13: 978-3642230981.
  2. Anders Logg, Garth Wells,
    DOLFIN: Automated finite element computing,
    ACM Transactions on Mathematical Software,
    Volume 37, Number 2, Article 20, April 2010.
  3. Martin Sandve Alnaes, Johan Hake, Robert Kirby, Hans Petter Langtangen, Anders Logg, Garth Wells,
    The FEniCS Manual,
  4. Hans Petter Langtangen,
    A FEniCS Tutorial,
  5. William Mitchell,
    Exact and numerical solutions for Stokes flows in glaciers,
    Master's thesis,
    Department of Mathematics and Statistics,
    University of Alaska, Fairbanks, Alasaka,
    August 2012.

Examples and Tests:

GLACIER solves the steady 2D Stokes equation that represent the motion of ice in a glacer,

You can go up one level to the Examples directory.

Last revised on 11 February 2015.