TCELL
Generate Navier Stokes solutions in TCELL region.

TCELL is a FORTRAN77 program which solves the time-dependent Navier Stokes equations in a two dimensional region called the "T Cell", by Hyung-Chun Lee.

As time increases, other flow parameters may be varied as well. For this code, the flow parameter that can vary is the magnitude of the inflow specified at the left boundary. The program stores the results of each timestep in a file. Each solution may be thought of as a "snapshot" of the gradual evolution of the flow field over time. In a typical run, 500 such files might be created.

A collection of snapshot files may be analyzed to seek a small underlying basis of fundamental flows. This basis may be used to approximate the entire set, or indeed to generate a reduced basis program that can quickly and cheaply compute new flow solutions whose flow parameters vary in small ways from those of the original set.

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
                  |+----------+|
                  +------------+
      

Usage:

tcell

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 GNU LGPL license.

Languages:

TCELL is available in a FORTRAN77 version and a MATLAB version.

Related Data and Programs:

BUMP, a FORTRAN90 program which seeks the appropriate channel inflow and bump shape which will cause the flow to most closely match a given downflow profile.

CHANNEL, a FORTRAN77 program which seeks the appropriate channel inflow which will cause the flow to most closely match a given downflow profile.

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

HCELL, a FORTRAN77 program which solves the Navier Stokes equations in the H-Cell region.

INOUT, a FORTRAN77 program which solves the Navier Stokes equations in the "In-Out" region.

MHD_CONTROL, a FORTRAN90 program which controls the evolution of an MHD system so that a particular state is achieved.

MHD_FLOW, a FORTRAN90 program which simulates the evolution of a system governed by the equations for magnetohydrodynamics.

NAST2D, a C++ program which solves the 2D Navier Stokes fluid flow equations

NAST2D_F90, a FORTRAN90 program which uses the finite volume method to set up and solve the 2D incompressible Navier Stokes equations with heat.

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

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

TCELL_FLOW_MOVIE, an animation of a time simulation of the flow computed by TCELL.

TCELL_MASS, a FORTRAN77 program which computes the mass matrix for the T-Cell problem.

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

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

• tcell.f, the source code.
• fuini.txt, an input file of initial solution coefficient values.

Examples and Tests:

RUN1 is a run of the program used to generated "snapshots", which are flow solutions at equally spaced time intervals. The program expects an input file, containing coefficients that define the solution at the initial time. During this run, the boundary conditions are varied in a way that is intended to make the system exhibit many varied behaviors. (Right now, to change the boundary conditions requires altering the source code of the program.)

• fuini_run1.txt, an input file of initial solution coefficient values, for the run that generates snapshots.
• tcell_run1_output.txt, the output file.
• elenode.txt, an output file containing the number of nodes and elements, the node coordinates, and the node indices that form each element.
• The 500 snapshot files created by this run are available in the TCELL snapshot dataset.

RUN2 is a run of the program used to generate a solution for comparison with one generated by a reduced order model. The program expects an input file, containing coefficients that define the solution at the initial time. During this run, the boundary conditions are varied in a new way. (Right now, to change the boundary conditions requires altering the source code of the program.)

• fuini_run2.txt, an input file of initial solution coefficient values, for the run that generates snapshots.

List of Routines:

• MAIN is the main program for TCELL.
• CH_IS_DIGIT returns TRUE if a character is a decimal digit.
• CH_TO_DIGIT returns the integer value of a base 10 digit.
• DAXPY adds a multiple of one vector to another.
• DDOT forms the dot product of two vectors.
• DGBFA factors a double precision band matrix by elimination.
• DGBSL solves a double precision band linear system.
• DIGIT_INC increments a decimal digit.
• DIGIT_TO_CH returns the character representation of a decimal digit.
• DSCAL scales a vector by a constant.
• FILE_NAME_INC generates the next filename in a series.
• GRID_WRITE writes the node and element information to a file.
• IDAMAX finds the vector element of largest magnitude.
• NSTOKE solves the Navier-Stokes equations using Taylor-Hood elements.
• REFBSP evaluates a linear basis function on the reference triangle.
• REFQBF evaluates a quadratic basis function on the reference triangle.
• SETGRD sets up the grid for the problem.
• SOLUTION_INIT returns an initial value for the solution coefficient vector.
• TRANS transforms data between the reference and physical elements.
• UBDRY evaluates the boundary conditions.

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