A Sample 2D Stokes Problem

FEM2D_STOKES_SPARSE_CHANNEL is a rectangular region that is 3 units wide and 1 unit high. It is used as a sample problem for the Stokes solver.


To run the problem directly, you only need the user-supplied routines in dirichlet_condition.m and rhs.m, the node data in nodes6.txt, and the element data in triangles6.txt.

You can run the program with the MATLAB command

        fem2d_stokes_sparse ( 'nodes6.txt', 'triangles6.txt' )


The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.


FEM2D_STOKES_SPARSE_CHANNEL is available in a MATLAB version.

Related Data and Programs:

FEM2D_STOKES_SPARSE, a MATLAB program which solves the steady (time independent) incompressible Stokes equations on an arbitrary triangulated region, using the finite element method and MATLAB's sparse facility.

Source Code:

Some of the files needed to run the problem include:

The pressures are associated with a subset of the full set of nodes, and with an order 3 triangulation of those nodes. Note that, in the order 3 triangulation, the nodes are renumbered, and do NOT inherit the labels used in the order 6 triangulation.

The velocities are associated with the full set of nodes, and with an order 6 triangulation of those nodes.

Last revised on 07 February 2006.