FEM2D_STOKES_CAVITY is a C++ library which defines the "cavity" problem for FEM2D_STOKES solver. The cavity is a square region that is 1 unit wide and 1 unit high. The tangential velocity is specified to be 1 along the top boundary, with a zero normal component. On all other parts of the boundary, the velocity is specified to be zero.
To run the problem directly, you only need the user-supplied routines in cavity.cpp, the node data in nodes6.txt, and the element data in triangles6.txt.
You compile and link the Stokes solver with cavity.cpp, using commands like:
g++ fem2d_stokes.cpp cavity.cpp
mv a.out cavity
and then run the program with the command
./cavity 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_CAVITY is available in a C++ version and a FORTRAN90 version and a MATLAB version.
FEM2D_STOKES, a C++ program which applies the finite element method to a Stokes flow problem, given user input data.
Some of the files needed to run the problem include:
The geometry is defined by sets of nodes and triangles. The velocities use the full set of nodes, and quadratic (6 node) triangles.
The pressures are associated with a subset of the nodes called "pressure nodes", and linear (3 node) triangles. Note that, in the order 3 triangulation, the nodes are renumbered, and do NOT inherit the labels used in the order 6 triangulation.
The pressures are a scalar quantity associated with the pressure nodes, the velocities are a vector quantity associated with the vector nodes.
