FEM2D_STOKES_CHANNEL
"Channel" Routines for FEM2D_STOKES

FEM2D_STOKES_CHANNEL is a FORTRAN90 library which supplies information defining a Stokes flow problem in a channel. The channel is a rectangular region that is 3 units wide and 1 unit high.

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

Usage:

You compile and link the solver with cavity.f90, using commands like:

        gfortran fem2d_stokes.f90 channel.f90
        mv a.out channel
      

        ./channel nodes6.txt triangles6.txt
      

Licensing:

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

Languages:

FEM2D_STOKES_CHANNEL is available in a FORTRAN90 version.

Related Data and Programs:

FEM2D_STOKES, a FORTRAN90 program which solves the 2D incompressible Stokes equations in an arbitrary triangulated region. In order to run, it requires user-supplied routines that define problem data.

FEM2D_STOKES_CAVITY, a FORTRAN90 library which contains the user-supplied routines necessary to run fem2d_stokes on the "cavity" problem.

FEM2D_STOKES_INOUT, a FORTRAN90 library which contains the user-supplied routines necessary to run fem2d_stokes on the "inout" problem.

Source Code:

• channel.f90, the user-supplied routines to evaluate the right hand side and boundary conditions;
• channel.txt, output from a run of the program;

The geometry is defined by sets of nodes and triangles. The velocities use the full set of nodes, and quadratic (6 node) triangles.

• 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;

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.

• 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;

The pressures are a scalar quantity associated with the pressure nodes, the velocities are a vector quantity associated with the vector nodes.

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

List of Routines:

• BOUNDARY_TYPE determines the type of boundary conditions imposed.
• DIRICHLET_CONDITION sets the value of a Dirichlet boundary condition.
• NEUMANN_CONDITION sets the value of any Neumann boundary conditions.
• RHS gives the right-hand side of the differential equation.

You can go up one level to the FORTRAN90 source codes.