NS3D_FEM
Steady Navier Stokes Equations in 3D
Finite Element Solution


NS3D_FEM is a MATLAB program which sets up and solves a finite element formulation of the steady 3D Navier Stokes equations in an arbitrary geometry defined by the user, by Jeff Borggaard.

The user defines the geometry as a collection of tetrahedrons. These tetrahedrons may be linear (using four nodes) or quadratic (using ten nodes, which are the four vertices and the six midside nodes).

Graphics Options:

The data from the program can be written out to a graphics file suitable for use with a visualization program. This is done by "turning on" the appropriate call in steady_navierstokes_3d.m There are currently several functions that can be called, to make an output file of various types:

Information about TECPLOT is available at http://www.tecplot.com/.

Information about PARAVIEW is available at http://www.paraview.org/ .

Information about VU is available at http://www.invisu.ca/e/index.php.

Licensing:

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

Languages:

NS3D_FEM is available in a MATLAB version.

Related Data and Programs:

PARAVIEW, examples which illustrate the use of the paraview graphics program.

TEC, a data directory which contains examples of TEC files, TECPLOT files used for storing 2D or 3D graphics with attributes;

TET_MESH, a MATLAB library which carries out various operations on tetrahedral meshes.

TET_MESH_TET_NEIGHBORS, a MATLAB program which reads information about a tet mesh and writes out a list of the tetrahedrons that are adjacent to a given tetrahedron;

THREED, a MATLAB library which contains functions useful for 3D finite element calculations.

THREED_TO_TECPLOT, a MATLAB program which accepts computational data from a 3D finite element fluid flow program and writes it to a PLT file for postprocessing by TECPLOT.

THREED_TO_VTK, a MATLAB program which accepts computational data from a 3D finite element fluid flow program and writes it to a VTK "legacy" file for postprocessing by PARAVIEW.

THREED_TO_VTU, a MATLAB program which accepts computational data from a 3D finite element fluid flow program and writes it to a VTK unstructured grid file for postprocessing by PARAVIEW.

THREED_TO_VU, a MATLAB program which accepts computational data from a 3D finite element fluid flow program and writes it to a file for postprocessing by VU.

TOMS866, a MATLAB library which is the Incompressible Flow Iterative Solution Software;
this library is commonly called IFISS;
this is ACM TOMS algorithm 866.

VTK, a data directory which contains examples of VTK files, an obsolete file format used by the Visual Toolkit, and which can be displayed by the PARAVIEW program;

VTU, a data directory which contains examples of VTU files, a modern XML-based file format for unstructured grid data, used by the Visual Toolkit, and which can be displayed by the PARAVIEW program;

Author:

NS3D_FEM was written by Jeff Borggaard, Mathematics Department, Virginia Tech.

Reference:

  1. Max Gunzburger,
    Finite Element Methods for Viscous Incompressible Flows,
    A Guide to Theory, Practice, and Algorithms,
    Academic Press, 1989,
    ISBN: 0-12-307350-2,
    LC: TA357.G86.

Source Code:

Examples and Tests:

CASE1 is a built-in example, which can be run by the command

        steady_navierstokes_3d ( 1 )
      
All data for this example is generated internally; no input files are needed.

CASE2 is a built-in example, which can be run by the command

        steady_navierstokes_3d ( 2 )
      
This example needs some large external input files.

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


Last modified on 21 December 2010.