NAVIER_STOKES_MESH3D
3D Meshes for Navier Stokes Problems


NAVIER_STOKES_MESH3D is a set of MATLAB data files defining meshes for several 3D test problems involving the Navier Stokes equations for fluid flow, provided by Leo Rebholz.

We are interested in these meshes as useful tests for a procedure in which we are able to redo the related Navier Stokes calculations using FENICS.

As part of that task, we need to take the meshes associated with the MATLAB computation and translate them into a format acceptable to DOLFIN / FENICS. Our proposed mechanism for this is either to translate the files to the Gmsh msh format, and then use DOLFIN-CONVERT to convert that to DOLFIN XML, or else to translate the files directly into DOLFIN XML format ourselves.

Licensing:

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

Languages:

NAVIER_STOKES_MESH3D is available in a MATLAB version.

Related Data and Programs:

DOLFIN, Python programs which are used with the DOLFIN library to set up and solve a wide variety of problems using the finite element method.

DOLFIN-CONVERT, a Python program which can convert mesh file from Gmsh, MEDIT, METIS or SCOTCH format to an XML format suitable for use by DOLFIN or FENICS, by Anders Logg.

DOLFIN_XML is a data directory which contains examples of XML files that describe 3D finite element meshes as used by DOLFIN and FENICS.

FENICS, programs which illustrate the use of a collection of free software with an extensive list of features for automated, efficient solution of differential equations.

GMSH, examples which illustrate the use of the Gmsh program, a 3D mesh generator for the finite element method (FEM).

NAVIER_STOKES_3D_EXACT, evaluates an exact solution to the incompressible time-dependent Navier-Stokes equations over an arbitrary domain in 3D.

NAVIER_STOKES_MESH2D, MATLAB data files defining meshes for several 2D test problems involving the Navier Stokes equations for flow flow, provided by Leo Rebholz.

TET_MESH_ORDER20, a dataset directory which contains examples of order 20 tetrahedral meshes.

TET_MESH_TO_GMSH, a MATLAB program which reads information defining a tet mesh, namely a file of node coordinates and a file of elements defined by node indices, and creates a corresponding Gmsh mesh file.

Author:

Leo Rebholz

Reference:

  1. Abigail Bowers, Leo Rebholz, Aziz Takhirov, Catalin Trenchea,
    Improved accuracy in regularization models of incompressible flow via adaptive nonlinear filtering,
    International Journal for Numerical Methods in Fluids,
    Volume 70, Number 7, 2012, pages 805-828.

Source Code:

The following MAT files each contain a mesh, defined by a structure whose components are nodeco, elnode, and bdynde:

The following functions operate on the mesh data:

Examples and Tests:

CYLINDER_3D is a tet mesh for 3D flow in a 0.41x2.5x0.41 channel, past a (square) cylinder, using 49,602 nodes and 9,792 elements, created from cylinder_3d.mat.

STEP_3D is a tet mesh associated with 3D flow in a 10x50x10 channel, past a forward and backward step, using 38,992 nodes and 7,680 elements, created from step_3d.mat.

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


Last revised on 03 June 2013.