NAST2D_F90
A Finite Volume Code for Fluid Flow

NAST2D_F90 is a FORTRAN90 program which implements the finite volume method to solve for the transient velocity, pressure, and temperature of an incompressible fluid in a variety of 2D flow regions.

NAST2D_F90 is a FORTRAN90 version of the NAST2D program, which was originally written in C.

NAST2D_F90 writes graphics files, in particular, a text file containing the stream function at every gridpoint. It can also create "streakline" files, containing the position of marker particles over time.

Usage:

nast2d_f90 < input.par > output

reads the parameter file input.par is a "parameter file", that specifies the problem to be solved and associated parameter settings, and runs the program.

Languages:

NAST2D_F90 is available in a FORTRAN90 version.

Related Data and Programs:

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

DIRECTION_ARROWS_GRID is a MATLAB program which reads files of node and velocity data, and, using interpolation, creates a velocity direction plot with arrows place on a uniform grid of the user's specification.

FREE_FEM_NAVIER_STOKES is a FORTRAN90 program which solves the 2D incompressible Navier Stokes equations in an arbitrary triangulated region.

FREE_FEM_STOKES is a FORTRAN90 program that solves the steady Stokes flow equations on a triangulated 2D region.

HCELL is a FORTRAN77 program which computes the pressure and velocity for a Navier Stokes flow in an "H"-shaped region.

INOUT is a FORTRAN77 program which computes the pressure and velocity for a Navier Stokes flow in a square region with an inlet and an outlet.

MHD_CONTROL is a FORTRAN90 program which tries to control the evolution of an MHD system so that a particular state is achieved.

MHD_FLOW is a FORTRAN90 program for the evolution of an MHD system.

NAST2D_CONTOUR is a MATLAB program which can produce contour images of the stream function data computed by NAST2D_F90.

NAST2D_F90 is also available in a C++ version.

NAST2D_STREAK_STILL is a MATLAB program which will display a single time frame of the streakline data;

NAST2D_STREAK_MOVIE is a MATLAB program which creates an animation of the streakline data over time.

TCELL is a FORTRAN77 program which computes the pressure and velocity for a Navier Stokes flow in a "T"-shaped region.

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

VECTOR_STREAM_GRID is a MATLAB program which reads node and vector data from a file, computes an interpolatory function, evaluates on a uniform grid of points specified by the user, and displays a streamline plot of the vector field.

VELOCITY_ARROWS_GRID is a MATLAB program which reads files of node and velocity data, and, using interpolation, creates a vector plot with arrows place on a uniform grid of the user's specification.

Reference:

1. Michael Griebel, Thomas Dornseifer, Tilman Neunhoeffer,
   Numerische Simulation in der Stroemungsmechanik,
   Vieweg 1995.
2. Michael Griebel, Thomas Dornseifer, Tilman Neunhoeffer,
   Numerical Simulation in Fluid Dynamics,
   SIAM, 1998,
   This book can be ordered directly from the SIAM web site.
3. http://www.math.jmu.edu/~dpruett/software/
   The NAST2D_F90 web site.

Source Code:

• main.f90, the main program;
• boundary.f90 handles boundary conditions;
• defs.h, an include file that contains the definitions of some constants;
• extras.f90, contains some extra routines to write data to files for visualization by MATLAB;
• init.f90, reads the user input file;
• interfaces.h, an include file that contains the subroutine interface definitions;
• nrtype.f90, reads the user input file;
• surface.f90 handles free boundaries;
• uvp.f90, handles the assembly and solution of the equations that define the state variables;
• visual.f90, routines that write out files for visualization;
• nast2d_f90.csh, commands to compile, link and load the source code;

Examples and Tests:

BACKSTEP is the backward-facing step:

• backstep.par, the parameter file;
• backstep.txt, the output file;

CIRCLE is flow past a circular obstacle:

• circle.par, the parameter file;
• circle.txt, the output file;
• circle.psi, the PSI file, containing the stream function values, which can be displayed by the NAST2D_CONTOUR MATLAB routine;

CONVECT is buoyancy flow with heated side walls:

• convect.par, the parameter file;
• convect.txt, the output file;

DAM is a free boundary problem with a breaking dam:

• dam.par, the parameter file;
• dam.txt, the output file;

DCAVITY is the driven cavity:

• dcavity.par, the parameter file;
• dcavity.out, the output file;

DROP follows the splash of a drop of liquid:

• drop.par, the parameter file;
• drop.txt, the output file;

FLTRAP is a buoyancy flow with obstacles:

• fltrap.par, the parameter file;
• fltrap.txt, the output file;

MOLDING is a flow based on injection molding:

• molding.par, the parameter file;
• molding.txt, the output file;

RAYLEIGH is a buoyancy flow with heated upper and lower walls:

• rayleigh.par, the parameter file;
• rayleigh.txt, the output file;

WAVE is a free surface flow over a step:

• wave.par, the parameter file;
• wave.txt, the output file;

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