tcell_mass


tcell_mass, a Fortran77 code which computes and saves mass matrix associated with a finite element approximation of the fluid velocity in the time-dependent Navier Stokes equations in a region called the "T-Cell", by Hyung-Chun Lee.

The actual simulation of successive solutions to the flow problem can be done by the TCELL program, which automatically builds and uses the mass matrix. However, what TCELL_MASS does is allow us to "export" the information about the mass matrix (and several other important quantities) so that a reduced order model of the TCELL flow problem can be solved.

In fact, enough information about the reduced order model will be saved so that the resulting system can be solved by a simple program that invokes a MATLAB ODE solver.

The "T-Cell" region is a T-shaped region contained in a subset of the unit square [01]x[01]. The top half of the region is a channel open to the left and right, with a rectangular recess in the bottom. The left channel opening extends from [0.0,0.5] to [0.0,1.0], and the right channel from [1.0,0.5] to [1.0,1.0]. The recess is a square, whose bottom extends from [0.25,0.0] to [0.75,0.0], and which extends upward to meet the channel.

             +--------------------+
             +--------------------+
            ->                      ->
      Inlet --->      Channel       --->  Outlet
            ->                      ->
             +-----+          +---+
             +----+|  Recess  |+--+
                  ||          ||
             Void ||          || Void
                  |+----------+|
                  +------------+
      
The grid can be thought of as being generated by a uniform 81 by 81 uniformly spaced grid of nodes, from which the nodes corresponding to the two voids have been deleted. This leaves 4,961 nodes. These nodes are organized into 2400 6-node triangular elements which are quadratic in velocity and linear in pressure.

Usage:

tcell_mass
runs the program

Author:

Hyung-Chun Lee,
Ajou University, Korea,
in collaboration with Max Gunzburger and John Burkardt.

Licensing:

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

Languages:

tcell_mass is available in a Fortran77 version.

Related Data and Programs:

tcell_mass_test

BUMP, a Fortran90 program which solves a fluid flow problem in a channel including a bump which obstructs and redirects the flow.

CVT_BASIS_FLOW, a Fortran90 program which extracts dominant solution modes from a set of flow solutions using CVT methods.

POD_BASIS_FLOW, a Fortran90 program which extracts dominant solution modes from a set of flow solutions using POD methods.

TCELL, a Fortran77 program which solves the time-dependent Navier Stokes fluid flow equations in a 2D T-shaped region.

TCELL, a dataset directory which contains 500 successive flow solutions of the T-Cell problem.

Reference:

  1. Max Gunzburger,
    Finite Element Methods for Viscous Incompressible Flows,
    A Guide to Theory, Practice, and Algorithms,
    Academic Press, 1989.
  2. John Burkardt, Max Gunzburger, Hyung-Chun Lee,
    Centroidal Voronoi Tessellation-Based Reduced-Order Modelling of Complex Systems,
    SIAM Journal on Scientific Computing,
    Volume 28, Number 2, 2006, pages 459-484.

Source Code:


Last revised on 15 December 2023.