INOUT_FLOW2 INOUT Flow Problem Solution Datasets

INOUT_FLOW2 is a dataset directory which contains velocity solutions to a time-dependent incompressible parameterized flow problem in in a cell with inflow and outflow.

The "INOUT" region is contained in the unit square, [0,1] x [0,1]. There is an inlet on the left side, from [0,0] to [0,0.2], and an outlet on the right side, from [1,0.8] to [1.0,1.0].

+-----------------------+
|+----------------------+
||                    --->
||                    ----->
||                    --->
||                     ++
||                     ||
||                     ||
||                     ||
||                     ||
++                     ||
--->                    ||
----->                  ||
--->                    ||
+----------------------+|
+-----------------------+

There are 3,721 nodes in a 61 by 61 evenly spaced grid. The nodes are organized into 1800 6-node triangular elements which are quadratic in velocity and linear in pressure.

The time-dependent incompressible 2D Navier-Stokes equations are to be applied in the region. No-slip and no-penetration conditions are applied along the walls. The inlet on the lower left is open, and parabolic velocity inflow is specified, whose strength is controlled by a parameter alpha. At the outlet on the upper right, the flow is specified to have zero vertical component, and zero horizontal derivative. The dynamic viscosity is taken to be nu=1/300.

The steady state solution is computed for the particular value alpha=1/3. Using the steady state solution as initial condition, the time-dependent problem is solved over the interval 0<=t<=40, using 4,000 timesteps of size dt=0.01. Only every 5th timestep is saved as a "snapshot". The parameter alpha is set to 5/3 for timesteps 1 through 400, and then alpha = 1/3 for timesteps 401 through 800. The two instantaneous changes to alpha induce violent transient behaviors.

The goal is to extract typical modes of behavior of the solution. Such a set of modes may then be used as a finite element basis that is highly tuned to the physics of the problem, so that a very small set of basis functions can be used to closely approximate the behavior of the solution over a range of values of alpha.

Related Data and Programs:

CASE1_FLOW, a dataset directory which 401 solutions of a flow problem in a channel;

CAVITY_FLOW, a dataset directory which contains 500 time steps of Navier-Stokes flow in a driven cavity;

CVT_BASIS_FLOW, a FORTRAN90 program which extracts dominant solution modes from a set of flow solutions using Centroidal Voronoi Tessellation (CVT) methods.

INOUT_FLOW, a dataset directory which contains 500 time steps of Navier-Stokes flow in a region with specified inflow and outflow;

INOUT_FLOW_DISPLAY, a MATLAB library which displays a single velocity field solution for the INOUT flow;

INOUT_FLOW_MOVIE, a MATLAB library which creates an animation of the velocity solutions for the INOUT cell;

POD_BASIS_FLOW, a FORTRAN90 program which extracts dominant solution modes from a set of flow solutions using Proper Orthogonal Decomposition (POD).

TCELL_FLOW, a dataset directory which contains 500 time steps of Navier-Stokes flow in a 2D T-shaped region.

TRIANGULATION_PLOT, a FORTRAN90 program which can be used to draw an image of the nodes and elements used in the triangulation of the region.

VECTOR_PLOT, a FORTRAN90 program which can be used to draw an image of the velocity vector field over the region.

Reference:

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

Example Datasets:

Grid files you may copy include:

Steady state velocity files you may copy include:

The data files are available in ZIP files:

• snap1.zip contains up001.dat through up200.dat.
• snap2.zip contains up201.dat through up400.dat.
• snap3.zip contains up401.dat through up600.dat.
• snap4.zip contains up601.dat through up800.dat.
• xy.zip, contains element.dat, steady.dat and xy.dat.

You can go up one level to the DATASETS directory.

Last revised on October 2005.