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.

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:

- dirichlet_condition.m, evaluates the right hand sides of the Dirichlet boundary conditions.
- rhs.m, evaluates the right hand sides (source terms) of the Stokes equations.
- channel_output.txt, output from a run of the program;

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.

- nodes3.txt, a text file containing the pressure nodes;
- nodes3.png, a PNG image of the pressure nodes;
- triangles3.txt, the linear triangulation;
- triangles3.png, a PNG image of the linear triangulation;
- pressure3.txt, a text file containing the pressure P at each linear node;
- pressure3.png, a PNG image of a contour plot of the pressure, produced by TRIANGULATION_ORDER3_CONTOUR.

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

- nodes6.txt, a text file containing the velocity nodes;
- nodes6.png, a PNG image of the velocity nodes;
- triangles6.txt, the quadratic triangulation;
- triangles6.png, a PNG image of the order 6 triangulation;
- velocity6.txt, a text file containing the velocity (U,V) at each velocity node;
- velocity6_dir.png, a PNG image of the velocity direction field, created by VECTOR_PLOT.
- velocity6_vec.png, a PNG image of the velocity field, created by VECTOR_PLOT.
- velocity6_arrows.png, a PNG image of the velocity field, created by VELOCITY_ARROWS.