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.
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.
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.